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Hydromagnetic waves, which were originally 
predicted theoretically by Alfvén,’ have been 
demonstrated in mercury or in liquid sodium in 
the laboratory.? A number of authors have dis- 
cussed the possible astrophysical® and geophysi- 
cal* importance of hydromagnetic waves. 

This Letter presents an observational result 
that may be considered as evidence of hydromag- 
netic waves propagated to the earth from its outer 
atmosphere. 

Dungey® and subsequently other authors® dis- 
cussed oscillations of the lines of magnetic force 
in the earth’s outer atmosphere, in order to ac- 
count for certain types of magnetic pulsations. 
For the sake of simplicity these authors assumed 
an isotropic conductivity. Piddington’ showed 
that because of the anisotropic conductivity of 
an ionized gas in the presence of a magnetic field, 
plane hydromagnetic waves of small amplitude, 
propagating along the magnetic field, are cir- 











FIG. 1. An example of damped hydromagnetic waves 
observed at College, Alaska, in the three orthogonal 
magnetic components, declination D, the horizontal 
force H, and the vertical force Z; 1 y =10~ gauss. 
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cularly polarized in opposite senses. For an 
oblique propagation the polarizations of the cor- 
responding waves become elliptical. 

Damped oscillations observed with the rapid- 
run magnetometers by the U. S. Coast and Geo- 
detic Survey Magnetic Observatory at College, 
Alaska, are now being analyzed. Figure 1 illus- 
trates a typical example. The curves, for the 
horizontal force H, declination D, and the verti- 
cal force Z, were traced from the rapid-run mag- 
netograms without smoothing. 

By combining H and D, changes in the hori- 
zontal magnetic vector were obtained at intervals 
of one minute. Figure 2 shows the locus of the 

















FIG. 2. The locus of the end point of the magnetic 


vector for the hydromagnetic waves shown in Fig. 1; 
projection onto the horizontal plane. 
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end point of the changing vector; circles indicate 
the epochs at which readings were actually made, 
and the attached numbers show the time elapsed, 
reckoned in minutes from the epoch 0459 U.T., 
April 4, 1960. 

The maximum ranges in D, H, and Z were 124, 
97, and 44 y (1 y=10°* gauss), respectively. The 
period varied to some extent from one cycle to 
the next, probably due to modifications of wave 
forms by various secondary effects; the average 
period was about 6 minutes. This is the typical 
period of damped oscillations, observed at Col- 
lege, of the kind discussed here. 

It is evident that the wave is very nearly ellip- 
tically polarized. The remarkably regular polari- 
zation of the magnetic waves strongly suggests 
that they are transverse hydromagnetic waves, 
generated in the outer atmosphere at an altitude 
of several earth radii, and propagated along the 
lines of magnetic force. This argument is sup- 
ported by two important related facts, to be 
demonstrated later: (a) Magnetic oscillations of 
this type are found to be mostly confined to the 
auroral zones. (b) They occur simultaneously 
in magnetically conjugate regions in the north- 
ern and southern hemispheres, namely at Col- 
lege (geomagnetic latitude 64°.7 N) and Macquarie 
Island (geomagnetic latitude 61°.1 S). The IGY 
records of the standard magnetometers at Mac- 
quarie Island have been compared with both the 
rapid-run and sensitive magnetograms taken at 
College; rapid-run magnetometers were not op- 
erated at Macquarie Island. The polarization at 
Macquarie Island is being studied, but the ac- 
curacy is limited owing to the slow speed of the 
standard magnetometers. 

The hydromagnetic waves are probably gener- 
ated by a sudden motion of the gas or by an elec- 
tric field suddenly imposed on the gas. There 
appear to be preferential local times at which 
the hydromagnetic waves are generated; this 
point will be discussed elsewhere. 

There is a possibility that these oscillations 
originate in the ionosphere, e. g., from transient 
phenomena such as those discussed by Cowling® 
in connection with the motion of solar material 
across magnetic fields, and by Piddington® in a 
more general form. Such oscillations might be 
propagated along the lines of magnetic force to 
the other hemisphere as in very low frequency 
radio noises of natural origin,’® known as whis- 
tlers, or in those generated artificially.“ It 
should be noted, however, that the wavelength 
of the hydromagnetic waves here considered is 
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extremely long; hence, the propagation must be 
quite different from that for very low frequency 
radio waves. A preliminary examination indicates 
that the above possibility is unlikely. 

Another evidence of propagation of hydromag- 
netic waves from the outer atmosphere is provid- 
ed by sudden commencements of magnetic storms. 
Examples of nearly elliptically polarized varia- 
tions have been found in some sudden commence- 
ments, and will be published later by C. R. Wil- 
son and the present author. 

The author is indebted to Professor S. Chapman 
for valuable discussions, to Mrs. Jean Chorbajian 
for her assistance in the analysis, and to the 
Coast and Geodetic Survey for making the data 
available. 
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X-RAY CONTINUA AND LINE SPECTRA FROM HIGHLY STRIPPED ATOMS 
IN A MAGNETICALLY COMPRESSED PLASMA* 
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A single-crystal x-ray spectrometer has been 
constructed and applied to the study of the x-ray 
spectrum emitted from the deuterium plasma in 
the Scylla fast magnetic compression experi- 
ment.’»? The Scylla magnetic field rises sinus- 
oidally to 55 kgauss in 1.25 usec and produces 
a deuterium plasma with an ion density of 5x10**/ 
cm® which emits neutrons and soft x rays for 
about 0.8 psec near the maximum of the second 
half-cycle of the magnetic field. 

The diffracting element of the spectrometer is 
a beryl crystal, ground for reflection from the 
1010 planes (2d =15.94 A), which allows examin- 
ation of the plasma emission spectrum between 
5Aand15A. The energy resolution of the spec- 


trometer is determined by the angular definition 
of the slit system, which is 10’ of arc for the 
measurements here reported. The diffracted 

x rays are recorded by a plastic scintillator — 
photomultiplier detector, and the time-resolved 
signal is displayed on an oscilloscope. The de- 
tector does not record individual x rays, but 
rather a pile-up signal whose amplitude indicates 
the time variation of the intensity of diffracted 
x rays produced during the 0.8-ysec emission 
interval at maximum plasma compression. 

The line spectra emitted by the deuterium dis- 
charge are shown in Fig. 1. The hydrogen-like 
spectrum of OVIII and the helium-like spectra 
of Na*, Mg¥!I, aiXIl and siXUI are well de- 
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FIG. 1. Spectrum of the deuterium discharge in Scylla. The arrows indicating the identified lines are at the 


predicted wavelength for each line. 
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veloped. The most energetic line observed is 
the SiXIII line corresponding to an excitation of 
1.84 kev. The impurities arise from the dis- 
charge tube wall which is 96% Al,O, and 4% 
oxides of Si, Mg, Na, etc. Linear Stark broad- 
ening in the Inglis-Teller theory of the merging 
of the levels at high quantum numbers’® is too 
small to account for the observed loss of OVIII 
Lyman lines at the eighth. It is more likely that 
the loss results from a combination of merging 
of the lines by Doppler broadening and instru- 
mental resolution and decreasing intensity of 
emission at large quantum numbers. An inter- 
pretation of the relative intensities of the im- 
purity lines from different elements is compli- 
cated by considerations of the relative refractor- 
iness of the constituent wall oxides, and hence 
the percentage constitution of the impurities in 
the discharge, and by the fact that the duration 
of the discharge is too short to produce steady 
state populations of the ionic species. 
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The line spectrum produced by the discharge 
when 10% Ne is added to the deuterium gas is 
shown in Fig. 2. These lines have not been re- 
corded previously, but their identification is 
certain since the wavelengths of the observed 
lines of the helium-like series can be predicted 
accurately by interpolation from other elements 
in the isoelectronic sequence.* The time-resolved 
x-ray detector, together with the pulsed Scylla 
source, provide a powerful technique for identi- 
fying the Ne species observed. The satellite 
lines labeled NeVIII peak earlier in time during 
the magnetic compression than the NeIX lines 
and are due to screened transitions of K-shell 
electrons in the more readily produced three- 
electron species.° The weak lines in the hydro- 
gen-like spectrum, NeX Ly and Ne* Lg, appear 
later than the NeIX lines. The intercombination 
line, 1s*S,-2p*P,, is nearly as strong as the 
corresponding line in the singlet series, as has 
been discussed by Edlén.* The line at 13.66 A is 
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FIG. 2. Spectrum of the Scylla discharge with 10% Ne added to the deuterium gas. The arrows 


are at the predicted wavelength for each of the identified lines. 


The probable origin of Ne lines 
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due to NelX, since it shows the same time his- 
tory as the other NeX lines. The most plausible 
assignment is 1s‘S,-2s*S,, where the normal 
selection rules are violated because of the strong 
ionic field. 

The spectral energy distribution of the con- 
tinuum during the second half-cycle is plotted 
in Fig. 3. Regions of the spectrum containing 
lines are omitted from this graph. Assuming 
that the electron energy spectrum is character- 
ized by a temperature T, and that the continuum 
consists of free-bound and free-free radiation, 
a plot of In(dE /dv) vs v should be a straight line 
of slope -h/kT,, except for small quantum- 
mechanical, Gaunt factor® corrections (less than 
10% in the frequency interval considered here) 
and possible discontinuities at the free-bound 
continuum limits of the various impurity ions. 
Changes in continuum intensity at the free-bound 
limits of oVill NaX, and MgX! are not observed, 
probably because the merging of the broadened 
spectral lines near each series limit obscures 
the change in intensity at its edge. The slope of 


the continuum in the case of a deuterium dis- 
charge with no added impurity corresponds to an 
electron temperature of 345+ 40 ev; the temper- 
ature for a deuterium discharge with 6% oxygen 
added is 295+30 ev. The greatest uncertainty 
in the temperature determination lies in the ex- 
perimental measurements of the beryl crystal 
reflectivity as a function of wavelength. A pre- 
vious analysis of the continuum spectrum by 
absorption methods indicated a somewhat lower 
electron temperature (240 ev).? This is because 
the newly discovered line radiation distorted the 
pure free-bound continuum below 14 A assumed 
for the previous temperature determination. 

A detailed description of these and related 
observations is being prepared for publication 
in The Physical Review. 





*This work was performed under the auspices of the 
U. S. Atomic Energy Commission. 

TConsultant; permanent address: Department of 
Physics, Cornell University, Ithaca, New York. 

1K, Boyer, W. E. Elmore, E. M. Little, W. E. 
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COSMIC -RAY FLARE OF NOVEMBER 20, 1960 


Richard T. Hansen 
High Altitude Observatory, University of Colorado, Boulder, Colorado 
(Received February 8, 1961; revised manuscript received February 27, 1961) 


In recent Letters to this journal, a solar flare, 
conjectured to have occurred 30° behind the sun’s 
west limb, was credited with production of a 
cosmic-ray neutron increase’ and a major 2800- 
Mc/sec solar noise burst.” We believe that at 
least the elevated portions of this flare, as well 
as associated prominences, were observed by 
us.* It is well known that flares are not re- 
stricted to the chromosphere, but typically occur 
elevated some thousands of kilometers above the 
chromosphere.*~* Based in the chromosphere 
30° beyond the sun’s limb, this flare would ap- 
pear if its radial height exceeded 93 000 km. 
While certainly rare, flares of comparably great 
heights have been reported in the literature.’»® 

The event first appeared at 1955 U.T. asa 
small mound on the limb of the sun and under- 
went a gradual increase in area and brightness 
until 2022 U.T. (see Fig. 1). The maximum 
brightness occurred between 2017 U.T. and 2024 
U.T., and it was only during this period that 
Lockheed and Sacramento Peak Observatories 
considered the event to be of flare brightness. 
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FIG. 1. The solar limb flare of November 20, 1960. 


West limb, 25° North. High Altitude Observatory, 
Climax, Colorado. 
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Abruptly at 2022 U.T. the feature began to as- 
cend, becoming detached from the limb by 2023 
U.T., and by 2032 U.T. appeared to have com- 
pletely disintegrated at a height of 300000 km 
above the sun. The transverse velocity was ap- 
proximately 900 km/sec. Because of the band 
pass of our filter, material having line-of-sight 
velocities in excess of about 100 km/sec would 
not have been observed. Another feature of 
flare brightness (described by the observer, 
Keith Watson, as a set of stationary, brilliant, 
coronal loops) appeared in the region at 2117 
U.T. and continued until the end of Climax’s ob- 
servations at 2257 U.T. 

The sequence of November 20, 1960 solar and 
associated events is given in the top of Table I. 
The following points seem most noteworthy: 

(a) The 2800-Mc/sec radio noise burst began 
with the start of the flare’s ascent (2022-2023 
U.T.) and continued long after the optical flare 
disappeared. 

(b) The ionospheric effects began between 2023- 
2028 U.T. and also continued long after the visi- 
ble flare. 

(c) The cosmic-ray increase began at 2055 
U.T. (£10 min), delayed about 33 minutes after 
the start of the disintegration. 

A strikingly similar sequence was also noted 
for another solar limb flare, believed to be the 
source of the May 4, 1960 cosmic-ray increase.** 
Some of the relevant observations are tabulated 
in the lower part of Table I. Kleczek and Krivsky 
reported a bright loop expanding with high trans- 
verse velocities (plus line-of-sight components 
of + 170 km/sec to -90 km/sec) at 1016 U.T. This 
was followed by the formation of intensely bright 
coronal loop prominences at 1040 U.T. Such 
prominences show the existence of regions of 
high density and strong, ordered magnetic fields 
in the sun’s corona, ingredients which are likely 
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Table I. Sequence of solar and associated events around the times of two cosmic-ray increases. 





Importance 


Start 


Maximum 


End 





November 20, 1960 


Flare (Climax) 3 
(Lockheed) * 1 
(Sacramento Peak)® 1 

Ejection phase (Climax) oe 

2800 Mc/sec” 300 % increase 

SCNA® (Boulder & Hawaii) 3 

sID* 1+ 

sIp4, d 3- 

SEA?) © 2+ 

Cosmic rays! eee 

Second flare (“bright loops,” Climax)® 3 

(Lockheed)* 1 
(Sacramento Peak)® 1- 


1955 
2017 
2019 
2022 
2023 
2028 
1945 
2023 
2023 
2055 +10 
2117 
2114 
2118 


2020 
2020 
2026.5 
2053 


2041 
2155 +10 
2145 
2135 & 2208 
2120 


2032 
2024 
2023 
2032 
after 2110 
2340 
2022 
2145 
2140 
after 2257 
2255 
2121 
after 2215 


(Sacramento Peak)® 2 
May 4, 1960 


Flare (NERA)* 3 

(Capetown)# 1 
Ejection (“expanding loops”’)& tee 
2980 Mc/sec (NERA)® 
SCNA® 
sip? 3 
SEA (Dunsink)# 1 
Cosmic rays eee 
Bright loops (Kleczek and Krivsk#)& 

(Capetown)* 


2126 2140& 2154 


1015 after 1105 
1000 a 1048 
1016 eee 


1700 % increase 1015 1118 


None reported 
1015 1050 
1016 1125 
1031.5 +1 eee 
1040 see after 1200 
1039 eee 1314 





2Reports submitted to IGY World Data Center A: Solar 
Activity, High Altitude Observatory and Central Radio 
Propagation Laboratory, Boulder, Colorado. 

DSee reference 2. 

Sudden cosmic noise absorption. 


necessary for the generation of cosmic rays by 
the sun. 

While the 1960 flares were observed only above 
the solar limb and may well have been a part of 
more extensive flares located behind the visible 
disk, two earlier cosmic-ray flares were ob- 
served near the limb (at about 70°-75° west lati- 
tude) but completely on the visible disk. Again 
these events were distinguished by high-velocity 
ejections at the start of the flare.’* For the 
flare of November 19, 1949, outward streaming 
prominence material was reported to a height 
of 560000 km above the limb and, in an early 
stage, a line-of-sight velocity component of 550 
km/sec was measured.** For the February 23, 
1956 event, Notuki et al. reported that “a surge 
of very high velocity was observed at the earlier 
stage of the flare.” ** Sketches of both of these 
ejections appear in the references cited. 

A more complete study of the characteristics 
of cosmic-ray-producing flares will be published 


dsudden ionospheric disturbance. 
©Sudden enhancement of atmospherics. 
See reference 1. 
&See reference 9. 
See references 10 and 11. 


in another journal. 

The superb observations of the sun on Novem- 
ber 20, 1960, are due to the Climax staff, 
D. Keith Watson, Robert James, Herbert Richey, 
and Robert Eddy. It is also a pleasure to ac- 
knowledge the strong contributions of Dr. James 
Warwick to this study. 





‘H, Carmichael, J. F. Steljes, D. C. Rose, and 
B. G. Wilson, Phys. Rev. Letters 6, 49 (1961). 

24. E. Covington and G. A. Harvey, Phys. Rev. 
Letters 6, 51 (1961). 

’The observations were made at the High Altitude 
Observatory— Climax with a Lyot-type H-alpha filter 
attached to a conventional coronagraph. A duplex 
optical system permits superposition on a single 35-mm 
frame of the solar disk observed at 0.75 A band pass 
with the over-the-limb features at 3 A band pass. The 
prominences are thereby photographed with a wider 
band pass than conventional solar flare patrol instru- 
ments. The instrument was designed by Dr. James 
Warwick and built at the HAO under support of the 
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U. S. Army’s Redstone Arsenal and the NASA. 
‘C. 8. Warwick, Astrophys. J. 121, 385 (1955). 
5J, W. Warwick, Astrophys. J. 121, 376 (1955). 
®R. G. Giovanelli and M. K. McCabe, Australian 
J. Phys. 11, 130 (1953). 
'H. W. Dodson and R. R. McMath, Astrophys. J. 
115, 78 (1952). 


~§p, H. Menzel, E. P. Smith, H. DeMastus, H. Ram- 


sey, G. Schnable, and R. Lawrence, Astron. J. 61, 
186 (1956). 

®J. Kleczek and L. KFivsky, Bull. Astron. Inst. 
Czechoslovakia 11, 165 (1960). 

1%. G. McCracken and R. A. R. Palmeira, J. 


Geophys. Research 65, 2673 (1960). 

"J, H. Trainor, M. A. Shea, and J. A. Lockwood, 
J. Geophys. Research 65, 3011 (1960). 

Similar high-velocity ejections have also been ob- 
served with flares for which there has been no report 
of cosmic-ray enhancements at the earth. For ex- 
ample: H. W. Dodson, E. R. Hedeman, and J. Cham- 
berlain, Astrophys. J. 117, 66 (1953); andJ. W. War- 
wick, Astrophys. J. 125, 811 (1957). 

43M’, A. Ellison and M. Conway, Observatory 70, 77 
(1950). 

“Mm. Notuki, Takeo Hatanaka, and W. Unno, Pubs. 
Astron. Soc. Japan 8, 52 (1956). 





DIFFUSION OF PLASMA ACROSS A MAGNETIC FIELD 


J. B. Taylor 
Atomic Weapon Research Establishment, Aldermaston, Berks, England 
(Received February 10, 1961) 


The diffusion of plasma across a magnetic field 
has been the subject of several calculations.'~* 
However, all these assume local thermal equilib- 
rium (i.e., the zero-order distribution is Max- 
wellian), although in some cases the ions and 


electrons may have different temperatures. These 


calculations all lead to a diffusion proportional 
to 1/B*. Bohm! suggested that there might be a 
form of diffusion in which fluctuating electric 
fields are present and produce a diffusion pro- 
portional to 1/B. In the course of work on the 


application of fluctuation theory to transport prob- 


lems in plasma,® a derivation of diffusion has 
been found which encompasses both classical and 
Bohm diffusion and indicates that, when suitably 


expressed, the Bohm formula gives the maximum 


value which the transverse diffusion can ever 
attain. 

We consider a situation in which a density gra- 
dient exists in one direction only, perpendicular 
to B, and where the magnetic pressure is domi- 
nant. The ion velocity will generally be much 
less than electron velocities so that the mean 
drag (dynamic friction) on an ion will be propor- 
tional to the ion velocity. Then we can write the 
Langevin equations for motion transverse to the 
magnetic field as 


u=-Bu+wu+X, 

b=-Bv-wu+Y, (1) 
where M5 is the coefficient of dynamic friction, 
w=eB/Mc, and MX/e, MY/e are the fluctuating 


electric field components. These fluctuations 
are assumed to be statistically independent of u, 
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v and to have a stationary distribution. 
The solution of these equations can be written 


u(t) ae Hy, cos wi+v, sinwt?) 


t 
. f 2®8-Dix (5) cosw(s-t) - ¥(s) sinw(s-t) Jas. 
o (2) 


Then, the mean square displacement per unit 
time is 


1 T ,T 
(Ax?) = a f f (u(t)u(t’))dt dt’, 
0 


which for T greater than a typical fluctuation 
time becomes 


T 
(A x*) = af (u,u(7)),dT. (3) 
0 


Evaluating the correlation coefficient from (2) 
and using the statistical independence of u,, vo, 
X, Y, one finds 


(A x?) = 2(ug?) B/(F? + w*). (4) 


Recalling that the contribution to the flux arising 
from (A x’) is 


F,=-4 2 [ax*ne)], 


we see that (A x”) is essentially the diffusion co- 
efficient. It can be verified that if we ascribe to 
6 the conventional value '** for a Maxwellian dis- 
tribution, then we recover the classical flux,’ 
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which for large magnetic fields is 


ria) See ©) 


However, the interesting feature of (4) is that, 
while 8 will depend on the electron distribution, 
no matter what value 6 may in fact take, (Ax) 
cannot exceed 


(Ax?) = (ug?)/w = §cW/eB, (6) 


where W is the mean ion energy. 

It will be recognized that, apart from a nu- 
merical factor, (6) leads to the value of diffusion 
suggested by Bohm, the value given by (6) being 


F = -i(c/eB)(ap/ ex), 


where we have taken p, the total pressure, to be 
twice the ion pressure. One concludes, therefore, 
that this diffusion is the maximum that can be 
attained with ions of a specified mean energy in 


a magnetic field of value B. 

It will be noted that we have not considered the 
effect of the first-order moment (Ax), which can- 
not easily be obtained from (1) since it depends 
on variations in 8. However, it can be demon- 
strated that, provided the mean ion and electron 
energies are equal, the first-order contribution 
vanishes for ions of Z =1. 
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CRYSTAL STRUCTURE OF THE § FORM OF He‘t 


R. L. Mills and A. F. Schuch 
Los Alamos Scientific Laboratory, University of California, Los Alamos, New Mexico 
(Received February 13, 1961) 


In 1953 Dugdale and Simon’ reported that upon 
investigating the thermal properties of solid He* 
they observed a first order transition whose 
equilibrium line cuts the melting curve at 14.9°K. 
The nature of the transition indicated a change 
in crystal structure. Since the calculated entropy 
and volume change associated with the transition 
were very small, they assumed that on crossing 
the equilibrium line from the low-temperature 
side the structure changed from the already 
known hexagonal closest packing, a-phase, to 
cubic closest packing, 8-phase. Whether helium 
solidifies in the cubic closest packed structure 
as the other inert gases do has been of theoreti- 
cal interest. Barron and Domb? from theoreti- 
cal considerations have shown that one might ex- 
pect a transition to the cubic form to occur in 
solid He* at an elevated temperature. We have 
therefore investigated by x-ray diffraction the 
structure of the new phase discovered by Dug- 
dale and Simon. 

The cryostat and camera arrangement pre- 
viously used*** for crystal structure studies at 
low temperatures was modified to permit oscilla- 


tion photographs to be taken. This was done by 
replacing the cylindrical Dewar by another which 
had the cross section of an annulus, thus per- 
mitting a shaft to be inserted down the geometric 
center of the Dewar. The upper end of the shaft 
projected through the top of the cryostat case 
where it connected to a cam and follower arrange- 
ment for oscillating it. As shown in Fig. 1, at 
the bottom of the shaft a beryllium cell was at- 
tached into which the helium was solidified. The 
cell had a 0.8-mm bore and a 0.4-mm wall. The 
open end of the cell was connected to the filling 
tube by a compression closure which consisted 

of a 59-degree cone tightened into a 60-degree 
seat with a brass gland nut. The differential in 
thermal contraction between brass and beryllium 
is such as to cause the joint to become tighter 
with lower temperature. The cell was cooled by 
the bath of liquid hydrogen through three braided 
copper straps. These and the filling and ther- 
mometer capillaries were attached with sufficient 
flexibility to allow the shaft to oscillate through 
an angle of 30 degrees. Not shown in Fig. 1 is 
the liquid-nitrogen-cooled copper radiation shield 
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FIG. 1. Arrangement of the x-ray diffraction cell. 
(a) Cam-oscillated shaft; (b) liquid H, bath; (c) vacuum 
space; (d) flexible copper braid; (e) H, vapor pressure 
bulb; (f) Be cell; (g) x-ray beam; (h) high-pressure 
conical seal. 


which surrounded the cell. The shield had win- 
dows for the entrance and exit of the x rays. 
The temperature at each end of the cell was 
measured with a hydrogen vapor-pressure ther- 
mometer. The two temperatures were averaged 
to give the temperature at the center of the cell. 
The pressure in the cell inlet capillary was meas- 
ured to + 10 psi with a 0-20000 psi Heise gauge 
previously calibrated with a free-piston gauge. 
The photographs were taken with a 4x5 inch flat 
plate camera placed 5 cm from the sample. The 
radiation used at different times was filtered Ky 
from an iron, copper, or molybdenum target. 
The usual procedure was to put helium at 1000 
psi into the cell and then to cool it to the neighbor - 
hood of 20°K by adding liquid hydrogen to the re- 
frigerant pot. The helium pressure was then 
raised to the amount desired for the experiment, 
18 440 psi, and the helium was next solidified by 
cooling the cell to 16°K. This latter cool-down 
required around 2§ hours. It was hoped that the 
slow cooling would promote the growth of strain- 
free crystals. In a single attempt to produce a 
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powder, the cell temperature was first lowered 
to 16°K, and then the helium pressure was sud- 
denly increased from just below the freezing 
pressure, 14000 psi, to 18440 psi. This proved 
to be unsuccessful because no more than a few 
large crystals were produced. During the x-ray 
exposures of from 23 to 7 hours, the pressure 
on the helium varied no more than + 20 psi and 
the cell temperature remained between 16.15 and 
15.45°K. 

In every exposure more than one crystal ap- 
peared in the beam. This had the advantage of 
giving a greater number of reflections in each 
photograph from which to derive the average 
value of the interplanar spacing, dppz, of each 
type of reflecting plane. , 

All of the reflections obtained could be indexed 
on the basis of a cubic lattice. The indices were 
always all even or all odd, indicating the face- 
centered cubic structure. 

The indices of the types of reflections obtained 
from one photograph, the interplanar spacings, 
and the corresponding length of the unit cell 
edge, A, are as follows: 












hkl dy py (A) A (A) 
111 2.447 4.288 
200 2.128 4.256 
220 1.487 4.205 
131 1.273 4.221 
222 1.230 4.251 


Average 4.232 


The average value of A from 7 photographs, 
weighted as to the number of types of reflections, 
is 4.240+0.016 A. The cell edge of 4.240 A cor- 
responds to a molar volume of solid helium of 
11.48 cc. This is 2.7% smaller than the value 
of 11.80 cc which can be computed from the data 
of Grilly and Mills® at a corresponding tempera- 
ture and pressure. 





Twork performed under the auspices of the U. S. 
Atomic Energy Commission. 
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(London) A227, 447 (1955). 

3a, F. Schuch, Atomic Energy Commission Report 
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‘a. F. Schuch, Proceedings of the Fifth International 
Conference on Low-Temperature Physics and Chemis- 
try, edited by J. R. Dillinger (University of Wisconsin 
Press, Madison, Wisconsin, 1958), p. 79. 

5g. R. Grilly and R. L. Mills, Ann. Phys. 8, 1 
(1959). 











wi 








VoLuME 6, NUMBER 6 





PHYSICAL REVIEW LETTERS 





Marcu 15, 1961 





NEW SOLID PHASE IN He‘t 


James H. Vignos and Henry A. Fairbank 
Sloane Physics Laboratory, Yale University, New Haven, Connecticut 
(Received February 20, 1961) 


Recent measurements of the longitudinal vel- 
ocity of first sound in solid He* have revealed a 
first order phase transition, and hence a new 
solid phase (designated the 7 phase), between 
1.45°K and 1.78°K in a narrow range of pressure 
adjoining the melting curve. The transition be- 
tween the previously established hexagonal close- 
packed a phase’ and the inferred y phase was 
detected by observing the discontinuity in the 
velocity of sound at the phase boundary. Subse- 
quent measurements of a change in molar volume 
across the phase boundary and a discontinuity in 
the slope of the melting curve at the lower liquid— 
solid a—solid y triple point have substantiated 
that this is indeed a first order phase transition. 

The standard ultrasonic pulse technique at a 
carrier frequency of 10 Mc/sec has been employed 
to measure the longitudinal velocity of sound in 
solid as well as liquid He*. First order phase 
transitions can also be explored since there is a 
discontinuous change in velocity on crossing the 
transition. The melting curve of He* has been 
investigated in this manner and is shown in Fig. 1 
(solid circles). Investigation of the velocity of 
sound in the solid region near the melting curve 
indicated the existence of another first order 
phase transition. The locus of the discontinuity 
in velocity observed in this region is shown in 
Fig. 1 (open circles). This discontinuity was 
found to be reversible in both temperature and 
pressure. The locations of the liquid—solid a— 
solid y triple points are given by T,, = 1.778 
+0.003°K, P,,=30.28+ 0.05 atm and T;=1.449 
+0.003°K, P;=26.18+0.05 atm. The velocity of 
sound in the solid a region, both along the melt- 
ing curve from 1°K to the lower triple point and 
along the a-y phase boundary, was found to be 
essentially constant and equal to 478 m/sec. 
Throughout the solid y region the velocity has 
values ranging from 520 m/sec to 545 m/sec, 
these two extremes being the predominant values. 
We attribute this variation to different orienta- 
tions of a single crystal or a few large crystals. 
This tendency for helium to form large crystals 
is consistent with previous observations.’~> For 
reference the velocity of sound in the liquid at 
the melting curve is 365 m/sec from 1°K to the 
A point and then climbs to 380 m/sec at 1.8°K. 
The lambda line shown in Fig. 1 (solid triangles) 





Solid a 


28 


Pressure (Atm) 


26 











24 
1.2 1.4 1.6 1.8 


Temperature (°K) 


FIG. 1. The phase diagram of He‘ showing the 
melting curve (closed circles), the solid a-solid y 
transition (open circles), and the A transition (tri- 
angles). 


was determined by noting points of maximum 
sound absorption in the liquid. The lambda line-— 
melting curve triple point is given by T, =1.765 

+ 0.003°K, and P, =29.90+0.05 atm. From 1°K 

to 2°K our melting curve and lambda line lie 
about 0.2 atm above those determined by Swen- 
son.* 

The gas-handling system employed allows a 
measure of the flow of gas into or out of the sam- 
ple chamber; consequently we were able to meas- 
ure, though very roughly, the change in molar 
volume in passing from one phase to the other. 
The volume of the y phase was larger than that 
of the a phase. At 1.50°K the change in volume 
at the transition was 0.07 cc/mole, which cor- 
responds to a 0.3% change in the molar volume. 
The slope of the a-y transition curve (Fig. 1) at 
1.50°K is 9.3 atm/°K and hence by the Clapeyron 
equation the difference in entropy of the phases 
(S, - Sq) is 0.016 cal/deg mole. 
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An enlarged graph of the phase diagram in the is immeasurably small. 
vicinity of the lower triple point is shown in The structure of the new y phase in solid He* 
Fig. 2. It is seen that at the lower triple point has not yet been established; however, there is 
there is an abrupt change in slope of the melting reason to believe that it is body-centered cubic. 
curve. From the Clapeyron equation applied to A phase transition has been observed in solid 
each transition at the triple point it is seen that He*, above about 100 atm,° which is similar to 


such a discontinuity is required, its magnitude the a - y transition in solid He* in its pressure 
and temperature dependence and in the change in 


being given by 
molar volume. In He* the higher pressure phase 
aP “ap (8) is hcp and the lower pressure phase (qa) is 
=) ( ) bec.® In the case of He* the higher pressure phase 
a,l yl (a) is also hcp, and hence one is led to believe 


by analogy that the lower pressure phase (y) is 
most likely bcc. Why a stable bcc lattice is 
{( ar) { ar) Juv /av_»), (i) found in He® (and possibly He*) and in no other 
aT , aT yt vb GY solid dielectric element is not yet clear.”»* It 
is very unlikely that nuclear spin effects are 


where the subscripts 1, a, andy refer to the responsible for the transition in He* if the new 


liquid, solid a, and solid y phases, respectively, |» phase in He* does prove to be bcc. 
and AV,. x,y* V,-V,. Using our values of slope The first order phase transition to the 6 phase 


for the two liquid- solid transition curves, our in solid He* above 15°K and 1100 atm, inferred 
value of AVq,,, at 1.50°K, and Swenson’ s* value from the specific heat measurements of Dugdale 
of AV,,,] at the triple point, we have calculated and Simon,” is probably a transition to a cubic 
from Ea. (1) the value of (dP/dT), y (shown by close-packed structure and is unrelated to the 
the dashed curve in Fig. 2) and find it to be con- low-pressure transitions in solid He* and He* 
sistent with the experimental value. At the upper discussed above. 

triple point Eq. (1) would predict a change in In summary, a new y phase has been shown to 
slope of the melting curve of only 2%, assuming exist in solid He*. It has a larger volume and 
AVa,y to be about that observed at 1.5°K. This entropy than the adjacent hexagonal close-packed 
is in agreement with the observed data of Fig. 1 a phase, as required by thermodynamics, and 
where the change in slope of the melting curve a larger longitudinal sound velocity. By analogy 
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with solid He* it seems likely that the new y phase 
is body-centered cubic. 
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COMPRESSIBILITY, ZERO- POINT ENERGY, AND SPECIFIC HEAT IN SUPERCONDUCTORS* 


J. G. Daunt and J. L. Olsen! 
Ohio State University, Columbus, Ohio 
(Received February 9, 1961) 


Keesom and Bryant’ showed that the specific 
heat of indium in the superconducting state at the 
lowest temperatures is considerably smaller than 
the lattice specific heat deduced from measure- 
ments in the-normal state. The work of Boorse, 
Hirschfeld, and Leupold and of Zavaritskii® in- 
dicates a similar effect in niobium and tin. Such 
a phenomenon would be easily understood if there 
were correspondingly large differences in the 
elastic constants and in the Debye characteristic 
temperatures, @p, of the two states. Work by 
Chandrasekhar and Rayne* has now confirmed the 
generally accepted view in ruling out the possi- 
bility of any differences in 0p greater than 1 part 
in 10*. We wish to show here that even the very 
much smaller changes in 6p which are known to 
exist can influence the zero-point energy enough 
to allow an explanation of the observations. 

The total energy, E, of a Debye solid at low 
temperatures may be written 


E=1.125 R60, +58.5RT*/0 *s (1) 


where R is the gas constant. The first of these 
terms is usually known as the zero-point energy 
of the lattice. If @p is constant, we obtain the 
well-known expression for the T° lattice specific 
heat, C Pa at low temperatures: 


C =8E/8T =234RT°?/0_°. (2) 
g D 


An attempt to explain the observed differences 
in specific heat using this formula requires com- 






plicated and unlikely changes in the phonon spec- 
trum. 

On the other hand, it is well known that the 
difference in compressibility between the two 
states is temperature dependent and given by® 


aV_ av oH * H #H 

(ee) ££ @) 

V\ 8p ap] 4n\ a | 42 ap?’ 
where V,, and V, are the volumes in the normal 
and superconducting states, H,, is the critical 
magnetic field, and p is the pressure. The com- 
pressibility in one of the states at least must 
then be temperature dependent, and @p must 
clearly be so, too. 

Such a temperature dependence of @p- affects 
the zero-point energy and contributes to the spe- 
cific heat, which now becomes 
Cc E* aE /8T 

= 1.125R(80 ,/8T) + 234RT°/0 ° 


x[1- (37/40 ,)(80 ,/2t)]. (4) 


It is obvious that even a very small value of 
80p/8T at the lowest temperatures may produce 
effects comparable with the value of C g given by 
Eq. (2). 

The expressions (3) and (4) and the known de- 
pendence of 6p on compressibility can be used 
to calculate the difference between the lattice 
specific heats, C,,, in the normal and C gs in the 
superconducting states, from the pressure de- 
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pendence of the critical field. Straightforward but tedious algebra leads to 


Cc -C = 
gs “gn ap 


where K is the bulk modulus, p is the pressure, 
H, is the critical field at T=0, y*T is the elec- 
tronic specific heat per unit volume, t=T7/T,, 
and T, is the transition temperature. 


F(t) = (2a- a?) f'(t) F(t) - °F - tf”) (ED 

+ (1-aP [3A F" (F(t) +P F(t) F(2)), (6a) 
G(t) =f’ (F(t) + F(t + "OF, (6b) 
H(t) =f" F(0) - (EP - "OF, (6c) 
J(t)=(1+ a) fF’) - (1- adler + FOF"), 


(6d) 
1/1 ay"\ /{ 1 ate 
e=2\y* a H, @/)’ 
and f(t) is defined by H, =H, f(t). Primes denote 
differentiation with respect to ?. 
These expressions are greatly simplified if it 


is assumed that f(t)=1-. Then 


where 


F(t) = 40°(3-2a)(1-2a), (7a) 
G(t) = 88 - 4t, (7b) 
H(t) =-4#, (7c) 
J(t) = -4t9(1-2a) - 4at, (7d) 


Table I. Data for calculation of the term B of Eq. (8) together with observed values of B for various supercon- 


-9RO 7-822)" ik | (=) F(t)+ 





H,2 &y* eH 2. oH 
7 cin, Eemo) +360, (5) ‘ 
and Eq. (5) reduces to the form: 
C -C_=AT+BT* (8) 
&s gn 


The coefficients of F(t) and J(¢) can be obtained 
since 4H,/%p and K are known. The value of a 
is somewhat less certain, but information about 
its value is available.*»’ The coefficient of H(t) ‘ 
which involves 8*H,/ 2p can be estimated for tin 
from work on the differences in the velocity of 
sound in the normal and superconducting states.°,’ 
No data are available at present allowing more 
than a guess at 8*)*/8p* which appears in the 
coefficient of G(t). Work on the difference in the 
modulus of rigidity’® seems to suggest that in tin 
this is about one-tenth of #H,/ 9p. 

Table I lists available information for various 
superconductors. Bac gives the coefficient of 
the term in 7° in Eq. (8), under the further sim- 
plifying assumption that 0H,/ep? = 9*)*/ap? = 0. 
Bobs is that obtained from the observed specific 
heats. 

The order of magnitude and sign of B,gj- en- 
courages the belief that the specific heat anoma- 
lies are, in fact, created by a temperature de- 
pendence of the zero-point energy of the lattice. 














ductors. 
Nb Pb Sn In Al ne 

T, CK) ~ 7.2 3.73 3.37 1.2 
@ (°K) 250 90 160 109 419 

a 1.5 40.5 0.10 0. 05 0 +0.03 0.25 
$K(8H,/8p)*x 10° 17 16 5 1.6 3.3 
(4K Hy*/y*)(8?-y */ap*) x 10° 0+3 ne one 
3K H,(8?H,/8p*) x10° +18 tee tee 
4Hy(8H,/8p) x10° -5 -4 -1.0 -0.6 -0.2 
Bobs (ujoule/mole deg) -22 -20( ?) -3002 ove 

B calc (ujoule/mole deg) -1.3 -4 -50 -12 -2500 





®The observations on In may perhaps better be described using A=-30 and B =-10. 
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For a further advance, more data on a, 9H,/@, 
#H,/ ep, and &*/op? for various directions of 
stress are required. 

Our calculation gives rise to a linear term, A, 


in C gs ~ C,,,- The sign and magnitude of this 
term depend upon the coefficients of G(t) and 
J(t) in (5). The experimental data indicate that 
the coefficient of J(t) in general makes a small 
positive contribution to A. The coefficient of 
G(t), which involves #*/ap?, is not known, but 
it might be negative in some cases. If so, and 
assuming that the parabolic form of f(t) remains 
valid, formally one deduces an apparent negative 
specific heat in the superconducting state at the 
lowest temperatures. It is not clear what modifi- 
cations or mechanisms should be invoked in order 
to avoid this conclusion, since existing data do 
not extend to low enough temperatures. Bryant 
and Keesom’s results indicate A <0, but in their 
case a positive nuclear quadrupole term keeps 
the total specific heat from becoming negative. 

We should point out that the role of the zero- 
point energy in causing the observed specific- 
heat anomalies was realized independently by 
Schrieffer who has approached the problem from 
a microscopical point of view. 

We are most grateful to Dr. J. R. Schrieffer 
and also to Dr. B. S. Chandrasekhar and Dr. J. A. 


Rayne for telling us of their results prior to pub- 
lication. We also wish to thank Dr. H. Rohrer 
for providing an estimate of 6H,/®p in niobium 
from his unpublished measurements. 
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OSCILLATIONS IN THE LONGITUDINAL TUNNEL CURRENT OF TUNNEL DIODES 


R. R. Haering and P. B. Miller 
International Business Machines Research Center, Yorktown Heights, New York 
(Received January 18, 1961) 


It has been observed by Chynoweth, Logan, and 
Wolff: that the tunnel current in InSb tunnel di- 
odes displays oscillations at low temperatures 
and large longitudinal magnetic fields (H |IF). The 
field dependence of the periodicity of these os- 
cillations suggests a relation to the de Haas- 
Schubnikow oscillations in galvanomagnetic pro- 
perties. These oscillations are due to density of 
state fluctuations in the highest occupied Landau 
level. Such oscillations are not to be expected 
for transverse magnetic fields,'»? and none were 
observed.! 

If the observed oscillations are of the de Haas- 
Schubnikow type, their magnitude suggests that 
a significant portion of the tunneling current is 
due to electrons with the highest possible mag- 
netic quantum number. This conclusion is sur- 





prising, since the tunneling process strongly 
favors the low quantum numbers. The depend- 
ence of the tunneling rate on quantum number is 
given by” 


exp[- 3(nhw/éa|, 


where 
d= &(2me)**/hmF, 


F =[u,+u,+¢€ - qV|/a=average junction field, 


a=junction width, n=quantum number, w =eH/ 
mc, €=energy gap, m =electron-hole reduced 
mass, V =applied voltage. For the particular 
InSb diode used in reference 1 we have \ ~ 30. 
Hence, the ratio of the contributions of the 
largest and smallest quantum numbers is 

~ exp(-45u,/e), where yp, is the electron Fermi 
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level. It is therefore clear that the observed 
oscillations are far too large in amplitude and 
in number to be explained in this manner. Fur- 
thermore, the above process only leads to cur- 
rent discontinuities and not oscillations. 

We wish to propose an alternative explanation 
of the observations based on the dependence of 
the electron Fermi level on magnetic field. In 
a magnetic field, the electron Fermi level os- 
cillates about its zero-field value. (Fluctuations 
of the hole Fermi level are much smaller be- 
cause of the heavy-hole states.) These small 
Fermi level fluctuations result in a small charg- 
ing or discharging of the junction capacitance 
and hence in a small oscillatory change of the 
average junction field F. Since the tunneling rate 
is a strong function of F, small fluctuations in F 
can result in large current oscillations. It is to 
be noted that whereas the origin of the Fermi 
level fluctuations is related to density of state 
fluctuations arising from the highest occupied 
quantum number, the resulting change in the 
junction field affects the tunneling rate of all 
quantum states. We shall show that the above ef- 
fect accounts for the observed period of the os- 
cillations, as well as for the magnitude and field 
dependence of the oscillation amplitudes. 

For an isotropic conduction band with constant 
effective mass m,, the Fermi level yu,(H) is de- 
termined by the condition: 


Ny 12 
${HO™ - > [4 -n-4-3 BH 
SLhw,) ~ pvol hw, "7 24th, 


hw, 281 hw, 
where @ is the Bohr magneton, g, is the conduc- 
tion band g factor, and N, is the maximum integer 
<([,(4)/hw,|-3*+42,RH/hw,. In the extreme quan- 
tum limit (u,/hw,< 1), u,@) is readily obtained 
directly from Eq. (1), whereas for p,/hw21 we 


N. v2 
+> [Re -n-a+8 st : (1) 
n=0 


obtain*® 
u, (A) 1] hw, |* u,(0) , 1, m 
ee a 


ved mt (2) 
0 


The function E(x) in Eq. (2) has approximately 
unit amplitude and is periodic with period one. 
For 0<x <1, we have 


E(x) =2(¢ + 3)%* - 3 +1) - 3x™. (3) 
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The oscillatory part of Eq. (2) is plotted vs 

ut, (0)/hw, in Fig. 1 for g, = 56 and m,/m,=0.013. 
For a junction for which a ~|p,+ pu, -qV+e)”%, 

a change Ay, in the electron Fermi level results 

in a change AF in the average junction field, 

given by 


#.(% *) Au, (4) 
F a /(u,+u,-QV+e) 


The number a is normally ~2, the value ob- 
tained for a step junction. For small biases qV, 
the tunneling exponential \ is therefore 


n=ro]t-(*3)( Fi FI (5) 
@ J\u,+u,+e) pv, 


Since the tunneling current density for small 
biases is 


j(A) = const ve * fr 2) csen (Me 2) 


2 8 
x cosh (E i): (6) 
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FIG. 1. Plot of oscillatory part of Eq. (2) with 
a1 = 56, m,/my = 0.013. 
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where E, = 3(g, - g,)BH, the oscillatory part of 
the tunneling current is given by 


JH) -5, 8) bs hw ap ” hw ao 
anepegice ONO <a> <a Se ae 
jay) 2 ve 2 

rn ar 

0 0 

So eT hee 





(Fe) (arin): © 
a Mit+bate)\ bw, 


Here Au,/p, is to be obtained from (2), and 
jay) is obtained from (6) with \ replaced by 
No 

Under normal conditions for observation of os- 
cillations, the coth and tanh terms are negligible 
and we find for InSb (with a = 2) 

. a . 92 
j(4) Jay) do “4 hw} 


Jay) 12 HytHotre]\ Hy 


«{ (fs o.82) += (fs o09] (8) 
W, hw, 


For p,/hw,~1.5, 5~30, u,/(u, + U2 +€)~ 3, the 
oscillations are therefore ~30% of the average 
current. For smaller fields the relative im- 
portance of the oscillations is seen to decrease 
like H**, Temperature and collision broadening 
of the sharp peaks in Fig. 1 are ignored in this 
estimate, which would tend to reduce the ampli- 


tudes. This damping effect also accounts for the 
appearance of only one current peak per cycle. 
The observations of Chynoweth et al. are in 
good agreement with the above results. 

It should be noted that the above effect is not 
to be expected for a transverse magnetic field. 
In this case the discontinuities in the density of 
electron states for states whose orbit center 
lies less than one orbit radius from the junction 
no longer occur at particular values of the mag- 
netic field, but are distributed more or less 
continuously and hence the average junction 
field remains unaltered. Density of states fluc- 
tuations due to states whose orbit center is far 
from the junction do result in the appearance of 
a weak space-charge region on the m side. How- 
ever, this space-charge region is not in the 
vicinity of the junction and hence does not affect 
the junction field. 

In previous theories it was customary to ignore 
the shift of the Fermi level with magnetic field.?>* 
The present analysis indicates that although this 
shift is small, it plays an important role in the 
tunneling process. 

The authors would like to thank Dr. E. N. 
Adams for several helpful discussions. 
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STIMULATED SPIN-ECHO MEASUREMENT OF CROSS-RELAXATION 
IN NEUTRON-IRRADIATED CALCITE* 


L. Kent Wanlass and J. Wakabayashi 
Electronic Research Laboratory, University of California, Berkeley, California 
(Received December 19, 1960; revised manuscript received February 9, 1961) 


In recent years, possibly due to the advent of 
microwave spectroscopy and maser applications, 
a great deal of research has been conducted on 
the relaxation times of spin systems. Presently 
there is a large amount of work being conducted 
on cross-relaxation decay’ and its subsequent 
effects. The majority of these experimental 
measurements have been conducted by saturation 
methods or pulse modifications of them, which 
lead to exponential decay with two or more time 





constants. The study of electron spin echo as an 
information storage medium for high-speed car- 
rier-type microwave computers has led the au- 
thors to some interesting measurements of re- 
laxation times in neutron-irradiated calcite. Of 
more significant interest, however, is the use 

of the stimulated spin-echo technique in cross- 
relaxation measurements. This method, in cases 
where the T, relaxation time is somewhat larger 
in magnitude than the cross-relaxation time, 
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FIG. 1. Oscilloscope photograph of a stimulated 
spin echo, Horizontal scale 5 usec/cm, vertical 0.01 
v/com, temperature=1.5°K. The pulses from the 
klystron show as breaks in the base line. 


gives a single exponential decay curve due only 
to the T,, time constant. 

The experimental arrangement used is very 
similar to that of Gordon and Bowers,” consisting 
of an X-band heterodyne spectrometer and a pro- 
grammed pulse generator. The pulsing system 
was made up of six Tektronix pulse generators 
connected into a single General Radio pulse am- 
plifier. The latter supplied negative pulses to 
the klystron repeller. The pulses were variable 
in width, separation, and also in repetition rate 
of the entire pulse group. The narrowest pulses 
obtainable from the programmed pulse generator 
were 1 usec long. This hindered the measure- 
ments to some extent, but in our case versatility 
was felt to be more important. The local oscil- 
lator was used in the early studies of the spin- 
spin and spin-lattice relaxation times, but was 


later eliminated for some of the tests. 

The stimulated echos obtained here at a fre- 
quency of 8807 Mc/sec are a paramagnetic analog 
of those obtained at nuclear frequencies by Hahn.* 
Figure 1 is a typical oscilloscope photograph of 
the stimulated spin-echo signals. The microwave 
pulses from the klystron drive the oscilloscope 
off scale, and are therefore merely indicated by 
breaks in the base line of the trace. All of these 
pulses were 1 psec long and of an intensity to 
give the spin system a rotation through an angle 
of approximately 7/2 radians. The first echo is 
the familiar 8-ball echo (so called because of the 
pattern of the spin vectors on a sphere’). This 
echo relaxes in a time characterized by spin-spin 
relaxation or T,. In our radiation- damaged cal- 
cite sample, this was approximately 18 psec at 
1.5°K.* However, the stimulated echo following 


the recollection pulse is of much greater interest. 


As pointed out by Fernbach and Proctor,° this 
signal is due to storage parallel to the direction 
of the Z axis and, consequently, its decay nor- 
mally depends only on spin-lattice relaxation. 

In our case, the spin-lattice relaxation time of 
the calcite sample was measured both by adia- 
batic fast passage and by saturation methods, 
giving a time constant of T,=5+2 seconds at 
1.5°K.* However, the stimulated echo decayed 
with a time constant of approximately T= 600 
psec or nearly 8000 times faster than predicted 
by T,. This rapid relaxation is undoubtedly due 
to cross-relaxation,' since the calcite spectrum 
consists of three spectral lines with a separation 
of about three gauss, when equally spaced. Fig- 
ure 2 shows a decay curve obtained by means of 
the stimulated spin-echo technique. The solid 
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line is a theoretical plot of an exponential decay 
exp(-t/6 x10-*), while the dashed curve repre- 
sents the level to which the energy stored in the 
line under observation decays by cross-relaxa- 
tion. The energy of the three lines then decays 
simultaneously from this level with a time con- 
stant equal to T,. No attempt was made to fit a 
curve to the experimental points of the echo sig- 
nal amplitudes, but the smaller values obtained 
at the start of the decay can be justified by the 
method used to measure the echos. 

It is of importance to note that the spin-echo 
decay continues at the cross-relaxation rate even 
after the three lines of calcite have come to en- 
ergy equilibrium. This is to be expected on the 
basis of the model proposed by Bloembergen 
et al. Since the balance of energy between the 
three lines has the characteristics of a dynamic 
equilibrium and the very nature of cross-relaxa- 
tion allows mutual spin flips between electrons 
with somewhat different resonant frequencies to 
transfer energy between spectral lines, the fre- 
quency storage of the spins is destroyed by the 
cross-relaxation process. The stimulated echo 
depends on the Z-axis storage of this frequency 
dependence and the echo signal consequently de- 


cays exponentially to zero with the cross-relaxa- 
tion time constant. 

This same process can occur within a single 
inhomogeneously broadened line and might make 
a very interesting study, since the stimulcied 
spin-echo technique effectively separates the de- 
cay due to cross-relaxation from that due to spin- 
lattice relaxation when the time constants are not 
comparable. . 

The authors would especially like to thank 
Professor J. R. Singer for suggesting this pro- 
ject and for his continuing support and helpful 
guidance. Thanks are also due to Professor 
E. L. Hahn as well as Dr. J. Kemp, S. Hartmann, 
and M. Weger for helpful discussions. 
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DISLOCATION RELAXATION SPECTRA o 
OF COLD-WORKED BODY-CENTERED CUBIC TRANSITION METALS 


R. H. Chambers and J. Schultz 
John Jay Hopkins Laboratory for Pure and Applied Science, 
General Atomic Division of General Dynamics Corportion, San Diego, California 
(Received December 5, 1960; revised manuscript received January 30, 1961) 


During the course of a low-temperature acous- 
tical study of the defect structure introduced in- 
to bec transition metals by cold work, peaks 
have been found in the internal friction vs tem- 
perature curves of Nb, Ta, Mo, andW. The 
peaks appear only with prior plastic deformation 
and have a detailed shape which is controlled by 
the annealing history of the specimen. 

Figure 1 shows the internal friction vs tem- 
perature curves for fine-grained vacuum- melted 
Nb, Ta, Mo, and W specimens given a 3% ten- 
sile prestrain at 300°K (~600°K for W). An in- 
ternal friction spectrum of a typical fcc metal 
(Cu) prestrained 3% is included for comparison. 
Most of the acoustical measurements were made 
by the longitudinal resonant bar method at 15 kc/ 
sec and oscillating strain amplitudes of ~ 10-’; 
several measurements on Mo were made at 6 
cps, at a strain amplitude of 2<10~" in a torsion 
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FIG. 1. Damping spectra of polycrystalline Cu, Nb, 
Ta, Mo, and W measured at ~15 kc/sec. All speci- 
mens were given ~3% tensile prestrain at 300°K 
(~600°K for W). 
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pendulum. It can be seen that at 15 kc/sec the 
Nb and Ta spectra have a large peak at 170°K 
and 180°K, respectively, and a much smaller 
one near 310°K. The Mo spectrum at 15 kc/sec 
is divided into two broad peaks, a low-tempera- 
ture peak (a) at 115°K and a higher temperature 
peak (8) near 400°K; the 8 peak is located at 
240°K for measurements made at 6 cps. The 
15-kc/sec W spectrum consists of an a peak at 
170°K and a 8 peak which appears to exist some- 
where above 450°K. 

Figure 2 shows evidence of structure in both 
the a and § peaks of fine-grained Mo prestrained 
in tension one percent. The observed structure 
resulted from a series of one-hour vacuum an- 
neals at successively increasing temperatures. 
The @ peak begins to break up into at least three 
subpeaks located at about 260°K, 320°K, and 
400°K for a 300°C, one-hour anneal, then vir- 
tually disappears after a 400°C, one-hour anneal. 
The a@ peak shifts to lower temperatures and 
shows evidence of structure with a slight shoul- 
der on its high-temperature side. The lower 
peaks in both Nb and Ta shift to lower tempera- 
tures on low-temperature annealing. 

There is additional evidence for believing that 
the observed internal friction spectrum in Mo is 
the sum of a number of subpeaks. The ratio of 
the quantity (Y - Y’)/Y,, where Y’ is the relaxed 
component, Y the elastic component of the 
Young’s modulus, and Y, the elastic component 
of Young’s modulus measured at 5°K after a 
800°C anneal, to Q,,"', where Q,,~' is the height 
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FIG. 2. Effect of successive one-hour vacuum an- 
neals on the damping spectra of polycrystalline Mo 
measured at ~ 15 kc/sec. The specimen had received 
~ one percent prestrain at 300°K. 
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of the internal friction peak, is about 8; where- 
as a single relaxation process would yield a ra- 
tio of 2.4 Also, careful examination of the 
Young’s modulus vs temperature curve shows 
the existence of a number of slight but distinct 
discontinuities in the slope of the curve at tem- 
peratures corresponding to the location of the 
suspected subpeaks. These discontinuities are 
similar to those found in cold-worked Cu single 
crystals by Thompson and Holmes.” 

All four of these metals also show the strong 
amplitude dependence in the internal friction and 
Young’s modulus characteristic of pinned dislo- 
cation motion.* It is observed that the annealing 
of the peaks is always accompanied by a strong 
(10- to 50-fold) increase in the breakaway ampli- 
tude. This effect was first noted in earlier work 
on Nb by the authors‘ and attributed to the pin- 
ning of dislocations by the diffusion of point de- 
fects (gaseous interstitials) to dislocations. The 
present data indicate that this mechanism can 
be extended to include Ta, Mo, and W. 

The peaks observed here are similar in sev- 
eral respects to those produced by cold-working 
fcc metals (the Bordoni peaks)*: Both possess 
strong dependence on cold work, have a wide 
distribution of relaxation times, and both occur, 
for a comparable measuring frequency in a sim- 
ilar reduced temperature range, i.e., tempera- 
ture normalized to the melting point of the met- 
al. Furthermore, the attempt frequency, /,,° 
in the relation f=/, exp(-H/RT) for the 8 peak 
in Mo is 10*%**°4 sec-!, which is comparable 
to that found for the Bordoni peaks in fcc metals. 
On the other hand, the peaks in bcc metals all 
anneal at a considerably lower reduced tempera- 
ture than do the peaks in fcc metals; and the 
bec a-peak positions show a marked shift to- 
ward lower temperatures with annealing in con- 
trast to a similar but smaller shift on annealing 
found in the Bordoni peaks of fcc Cu. 

It has also been observed that variations in the 
concentration of gaseous impurities affect princi- 
pally the annealing characteristics of a given 
peak, while the peak temperature and height re- 
main essentially unchanged. A more complete 
account of these impurity effects and their impli- 
cations regarding the nature of the defect mech- 
anisms operating here will be published else- 
where. 

On the basis of the above observations, the 
authors consider that the defects responsible for 
the relaxation spectra (especially the 8 peaks) in 
cold-worked bcc metals are similar to those that 











VoLUME 6, NUMBER 6 





PHYSICAL REVIEW LETTERS 


Marcu 15, 1961 





produce the Bordoni peaks in fcc metals. Fur- 
thermore, we feel that the different response of 
the bcc metals to annealing can be attributed to 
the action of efficient dislocation pinning points 
(gaseous interstitials and/or carbon) known to 
exist in bcc metals. 

We wish to thank P. H. Miller, Jr., for en- 
couragement and support and E. Moxley for her 
technical assistance. 
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ABSENCE OF AN ISOTOPE EFFECT IN SUPERCONDUCTING RUTHENIUM 


T. H. Geballe, B. T. Matthias, G. W. Hull, Jr., 


and E. Corenzwit 


Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received February 17, 1961) 


Within the accuracy of our measurements, 
(0.001°), naturally occurring ruthenium and its 
isotopes of mass 99 and 104 have identical super- 
conducting transition temperatures, T,,. Our 
results are presented in Fig. 1. This lack of 
dependence of T,,. on atomic mass is quite dif- 
ferent from the results obtained for the isotopes 
of other superconducting elements.’ All elements 
which have been previously measured, Sn, Hg, 
Tl, and Pb, show an isotope shift of T,, which 
is essentially proportional to M~”*. However, 
they are all nontransition elements, having s 
and p electrons only. Ruthenium is the first of 
the transition elements to be checked with re- 
spect to an isotope shift. The reason for this 
lack of data may have been-—in part, at least— 
experimental difficulties. The superconducting 
properties of ruthenium and osmium have been 
discovered and reported by Goodman’ and Hulm 
and Goodman.* With the exception of these two 
metals, transition element superconductors have 
critical temperatures which change strongly 
with small amounts of dissolved nitrogen or 
oxygen. These variations exceed the shift that 
may be expected from an isotope effect. There 
is no good way at present to get rid of these gases 
in a quantitative way with the small amounts of 
material available. 

Ruthenium isotopes were obtained in powder 
form from Oak Ridge National Laboratory through 
the courtesy of P. V. Arow. The isotope analysis 
furnished is given in Table I. After our initial 
results, a second analysis was made at the Oak 
Ridge Laboratory which was in good agreement 
with the first. Spectrographic analysis at Oak 
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FIG. 1. Change in frequency of oscillator circuit as 
a function of temperature for arc-melted pellets of 
ruthenium. Samples Ru® and Ru! were made from 
enriched powders of Ru® and Ru’; sample J&M was 
made from Johnson-Matthey powder and sample E 
from Englehard Industries’ natural powder. 


Ridge indicated an absence of all transition metal 
impurities which might falsify the results, with 
the possible exception of iron. However, a fur- 
ther analysis indicates less than 10 parts per 
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Table I. Isotopic analysis of Oak Ridge ruthenium. 

Atomic Atomic 

Ruthenium Isotope percent Ruthenium Isotope percent 
99 96 <0.02 104 96 <0.03 
Series IL 98 0.38 Series IL 98 <0.03 
Sample 1017a 99 89,82 Sample 1021la 99 0.09 
100 6.71 100 0.11 
101 1.49 101 0.28 
102 1.23 102 3.18 
104 0.34 104 96 .35 








million of iron.* Such a concentration can only 
reduce Tg¢ by about two-tenths of a millidegree. 

The powders were measured as received and 
showed broad (~ 0.08°) transitions commencing 
near 0.48°. Somewhat sharper transitions near 
0.48° were obtained with natural powders. Elec- 
tron diffraction pictures of natural ruthenium 
showed a size distribution of the powder particles 
ranging down to 500 A.® This result indicates 
that the broad transitions were due to the small 
particle size, since 500 A is in the order of the 
penetration depth. Therefore, the powder was 
melted into a pellet in an arc furnace. The pre- 
liminary results with natural ruthenium imme- 
diately gave sharp transitions with widths of the 
order of a few millidegrees. Different methods 
of melting were tried and possible contamination 
was checked for by comparison of the measured 
Tsc, but no significant differences were noted. 
The following procedure was finally adopted. 
About 0.35 gram of powder was pressed, out- 
gassed, and then melted in an argon arc furnace 
using a ruthenium cathode instead of the cus- 
tomary wolfram tip. The sample, resting ona 
water-cooled copper hearth, was melted in just 
one pass, lasting only seconds. The samples 
lost varying amounts of weight upon melting, 
ranging between 0.5 and 10%. 

The transitions were measured by immersing 
the samples in a bath of helium-3. The large 
temperature gradients which exist in helium-3 
made it necessary to mount a thermometer in 
good thermal contact with the sample. A non- 
encapsulated four-terminal germanium thermom- 
eter was used. It was calibrated against the 
vapor pressure of helium-3 under conditions 
where thermal gradients were avoided by stir- 
ring; thermomolecular pressure corrections 
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were made using the data of Roberts and 
Sydoriak.® 

The superconductivity was detected using the 
oscillator circuit of Schawlow and Devlin,’ in 
which the sample is in the center of the tank 
circuit coil. The natural frequency of the os- 
cillator circuit, 30 kc/sec, increased about 12 
cps when the sample became superconducting. 
This inductance change of the coil is due to the 
superconducting penetration depth being smaller 
than the skin depth which is about 10° A. This 
causes the small frequency change when the 
sample changes from normal to superconducting. 
There was, of course, a much larger change in 
the case of the powders. The oscillator current 
produced a magnetic field of about 0.02 gauss; 
the earth’s field was reduced by a mu- metal 
shield to about 0.003 gauss. 

The change in oscillator frequency as a function 
of temperature for pellets prepared from four 
samples of ruthenium powder is plotted in Fig. 1. 
Several warming and cooling curves were taken 
for each sample and the points are included. 

The frequency was read to the nearest cycle. 
The Johnson- Matthey powder is of spectroscopic 
grade in which detectable impurities are in the 
parts per million range or less. It is believed 
that chemical impurities are present in such 
small concentrations in all four samples that 
they have no measurable effect on Tg-. A sec- 
ond pellet of Ru®® (not shown in Fig. 1) was less 
pure. This sample had roughly twice the skin 
depth in the normal state and gave a broader 
transition curve about 0.004° below the others. 
Repeated arc-melting (three times) apparently 
boiled out the impurity somewhat and brought the 
transition temperature within 0.0015° of the 
others. 
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A shift for the Ru’™ curve of 0.012° below the 
Ru® curve, which would be expected if Tg, were 
proportional to the inverse square root of atomic 
mass, is clearly nonexistent. In fact, the Oak 
Ridge samples are within 0.001° of each other. 
From the agreement of these isotopes within 
themselves and with the natural ruthenium, we 
conclude that less than 10% of the normal isotope 
shift is present. This is considered as an indi- 
cation that mechanisms other than electron- 
phonon interaction may also lead to supercon- 
ductivity. Mechanisms such as s-d interactions, 
for instance, have been already recognized in 
the past for their importance to the supercon- 
ductivity of the transition elements.® 

We wish to thank P. P. Cioffi for measuring 
the residual earth’s field. We are grateful to 
J. P. Remeika for help in providing answers to 
questions of sample purity. We are particularly 
happy to acknowledge the work of L. D. Longinotti 
in the design and construction of the helium-3 
cryostat, the details of which will be published 


elsewhere. We are indebted to R. L. Garwin for 
sharing his knowledge of helium-3 with us and 
to P. W. Anderson, H. Fréhlich, C. Herring, 

H. W. Lewis, A. L. Schawlow, and P. A. Wolff 
for many discussions over the last six years of 
isotope effects in transition elements. 





1B. Serin, in Handbuch der Physik, edited by 
S. Fliigge (Springer-Verlag, Berlin, 1956), Vol. 15; 
J. C. Swihart, Phys. Rev. 116, 45 (1959). 

2B. B. Goodman, Nature 167, 111 (1951). 

3J. K. Hulm and B. B. Goodman, Phys. Rev. 106, 
659 (1957). 

“We are indebted to Dr. Bell of Lucius Pitkin, Inc., 
for carrying out the analysis. 

5We are indebted to Miss Koonce of our Laboratories 
for taking the electron microscope pictures. 

‘T. R. Roberts and S. G. Sydoriak, Phys. Rev. 102, 
304 (1956). 

"A, L. Schawlow and G. E. Devlin, Phys. Rev. 113, 
120 (1959). 

8H. Suhl and B. T. Matthias, Phys. Rev. 114, 977 
(1959); B. T. Matthias, H. Suhl, and E. Corenzwit, 
J. Phys. Chem. Solids 13, 156 (1960). 








EXCHANGE POLARIZATION AND THE MAGNETIC INTERACTIONS OF RARE EARTH IONS* 


R. E. Watson? 
Avco, RAD, Wilmington, Massachusetts 


and 


A. J. Freemant 
Materials Research Laboratory, Ordnance Materials Research Office, Watertown, Massachusetts 
(Received February 13, 1961) 


We are reporting results of a study of the con- 
tribution of spin or exchange polarization’ to the 
magnetic interaction of a rare earth ion with its 
neighbors and with its own conduction electrons. 
This investigation was carried out by means of 
conventional analytic Hartree- Fock (H- F) calcu- 
lations? for the Gd** and Gd* ions and a spin- 
polarized H-F calculation? for the Gd*® ion. 
While the conventional H-F results are of inter- 
est in themselves, since they represent the first 
time such information has been available for the 
rare earths, we shall only discuss a result of the 
Spin- polarized calculation which may have im- 
portant consequences for the magnetic behavior 
of these elements. Our results suggest that the 
rare earth ions carry a “paired” electron spin 
density which is negative in their outer reaches. 
Since the unfilled 4f shell electrons are imbedded 
in the interior regions of the ions, this outer 





spin distribution can play an important (and per- 
haps even dominant) role in the magnetic inter- 
actions of rare earth ions. In fact, as will be 
shown, these ions may even appear to their 
neighbors as having negative spins (i. e., anti- 
parallel to the 4f spin direction) and some ex- 
perimental results are discussed on this basis. 
The contribution of “paired” electrons to the 
magnetic properties of solids has been the object 
of recent studies*»* which have utilized the spin- 
polarized Hartree-Fock method. Two effects 
have been studied. The first of these is the con- 
tribution of the core s electrons, via the Fermi 
contact term, to the effective magnetic field at 
a nucleus in a ferromagnet.* This contribution 
appears to be the dominant source of the recently 
observed negative effective* fields in ferromag- 
nets. Secondly, computations’»* for the contribu- 
tion of the “paired” electron spin density® to an 
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ion’s neutron form factor suggest that the con- 
tribution is observable, a suggestion which has 
been borne out by a recent neutron diffraction 
investigation’ for NiO. This Letter discusses a 
third aspect of such spin polarization in solids, 
namely, its contribution to an ion’s magnetic 
interaction with neighboring ions and conduction 
electrons. The Gd*® ion’s unpaired 4f electrons, 
which are well in its interior, supply an extreme 
test of the effect. In what follows we presume 
that the free-ion results are a reasonable descrip- 
tion of the Gd “core” electrons in a solid, a not 
uncommon assumption. The neglect of important 
relativistic effects, limitations in the spin-polar- 
ized formalism,’ and limitations in computational 
accuracy for an ion of so high an atomic number 
make it advisable to view the results which follow 
for qualitative rather than detailed quantitative 
behavior. 

Figure 1 shows the computed “core” electron 
spin density (p, - p,) for all electrons other than 
the 4f shell and, for comparison, the 4f density 
as well. [Note the change of scales at r=3.0 
atomic units (a.u.) to a common scale for both 
Py-Py and the 4f density.] The two negative 
regions indicate densities associated with a spin 
antiparallel to the net spin of the ion. The region 
near the nucleus produces the negative effective 
fields of the type already discussed for iron. 

The outer region® is important for interactions 
with neighboring atoms and two examples are dis- 
cussed below. 

Nuclear magnetic resonance (NMR) studies® in 


magnetic crystals have revealed large internal 
magnetic fields at the nuclei of normally diamag- 
netic atoms like F~. The hyperfine interactions 
between the fluorine nucleus and the 3d magnetic 
electrons of the transition metal ion have been 
interpreted as arising from unpaired spins in the 
fluoride ion orbitals. This unpairing is thought 
to arise in either of two distinct ways: (1) be- 
cause of an admixture of covalent bonding into 
the purely ionic configuration,’® or (2) by the un- 
pairing action of the Pauli principle (“Pauli dis- 
tortion”)"' which affects those fluorine orbitals 
which have the same spin as the cation 3d orbital 
differently from those orbitals which have oppo- 
site spin. In either method the unpaired 2s elec- 
trons produce an isotropic hyperfine interaction 
whereas the 2p electrons produce an anisotropic 
interaction, the degree of unpairing being deter- 
mined by the squares of the overlap integrals, 
(S?), between the free F~ orbitals and the 3d elec- 
trons on the cation. Let us consider the applica- 
tion of this type of analysis to a rare earth salt 
such as GdF, in order to explore the consequences 
of the negative spin density in the outer region of 
the Gd** ion. Since the “paired” orbitals (5s and 
5p) have different radial distributions, they will 
overlap the fluoride orbitals differently; and it is 
this difference in their interactions which gives 
rise to a hyperfine interaction with the fluorine 
nucleus. Denoting the Gdt® 4f shell spin as up 
(t), then the extent of the interaction is meas- 
ured by S;?(t) -S;?(+) where i denotes some pairs 
of electrons. In Table I we list the squares of the 
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density. (See Note added 
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Table I. The squares of the overlap integrals, S?, 
between the Gd** outer orbitals and the F~ 2s orbital 
at the nearest neighbor distance of 4.4 a.u., and the 
net contribution to the internal field at the F~ nucleus 


(DS4?-27 Sy”). 








Sy? S;? 
4fy-2s  0.012x10-* 
5py-2s 0.365107? 5py-2s 0.483x10~ 
5sy-2s 0.085x107-? 5sy-2s 0.107x10~ 


27S4?=0.462 x10~? 2754? =0.590 x10~ 


Net effect = 2)S,?- 2)S,?=-0.128 x 107? (equivalent to 
-7200 gauss at the F~ nucleus). Contribution from 
the 4f shell alone=+ 0.012 x10~? (equivalent to + 700 
gauss). 





overlap integrals between the Gd** 4f(t), 5s(t 
and +), and 5p(¢ and +) and the F~ Hartree-Fock” 
2s orbital at the observed’ nearest neighbor dis- 
tance of 4.4 a.u. [For lack of space the aniso- 
tropic (2p) interaction will not be discussed here. | 
We see that the sum of S,’ is greater than the 
sum of S,” and that S,? -S,’ gives an effect which 
is ten times as large as, and opposite in sign to 
that obtained by considering the 4f overlap alone. 
In other words, in the context of either the Pauli 
distortion or covalency mechanisms (2 and 1 
above), our spin-polarized functions predict that 
a Gd** ion, as seen by a nearest neighbor F™~ ion, 
appears (in its interactions) to have a spin which 
is antiparallel to the actual Gd** spin." 

Jaccarino et al.’* have reported the NMR of Al?’ 
in the rare-earth intermetallic compounds and 
have interpreted the observed negative Knight 
shift as arising via the Ruderman- Kittel- Kasuya- 
Yosida’® mechanism from a negative exchange 
interaction between the localized 4f electrons and 
conduction electrons. We have calculated over- 
lap integrals and simple electrostatic exchange 
integrals (1/7,,) (which are important for direct 
magnetic interactions) between the Gd*® spin- 
polarized functions and Hartree- Fock Al wave 
functions” at the observed GdAlL, internuclear 
distance of 6.2 a.u. The outer region of negative 
spins again dominates over the 4f interaction for 
both the S*’s and the (1/r,,) integrals, i.e., the 
Gd** ion again behaves as if it had a negative 
Spin. In fact, these outer “core” electrons ac- 
count for about one-tenth of the observed Knight 
shift’ and their contribution is of the correct 
(i.e., negative) sign. In the metal the conduction 
electrons (which have not been included in our 
calculations) would also have a “paired” negative 
spin density in their outer regions; and since 








their overlap with the Al functions would be 
greater than that of the “core” functions, their 
contribution to the hyperfine field at the Al nu- 
cleus would also be correspondingly greater. It 
is, therefore, possible that the observed nega- 
tive Knight shift could be interpreted as arising 
from exchange polarization without invoking a 
negative exchange interaction between the local- 
ized 4f electrons and the conduction electrons. 

We also studied the exchange interaction be- 
tween a spin-polarized Gd*® ion and conduction 
electrons in the metal. The “conduction” elec- 
trons were assumed to be 6s electrons (which 
we obtained from a conventional and more accu- 
rate H-F calculation for the free Gd* ion) and 
exchange integrals were calculated between this 
6s electron and the spin-polarized Gd** orbitals. 
In this case the 4f exchange interaction (i.e., be- 
tween the 4/t and 6st electrons) was found to 
dominate over the net exchange interaction (i.e., 
the difference in exchange between the “core” t 
and + spins with the 6s electrons) which was nega- 
tive and about one-fifth that of the 4f electrons, 
in agreement with Hund’s rule for localized inter- 
actions. 

These investigations and related ones on the 
magnetic interactions of rare earth ions will be 
discussed at length in a forthcoming publication. 

We are grateful for the assistance of R. K. 
Nesbet, stimulating conversations with A. M. 
Clogston, V. Jaccarino, and M. Peter, and the 
help of A. Switendick with his computer programs. 


Note added in proof. Through an error the 5p spin 
density was not included in the (p, - py) plot shown in 
Fig. 1. The inclusion of the 5p spin density shows a 
much larger and more extensive negative spin density 
for large vy than that shown in the figure. 
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OPTICALLY DRIVEN SPIN PRECESSION 


William E. Bell and Arnold L. Bloom 
Instrument Division, Varian Associates, Palo Alto, California 
(Received February 14, 1961) 


The possibility of modulating a light beam at 
rf resonant frequencies by applying rf energy to 
an optically pumped spin system has been amply 
demonstrated theoretically’»? and experimental- 
ly.*»* We report here an interesting converse— 
that of modulating the light source at frequencies 
close to the Larmor frequency and observing a 
precessing spin polarization induced by the 
modulated light. The resulting resonance has 
quite different characteristics from the usual 
rf magnetic resonance; in particular, it does 
not show saturation effects although the line- 
width does depend in part on the average inten- 
sity of the incident light. This is, to our knowl- 
edge, the first time that an rf resonance effect 
has been produced by means that do not involve 
variation of a field directly coupling the rf- 
separated energy states. 

The effect may be observed in optically pump- 
able spin systems such as alkali metal vapors’»® 
or metastable helium’ by applying circularly 
polarized resonance radiation at right angles to 
the magnetic field H, (Fig. 1). In an earlier 
paper® we developed a phenomenological treat - 
ment of optically pumped spin systems that in- 
cluded the effects of several light beams. If we 
adhere to the notation of reference 3, and spe- 
cialize to a single light beam perpendicular to 
H,, and no rf field, we obtain the following equa- 
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tions for the spin system: 
dF /at +tyH \F +F/S, =P ml 


aM , /at +M_,/S, =0, 


where F=M,+iMy, Px is the pumping rate of 
the light, S, and S, are relaxation times including 
effects of the light, and 9m’ is the equilibrium 
polarization that would be produced by the light 
in the absence of a magnetic field and thermal 
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FIG. 1. Block diagram of apparatus. The inter- 
ference filter is employed in alkali vapor to remove 
the Sy.— Py, resonance radiation. In the metastable 
helium experiment no interference filter is used, and 
a high-frequency electrodeless discharge must be 
applied across the absorption cell to generate meta- 
stables. 
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relaxation effects. Normally (with no modula- 
tion) the steady-state solution is very small 
(Mz=0, F <1/yH,). However, let the light be 
modulated so that 


P =a+b coswt, 
x 


1/S,=R+b coswt, 
then a large steady-state solution exists of the 
form 
patho -iwt 


~ R+idw 
where Aw =yH, -w. 

The effect can be observed most easily as a 
change in average intensity of the transmitted 
light. The absorption of light by the cell is pro- 
portional to P,(1-M,), and the resulting signal 
in the photocell is then proportional to 


;RE? IM’ 
ig ee 
, ry (terms at w and 2w), 


which represents a secular increase in the trans- 
mitted light at resonance. 

One can picture the effect as follows. If one 
takes a “snapshot” of the system at ‘=O ina 
time short compared with a Larmor cycle, the 
effect of H, is a negligible perturbation and the 
light “pumps” a small amount of polarization in 
the x direction. This polarization immediately 
starts to precess in the xy plane at the Larmor 
frequency and is added vectorially to the polar- 
ization produced at all other instants of time. 
For unmodulated light this vector sum is clear- 
ly very small or zero. If the light is modulated 


at the Larmor frequency, however, then the 
polarization produced at ¢=0 is not exactly can- 
celled by that produced at ¢ =7/w and is rein- 
forced by the polarization produced at t=27/w, 
47/w, etc. An equivalent point of view is to con- 
sider the system in a coordinate system rotating 
at w; since the “rotating” lamp is turned on only 
when it is pointing in a certain preferred direc- 
tion, it is effectively stationary in a frame of 
reference for which the effective magnetic field 
is very weak. Optical pumping then takes place 
as if the effective field were parallel to the light 
beam. 

We have observed the effect in Cs and Rb vapor 
at Larmor frequencies between 175 and 350 kc/sec, 
and in He at 1.4 Mc/sec. Helium provides the 
best demonstration of the effect because a great- 
er index of modulation can apparently be achieved 
than with alkali vapor lamps, and the signal is 
proportional to the square of the modulation in- 
dex. In He with 60-cycle sweep and about 60% 
modulation of the lamp exciter carrier, the 
signal is comparable in amplitude to the stand- 
ard magnetic resonance signal. 

We are investigating the possibility of extend- 
ing this technique to higher resonant frequencies. 
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LEPTONIC DECAY OF A =~ HYPERON* 


Paolo Franzini! and Jack Steinberger? 
Columbia University, New York, New York 
(Received February 21, 1961) 


The search for the leptonic decay modes of the 
~ hyperons has not been entirely successful until 
now. The only previous positive evidence’ is the 
decay in flight of a = particle of unknown charge. 
If the decay particle is taken as a pion, the Q 
value would be too high by two and one-half times 
the quoted error. Recent systematic search for 
the leptonic decay in the bubble chamber among 
a@ much larger sample of decays has yielded no 
further confirmation of these decays.’ 


We report here a clear example of the decay 
=~ +e~+v+n. It was obtained in an exposure of 
the Columbia-Brookhaven 30-in. propane chamber 
to a 2-Bev 1~ beam at the Brookhaven Cosmo- 
tron. The event is reproduced in Fig. 1. The 
incident pion produces one long positive prong, 
a V°, and a prong which breaks sharply after 
7mm. Measurements on these tracks are given 
in Table I. There is no question that the track 
which is emitted almost straight back from the 
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FIG. 1. 1 +p 2~+K°+n"; B-+e-+n+v. Fora 
discussion of the various tracks see text. 


short prong is an electron. It leaves the chamber 
with a momentum of only about 25 Mev/c; all 
other known particles are densely ionizing and 
have, at most, 1 mm of range with such low 
momentum. The V° is a K° meson of 550+ 50 
Mev/c momentum. The positive prong with a 
momentum of 690+ 80 Mev/c is most probably 

a m* meson. No hyperon to carry the missing 
strangeness other than the short track is seen, 
and there is then the presumption that this track 
represents a decaying =~ meson. It is a priori 
fairly likely that the collision is a free-proton 
collision. About 40% of the collisions in pro- 


pane are of this type, and the carbon collisions 
normally exhibit evaporation prongs at such high 
energy. If the reaction is analyzed on this basis, 
1” +p-~X~+K°+n*, the unknown mass can be 
calculated from the measured K°, m*, and 1~ 
momenta. We find then, my =1160+50 Mev in 
good agreement with the known =~ mass.° The 
=~ hyperon has then a momentum of 860+ 100 
Mev/c. Transformed to the rest system of the 
=~, the electron energy is 163+ 20 Mev, well 
within the maximum energy of 257 Mev of the 

8 spectrum. 

The only rival interpretation on the basis of 
known interactions we have been able to think of 
is the following: A =~ is indeed produced, but 
suffers a collision: 


©” +p (free or bound) + A°+n + (n° or y). 


The A° is not observed; the 7° or y gives rise to 
a Dalitz pair whose positron escapes detection 
because its energy is less than ~ 100 kev or be- 
cause it annihilates in the initial millimeter of 
flight. Such an event can be expected about once 
for every 10° K°-° production events. We feel 
that this is small enough compared to the number 
of observed K°-=~ productions to be ignored. 
Another possible explanation, the reaction 


K~ +p (free or in C) +~K*+K°+K~ +n, 


produced by the K contamination of the beam and 
followed by Kz3° decay in flight, is not possible 
on energetic grounds. 

The only reasonable interpretation of the event 
seems, therefore, that of the 8 decay of the =". 








Table I. Measured and calculated angles and momenta. 
Particle 
number Interpretation Dip angle Azimuth Momenta (Mev/c) 
1 Incoming 1 0° 0° 2000+ 50 
3 i. -7° +7° (meas) -5.5° +2° (meas) P 
-4° (calc) -7° (calc) 860 +100 (calc) 
3 K® -33.6°+4° (meas) 1° (meas) 
-29° (calc)® 550+ 50 (calc) 
4 a 26° +3° - «-* 690+ 80 
5 e~ -29.5°+3° 155° +1° 87+ 10 
6 a* from K° decay -2.5°+3° -20° +1° 500 +150 (meas) 
375 +. 30 (calc)® 
7 a from K° decay -54° +3° 47” +1° 320 +100 (meas) 


315+ 30 (calc)* 





“Fitting the K°-decay kinematics. 
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The absolute value of the rate of this new proc- 
ess is important for the understanding of the 
Fermi couplings of hyperons.* Clearly, an iso- 
lated event does not permit an answer to this 
question. We are pursuing a systematic analysis 
of our pictures towards a determination of the 
absolute 5* g-decay rates. 

It is a pleasure to acknowledge the generous 
help of Dr. Albert Prodell who was instrumental 
in the design of this chamber and has brought it 
into operation. 
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K’*+=~, From the K° and ° momenta, we find mg,c* 
= 933 +33 Mev for this event. 

‘See, for instance, R. P. Feynman and M. Gell-Mann, 
Phys. Rev. 109, 193 (1958); L. Okun’, Annual Review 
of Nuclear Science (Annual Reviews, Inc., Palo Alto, 
California, 1959), Vol. 9, p. 61. 








PROPERTIES OF THE Y* AS OBSERVED IN THE INTERACTION K,°+p + A°+2*+7°* 
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and 
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(Received February 14, 1961) 


An investigation by Alston et al.’ of the re- 
action K~ +p~-A°+n*+77, at a K~ momentum 
of 1120 Mev/c, has shown that strong momen- 
tum correlations exist among the reaction pro- 
ducts. They concluded that these correlations 
were due to the existence of a strangeness -1, 
isotopic spin 1, baryon state, Y ae which decays 
by strong interactions into a A°+z. Analyses of 
their data indicated to them that the mass of this 
state was 1384 Mev and the width was 63 Mev. 

We are reporting on similar analyses of the 


products of the related reaction K,°+p—~ A°+2*+7°. 


The K,° mesons were produced by a chain of re- 
actions initiated by the external beam of the 
Brookhaven Cosmotron. Three-Bev protons 

from the external beam were focused on a steel 
target to produce 7~ mesons which in turn were 
momentum-analyzed and focused on a polyethylene 
target. The K,° mesons produced in the forward 
direction by reactions such as 7~ +(p or n) + 

(A° or £)+K,°+(0,1, or 2)”, passed through a 
lead absorber and a magnetic clearing field 
designed to eliminate charged particles and y 
rays, and into a 14-in., 20-liter, liquid hydro- 
gen bubble chamber 10 ft from the target. An 





average flux of 710° protons per pulse incident 
upon the meson target produced an intensity of 
about 3 x10° 7~ mesons with a momentum of 1.60 
Bev/c and a spread of 7%. At these z fluxes, 

we estimate, from consideration of the number 
of K,° decays observed and the measured life- 
time, that the K,° flux through the chamber was 
about 1.5 per pulse. Maximum intensities about 
three times these values were achieved. For 
each event of the type K,°+p—A°+7°+m*, meas- 
urements of the momentum of the associated 1* 
and the A° produced in the chamber, together 
with the knowledge of the direction of the incident 
K,°, allowed the determination of the momentum 
of the 7° and the K,°. This momentum distribu- 
tion for the 60 analyzable events reported here 
can be described adequately as a Gaussian dis- 
tribution centered at 975 Mev/c with a half-width 
of 100 Mev/c. From charged = and 3-prong 
events we estimate that about 15% of these events 
are £°+2m and A°+3n events. 

Any interaction of a K,° with a proton which 
produces a hyperon and 7 mesons must be initi- 
ated by the strangeness -1, K° part of the K,"°. 
Since the third component of the isotopic spin of 
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the K° is +4, the interaction with the proton takes 
place wholly in the T=1 state. As the A° has an 
isotopic spin of zero, the two 7 mesons, of the 
reaction products A°+7*+7°, will have isotopic 
spin one, and the final space wave function will 


only change sign upon exchange of the two mesons. 


Therefore, if the interaction is charge independ- 
ent, any description of the distribution of reac- 
tion products must not depend upon the charge 

of the 7 meson. Conversely, an observed sym- 
metry of this kind is evidence for charge inde- 
pendence in interactions involving strange parti- 
cles. 

When the observation of final-state interactions 
is important, it is useful to consider the dis- 
tribution of Q values as calculated between pairs 
of particles. This is illustrated in Fig. 1, where 
the distribution of Q values is plotted for (A°, 7*), 
(A°, 7°), and for the lesser of the two values for 
each event. Independent of any detailed assump- 
tions concerning the origin of the peaks near 
130 Mev, the similarity of the (A°, 7°) and (A°, 7*) 
distributions is evidence of the charge independ- 
ence of the interaction. This sharp Q-value 
peak suggests again the existence of the Y* state 
previously postulated. About 40 events can be 
assigned to the peak, 22 Y* and 18 Y°, the Y° 
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FIG. 1. Histograms representing Q-value distribu- 
tions. The lower histogram shows the distribution of 
@ values calculated for the A° and 2°, the center histo- 
gram represents the same plot for A°+n°, and the 
upper histogram shows a plot of the lower of the two 
@ values plotted for each event. Sixty events are 
represented. 
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not having previously been observed. While 
alternative explanations of this anomaly can be 
offered,” we shall proceed in the rest of this 
discussion on the assumption that the explanation 
in terms of the existence of a state, y”. is cor- 
rect. 

The experimental width of the state, as seen 
in Fig. 1, is about 29 Mev. This measured width 
will be approximately equal to the square root of 
the sum of the squares of the natural width and 
the experimental resolution. From a study of 
our errors we conclude that the full width of the 
resolution function cannot be less than 20 Mev 
and is more likely near 30 Mev. The value of 
the natural width is then, 0<I <20 Mev, which 
is quite different from the value of 63 Mev quoted 
by Alston et al.’ 

Information concerning the spin and parity of 
the state may be derived from analysis of angular 
distributions of Y* and A° decays. The Y* pro- 
duced at small angles with respect to the K beam 
will decay preferentially in the beam direction® 
if the spin is greater than 3. Of the 17 events 
such that |cos@|>0.7, where @ is the production 
angle, the fraction 0.50+ 0.15 decayed such that 
|cos¢|>0.5, where ¢ was the decay angle, in 
agreement with the value 0.5 to be expected for 
J=4, and about 1.3 standard deviations from the 
value 0.69 for J= 4. For any angle of production 
there will be some alignment, if the spin is 
greater than 4, and the direction of the align- 
ment will vary slowly with angle. Therefore, 
the examination of correlations between the direc- 
tions of decay of Y* states produced at nearly 
the same production angle will be informative.’ 
The quantity C =2;2;P2(cos6;;) was evaluated, 
where P, is the Legendre polynomial normalized 
over the interval 0 to 1, and 64; is the angle be- 
tween the decay of the Y* event i and event j. The 
sum i is over all events, the sum j over all 
events such that |cos¢; -cos¢;| <0.3, in this 
case, where ¢ is the production angle of the Y*. 
The value of C from this experiment was -0.03 
+ 0.10 which is to be compared to expected values 
of 0 for j=4, and 0.09 for j= 3. 

Information concerning the parity of the Y* can 
be obtained if the Y* is produced such that it is 
polarized. Assuming that the spin is 3, if the 
Y* decays by S wave the A will be polarized in 
a direction perpendicular to the plane of produc- 
tion and will decay with an angular distribution 
of 1+A,cos@,, where @, is the angle between the 
decay direction and the normal to the plane of 
production (Py x Py)/(|Px | |Py|) and As = Pa; 
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P is the degree of polarization of the Y* and a 

is the usual asymmetry parameter in the decay 
of the A which is about +1. If the decay is by P 
wave, the distribution will take the similar form 
1+Ap cos@p, where Ap is again Pa but 6» is the 
angle between the A-decay direction and a vector 
which lies in the plane formed by the production- 
plane normal and the Y*-decay direction, and 
whose angle with the normal is twice that of the 
Y* decay. The measured values of the quantities 
Ap and A, are correlated; for an S decay, Ap 
=-3Ag, while for a P decay, Ag=-3Ap. Our re- 
sults, using 37 events, are A,=-0.38+0.25, Ay 
=0.19+ 0.25, which suggests that the decay is 
S-wave. Statistically more significant results, 
in accord with these, have been reported by 

] Ferro-Luzzi et al.* 

If these spin and angular momentum assign- 
ments are valid, and if the K-meson parity is 
odd, this state can be expected to be the state 
related to the (a-) K~-p scattering length solu- 
tion of Dalitz and Tuan,° a = (-1.0+0.18i) fermi. 
It is then interesting to calculate the threshold 
behavior of the K~-p system under the assump- 
tions that a single-level resonance formula® 
dominates the situation, and the -A parity is 
odd. The scattering length will then be equal to 
yk’ /(E, - #T)-R, where I; =2kiy;’ is the width 
for channel i for S wave and I; =2ky*(kR)?/[1+(kR)’] 
for P wave, I is the total width equal to [py +I, 
+T's, where K, A, and = represent the A+z, 
o+m, and K~ +p channels, respectively, E, is 
the Y* resonance energy, and R is an interaction 
radius, taken here as h/Myc. From the Y* 
widths taken as I’, =16 Mev and I's =4 Mev,’ 
we find y,*=7 Mev f, and y>*=45 Mevf. We 











choose yx” =30 Mev f to fit the real part of the 
scattering length. With these numbers, a 
= (-1.0+0.322) f and the ratio of = to A produc- 
tion at threshold in the T =1 state is about 0.75, 
not greatly different than results of the analysis 
of Dalitz and Tuan.® Even A-= parity results 
in a much poorer fit to the 2-A production ratio. 
It is interesting to note that the sum of the 
reduced widths, about 80 Mev f, is, within the 
very rough estimates made here, consistent with 
the Wigner sum rule,® Dgys* sf? /MR =200 Mev f, 
where the sum is over all channels open or closed. 





*Work performed under the auspices of the U. S. 
Atomic Energy Commission. 
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ton, Indiana. 

ton leave at the University of California, Berkeley, 
California. 

Now at Centre d’Etudes Nucléaires de Saclay, Saclay 
(Seine et Oise), France. 

'M. Alston, L. W. Alvarez, P. Eberhard, M. L. 
Good, W. Graziano, H. K. Ticho, and S. Wojcicki, 
Phys. Rev. Letters 5, 520 (1960). 

2R. K. Adair, Conference on Strong Interactions, 
Berkeley, 1960 [Revs. Modern Phys. (to be published)]. 

3R. K. Adair, Phys. Rev. 100, 1540 (1955). 

‘mM. Ferro-Luzzi, Conference on Strong Interactions, 
Berkeley, 1960 [Revs. Modern Phys. (to be published)]. 
Dr. Ferro-Luzzi and his collaborators independently 
deduced the method of determining parity presented 
here. 

°R. Dalitz and S. Tuan, Ann. Phys. 10, 307 (1960). 

SE.g., E. P. Wigner and L. Eisenbud, Phys. Rev. 
72, 29 (1947). 

Preliminary investigations of the +27 production 
from this experiment, by B. Musgrave, suggest that 
T's/T, may be about 1/4. 

SE. P. Wigner, Am. J. Phys. 17, 99 (1949). 








VoLuME 6, NUMBER 6 


PHYSICAL REVIEW LETTERS 


Marcu 15, 1961 





ELECTROMAGNETIC PROPERTIES OF THE PROTON AND NEUTRON* 


D. N. Olson, H. F. Schopper,? and R. R. Wilson 
Laboratory of Nuclear Studies, Cornell University, Ithaca, New York 
(Received February 15, 1961) 


Measurements of electron-proton scattering 
made at Cornell’ and Stanford? below 1000 Mev 
have shown that Fj, and F 2p» the electric and 
magnetic form factors of the proton, are not 
equal; indeed that F 2p is nearly zero for values 
of g, the momentum-energy transfer, of about 
5 fermi~’, while F1, remains in the vicinity of 
0.4—all this is indicative of a core of charge of 
about 0.4 e. 

New measurements of scattering cross sec- 
tions from hydrogen and deuterium at higher 
values of g appear to permit us to determine the 
extent and charge of this core of the nucleon and 
also to deduce similar properties of the sur- 
rounding mesonic cloud. The Cornell 1.3-Bev 
electron synchrotron has been used in extending 
our earlier measurements in exactly the manner 
described in Nature.’ The new results are shown 
in Table I and Figs. 1 and 3. It is of striking 
physical significance that the cross section falls 
so low at 1.2 Bev, i.e., to 107 cm?/sr. 

Now one of the simplest models of the nucleon 
has a point core of charge +$e surrounded by an 
extended mesonic cloud of +$e for the proton or 
-4e for the neutron. Indeed, some aspects of 
such a model are reflected in the experimental 
data, but the model departs significantly from 


Table I. 
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FIG. 1. Experimental results for electron-proton 


scattering at 112°. 


the facts at a number of points. In the first place, 
from Foldy’s interpretation of neutron scatter- 
ing by atomic electrons® we know that the rms 
radius of the charge distribution of the neutron 

is zero and this implies that at q=0, Fj, £0.002q’; 
the above picture implies that Fj, =3(0.8f¢, ¢ 
being measured in fermi~. In the second place, 


Experimental results. Columns 2 and 3 give the experimentally determined cross sections and peak 


cross sections for electrons of incident energy E, scattered through an angle of 112°. Radiative corrections have 
been applied. In order to calculate the integrated deuteron cross sections, allowance must be made for electrons 
which do not scatter into the detection channel because of internal motion of the nucleons in the deuteron. Cal- 
culating the energy spread of scattered electrons using the impulse approximation and including the resolution of 
the apparatus as determined from electron-proton scattering, we obtain the values for AE/E listed in column 4. 
Then, subtracting the proton cross section from column 2, we arrive at the integrated cross section for electron- 


neutron scattering given in the last column. 














Ey do/dw) (d*a/dE dw) 4 (do/dw) , (do/dw) 
(Mev) (107°? cm?/sr) (107*3 cm?/Mev sr) AE/E (107? cm?/sr) (107%? cm?/sr) 
200 6.6 +0.6 2.30 +0.33 0.154 5.35 +0.77 -1.2 +1.2 
325 2.40 +0.25 0.85 +0.076 0.133 2.49 +0.23 0.09 +0.34 
502 0.76 +0.06 0.250 +0.016 0.112 0.78 +0.11 0.15 +0.12 

640 0.270 +0.023 
761 0.119 +0.008 0.034 +0.008 0.093 0.114 +0. 027 0.005 +0. 030 
900 0.059 +0.006 
1000 0.037 +0.005 0.0158 +0. 0028 0.080 0.051 +0. 009 0.014 +0.010 
1100 0.023 +0.004 
1200 0.0123 +0. 0025 0.0059 +0. 0028 0.074 0.019 +0. 009 0.007 +0. 009 
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the form factor F;, should approach 0.5 for 
large values of q*; in fact, it seems to be falling 
below 0.4. Nevertheless, we have interpreted 
our experiment in the spirit of this very simple 
model; to fit the facts we have found it necessary 
to assign a radius of about 0.2 f to the core and 
to put a small part of the positive charge of the 
core in an extensive cloud that is the same for 
neutron and proton. 

If we assume that Rosenbluth’s formula is still 
valid at our energies, i.e., 


do_ et 


_ e* cos?(6/2) 2E\ in 
To seems  **(Gz) era 
x [+ gles AP, Ptan(e/2) +P. , (1) 


where E is the incident electron energy, @ the 
laboratory angle, M the nucleon mass, and k the 
anomalous magnetic moment, then we want to 
find the form factors F, and F, which are func- 
tions of g*, and for this at least two measure- 
ments of do/dw at the same g value but at differ- 
ent angles are necessary. The data have been 
resolved in this manner for values of q*< 25 f~? 
where cross sections have been measured at 112° 
and 66°, and the values of F, and F, that have 
been obtained are plotted in Fig. 2. At values of 
¢ >25 we have only the data at 112°; at smaller 
angles the large values of q* cannot be reached 
with our available energy nor can we easily go 

to larger angles for then the counting rate be- 
comes prohibitively small. We are preparing 
new experimental equipment, however, to do 

just this. 

On the other hand, by measuring electron- 
deuteron scattering, we have been able to infer 
values of F, and F, for both the proton and neu- 
tron as follows: Equation (1) can be rewritten 
in the form 


do/dw @.4F; +41 oF Fo +4a9F 5 os? (2) 


so that we may use the values of the constants 
and o,.., the point “no-spin” electric cross sec- 
tion, which have been conveniently tabulated.‘ 
Thus an ellipse is obtained for a particular value 
of the cross section when F, is plotted against 
F,. Such an ellipse is shown in Fig. 2 for the 
highest energy measured (g*=37). Depending on 
the value of F,, F, can take on values between 
plus and minus 0.3. Thus it may be that our 
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FIG. 2. (a) Experimentally determined form factors 
for the proton, as functions of the momentum transfer 
q. Circles (@) are used for the form factor F,, and 
squares (m) for the magnetic form factor Fo as de- 
termined from measurements at two different scatter- 
ing angles (points with flags: Cornell data; points 
without flags: Stanford data). For g?>25, only meas- 
urements at 112° were performed (Fig. 1). These 
determine ellipses in the F,-F, plane, as shown in the 
insertion, from which the extreme values of Fip and 
Fo» can be inferred. These are shown by the arrows 
but are not, of course, simultaneously realized. As- 
suming that Fo» is given by curves I, Il, or III, one 
can derive respective values for F,,. The results 
presented in Fig. 3 exclude curves II and III, so that 
curve I (4) gives the proper values for F,, and Fop- 
(b) F; can be decomposed into partial form factors 
according to Eqs. (3) and (4). Three fits using the 
core charges and radii given in the insertion are 
shown. The experimental data are best reproduced 
by case (b). Only for this case are the partial form 
factors displayed. Fj, computed by using them is also 
included in the figure. 
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measured cross section, ~10~* cm?/sr, is 
small because of destructive interference be- 
tween large values of a positive Fj, and a nega- 
tive F9,. In the case of the neutron, however, 
the magnetic form factor would be just reversed 
so that the interference should then be construc- 
tive and a large value of the cross section would 
result. That the magnetic form factors of the 
proton and neutron are very nearly equal and 
opposite follows from the fact that the anomalous 
magnetic moments are equal and opposite to 
within 5%.° The curves of Fig. 3 show the ratio 
of neutron and proton cross sections based on 
calculations using expression (2) for the neutron 
cross section with various assumptions for F,, 
but assuming that F, is due to a core and is the 
same for neutron and proton-—essentially true 
for g’>25. The experimental points, obtained by 
replacing our CH, targets with CD,, using the 
peak value method,* and then making the appro- 
priate subtractions, are also shown. They are 
consistent with curves having values of F9 
=0+0.03. We conclude, therefore, that F2,~0 
and is, in fact, given as a function of ¢ by Curve 
lof Fig. 2. This can then be used to obtain the 
corresponding values of F 1, from ellipses simi- 
lar to that of Fig. 2 but corresponding to meas- 
urements at ¢* = 28.7 and 32.7. These are plotted 
in Fig. 2. Curve II of Fig. 3 and the correspond- 
ing curve of Fig. 2(a) indicate the accuracy of 








6r pA 
/ 
/ 
5+ / 
/ 
al / 
On / 
r / 
Pp 3+ / 
/ 
f 
2 
l J a 
ir / A -—t 
a a ee 
Vv 
5 * 10 20 30 aot? 
2 
q==.8 


FIG. 3. The ratio of the neutron and proton cross 
sections as a function of the momentum transfer q. 
Curves I, II, and II were calculated with form factors 
Fip, and F2p =F 2 given by the curves I, I, and I] 
in Fig. 2(a), and with Fj, from Fig. 2(b). The experi- 
mental points were computed from the data in Table I. 
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this procedure, i.e., about + 0.03 in the value of 
F2p and + 0.05 in Fy. Note the pronounced min- 
imum in the curve of electron-neutron scattering 
at g’=15f7*, i.e., about 700 Mev. 

Now let us see how the simple core model must 
be changed to conform with these form factors. 
We will separate each form factor into an iso- 
scalar and an isovector partial form factor—the 
vector part changing its sign when we change 
from a proton to a neutron, the scalar part re- 
maining unchanged. We can write 


Pip Fistfiv Fin**is-F iv (3) 


+F F, =F,.-F 


Pop FostFaw Fan=FosFoy = @) 


As the charges of the proton and the neutron are 
different, F, must consist of an isoscalar and 
an isovector part. F,, on the other hand, will 
be determined predominantly by Fay as the 
anomalous magnetic moments of proton and neu- 
tron agree within 5%. We believe that the ex- 
periments are not accurate enough to allow the 
detection of an F9g contribution of the order of 
a few percent and therefore we assume Fo5=0.° 
Then in order to get a vanishing rms radius for 
the neutron charge distribution, F, has to be de- 
composed in at least three terms. We shall 
write Fys5 as a sum of two terms and shall asso- 
ciate F;>°, the term with the smaller radius, 
with a core, the other one Fy," with a cloud. 

It will be shown that all experimental results 
can be reproduced without assuming that Fy y 
has a core, too. Again, however, we cannot 
definitely exclude such a core because of exper- 
imental uncertainties. With each form factor 
we can associate a charge, i.e., the value of the 
form factor at q=0, and for Fy s° Fis", and 
Fyy we have, respectively, e,, e wp and éy. 
The charge on the proton requires that e+e m 
+ey=1, and that of the neutron requires that 

eo +e, -ey=0; from this it follows that e,+e 

= 4 and that ey=4, a unit charge being that of 

the electron, of course. We can also write some 
relations for the mean square radii associated 
with each of the above charge distributions, i.e., 


te. 2 2 Ps 

a, ef. +¢ a, +e yay, (0.8 ff, (5) 
and “ 3 2. 
a. re ate a, -e a," =0, (6) 


where a, and a, are the measured rms radii of 
charge already determined at low q’. It follows 
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that 


HH 


and hence that ay =ap. 

The core model implies that a,”?«<a,,* in which 
case (7) gives a *=ay' /2e . Weare left with 
the problem of determining e,, ande,. This we 
do by fitting our experimental curve of Fj, at 
very large values of ¢’, i.e., between 20 and 37, 
at which place both Fy5" and Fy have become 
small compared to F;>°. For example, let us 
assume most simply that the core is a point 
charge; then, quite independently of what we 
choose for the form of Fs or Fy, we obtain 
e,=0.25 and it follows that ay, =1.1. However, 
such an assumption does not give a good fit of 
the experimental data. See curve (a) of Fig. 2(b). 
We can do better if we assume an extended dis- 
tribution for the core charge-—for simplicity a 
Gaussian form is used. Now although we know 
that ay =ap= 0.8 f, we do not know the shape of 
the charge distribution. We assume rather 
arbitrarily the same Gaussian shape that fits 
Foy. This becomes quite negligible for g>25 f~? 
where only the form factor of the core remains 
and which we then fit in this region. We could 
arbitrarily make e,=0.5; then a, would be 0.37 f, 
but the fit at intermediate values of q* would be 
poor. The best fit gives e,=0.35+ 0.1 and a, 
=0.2+0.1f. This choice determines that e Mm 
=0.15+0.1 and a, =1.4+0.4f. (See Table II.) 
The corresponding partial form factors are plot- 
ted in Fig. 2(b) together with the total form fac- 
tor for the neutron and proton given by (3). It 
should be emphasized that the model leaves only 
one free parameter, i.e., @-, if a point core is 
assumed, and only two parameters for an extended 
core. The shapes of the distributions can also be 
considered as free parameters, but we have found 
by trial that our results are insensitive to whether 
Gaussian or exponential forms are used. Yukawa- 
type distributions give somewhat different results, 


ea, +e a “= ta," (7) 


Table Il. Parameters giving the best fit to all ex- 
perimental data with F;=e; exp(-a;7q?/6). 





€- = 0.35 40.1 a, =0.2+0.1f 
ey =0.15 20.1 a, =1.4440.5f 
€yy=0.5 a,y=0.8 f 
egy=1.0 @oy=0.8f 








but such distributions have a built-in core. How- 
ever, the spirit of our approach and underlying 
it has been the core model; hence, we have 
avoided, but by no means dismissed, Yukawa 
distributions.® 

The core form factor and its range seem quite 
reasonable, but the reason for the long-range 
isoscalar charge form factor, F; gs, is not so 
obvious. If F;y is due to two-meson processes, 
then F;¢ may be due to three-meson processes® 
but why, except for a very strong resonance, 
should this have a longer range? Again, it may 
have its origin in some obscure one-meson proc- 
ess. Is it possible that F;¢ is of electromag- 
netic origin, i.e., is the meson cloud polarized 
by the core charge? The sign of the effect is 
right, although the Coulomb forces would seem 
to be too weak. 

Taking the Fourier transforms of the form 
factors, one can compute the charge distribu- 
tions of the proton and neutron. These are 
plotted in Fig. 4. Of course, these distributions 
and the form factors themselves will be changed 
by relativistic effects which have not been taken 
into account.” Also, by choosing somewhat dif- 
ferent models one might obtain somewhat differ - 
ent parameters. We believe, however, that the 
charge distributions shown in Fig. 4 give a quali- 
tative picture of the nucleons. Its main features 
are a positive core for both particles with a 
charge of about 30% and a radius of approxi- 
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FIG. 4. Charge distribution for the proton and the 
neutron implied by the form factors shown for the 
fit (b) in Fig. 2(b). 
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mately the nucleon Compton wavelength. The 
proton core is surrounded by a positive cloud, 

the neutron by a negative one. The neutron has 

in addition a positive shell at its outside that con- 
tains a few percent of the elementary charge. 

The distributions of the anomalous magnetic mo- 
ments are spread out with rms radii of about 

0.8 f and it is not necessary, with present ac- 
curacy, to assign a magnetic core.® 
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DIRAC AND PAULI FORM FACTORS OF THE NEUTRON* 
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Recent work on electron scattering at Stan- 
ford'~* has shown that the electromagnetic form 
factors of the proton were split apart at large 
values of the momentum transfer (q) and the de- 
tailed behavior of the Dirac 1p) and Pauli F 2p) 
form factors was reported. These studies showed 
also that F2p is approaching zero and that the 
electron-proton scattering cross section exhibits 
a diffraction dip at ¢ =25 f~* which is associated 
with the behavior of F2, at that value of the mo- 
mentum transfer. Some information concerning 
the proton form factors has also been reported 
by the Cornell group.* The information in refer- 
ences 1-3 was used by Herman and Hofstadter,° 
who deduced values of the Dirac and Pauli form 
factors of the neutron (Fj, and F9,, respective- 
ly) from the above data by making the assump- 
tion that F9,,=F2, which was known from earlier 
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measurements®’’ to be roughly true at low values 
of g’. In this way the work of reference 5 showed 
that F;,,#0. Although there is an ambiguity in 
the sign of Fj,, Herman and Hofstadter chose 
the negative sign because it has been commonly 
accepted that the charge cloud of the neutron is 
due primarily to the presence of negative mesons. 
The chief result of the present communication is 
the independent experimental determination of 
the two form factors of the neutron (Fj,, Fay) 
and a verification that F},,#0. In another com- 
munication® we attempt to resolve the ambiguity 
of sign in Fin- 

The above results were obtained by combining 
measurements of the inelastic electron scatter- 
ing cross section of the deuteron at two sets of 
values of energy (E) and angle (6) of the scattered 
electron for the same value of g*. In essence 
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this method is the same as the method of inter- 
secting ellipses used in determining Fp and 

F 2p of the proton.*?° The application of this 
method to the neutron was given in reference 9. 
This method of determining Fj, and F2, elimin- 
ates many errors. In the present work the modi- 
fied Jankus theory of electrodisintegration of the 
deuteron was used to evaluate Fj, and F2y from 
the value of the cross section at the peak of the 
inelastic continuum.*®*"*»!* The modified Jankus 
theory was employed in an extended form pro- 
vided by Goldberg"* which takes account of finite 
nonzero values of Fj,. Calculations made by 
Durand" show that the modified Jankus theory’® 
is quite accurate at the peak of the inelastic con- 
tinuum. 

In all cases the measurement of the deuteron 
peak was accompanied by a corresponding meas- 
urement of the cross section of the proton peak. 
This procedure minimizes many possible exper- 
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FIG. 1. A pair of deuteron inelastic continua is 
shown for g?©16.7 f~?, along with the comparison 
proton peaks used for absolute calibration. The very 
small x* background at 858 Mev and 75° was meas- 
ured in order to estimate the 1 background under the 
deuteron peak (see references 7 and 11). 


imental errors. Thus, an absolute cross section 
of the deuteron peak could be obtained by using 
the absolute data in references 1-3. Radiative 
corrections were calculated for the deuteron and 
proton peaks by Sobottka’s method’ and the mag- 
nitude of the radiative correction applied in find- 
ing the final value of the deuteron cross section 
was nearly always small and constant: ~10%. 
Furthermore, the differences in the corrections 
for the interaction in the final state’* were esti- 
mated and found to be <~4% when calculated for 
the two members of the deuteron cross-section 
pair. The final-state corrections were not ap- 
plied in the above evaluation of the deuteron 
cross sections since they are small (<~4%) and 
not particularly well known at the present time. 
Improved calculations of these interactions are 
now in progress.*® 

In Fig. 1 we show a pair of deuteron inelastic 
peaks at a value of g* =16.7 f~? with the accom- 
panying proton peaks. The data were taken with 
targets of liquid hydrogen and liquid deuterium 
with the new 72-in. magnetic spectrometer.’ 
We present in Table I the values of the deuteron 


Table I. Experimental electron-deuteron scattering 
cross sections. 








(70 /d24E) max*10™ 
gq’ (f-2) E,(Mev) 0° (cm*/sr Mev) 
5.1 300 135 8.32* 
500 60 35.9* 
7.0 362 135 4.6** 
600 60 17.8* 
8.6 414 135 3.07** 
675 60 11.1 
9.1 429 135 2.65** 
700 60 8.9 
11.5 500 136 1.59 
800 60 5.44 
14.7 590 135 0.82 
850 67.5 2.6 
16.7 646 135 0.51 
858 75 1.31 
18.0 675 135 0.46 
900 75 1.32 
21.0 750 141.5 0.25 
900 90 0.53 





®The experimental cross-section values with one 
asterisk have been taken from Sobottka, reference 7. 
Those with two asterisks have been interpolated by 
using Sobottka’s values for 6 = 135°. 
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FIG. 2. An example of intersecting ellipses accord- 
ing to the modified Jankus theory for the pair 858 Mev, 
75° and 646 Mev, 135° at g?=16.7 f*. Note the effect 
that would be caused by +10% experimental errors in 
the deuteron cross section. The value of the cross 
section at 646 Mev and 135° used in the text is slight- 
ly different from the value used in this example and 
is 5.1x 10~* cm?/(sr Mev). 


cross sections obtained in the above manner. 

The convenient method of intersecting ellipses 
used in obtaining the form factors of the neutron 
is illustrated in Fig. 2 at g=16.7f~*. The effect 
of possible errors in the measurement of the 
deuteron cross sections is also shown in the 
figure. The errors in the neutron form factors 
arising from all sources other than possible sys- 
tematic errors are believed to be approximately 
of the order of the spread in the final values. 

The errors in the deuteron cross section are 
believed to be less than about + 10%. 

In studying the neutron problem, one finds 
four possible intersections of ellipses for a given 
value of g which determine the Fj, F 2, pair. 

Of these we have chosen the one set which seems 


to be physically reasonable. This set is given in 
Table II together with the known proton form 
factors.’~* Further details concerning the other 
possible solutions will be described subsequent - 
ly. It should be noted that in all cases F2,#F 25 
although at small values of q’ the assumption of 
equality used in reference 5 is satisfied approxi- 
mately. 

From the chosen set of values of Fj», Fiy, 
F 2p» Fn, we may now form the corresponding 
set of values of the isotopic form factors which 
are defined as follows: 


Pig hy TF in? (1) 
Piy ip Fin’ (2) 
Fo,=[1.79F, -1.91F, ]/(-0.12), (3) 
Foy =[1-19F,, +1.91F,, J/(8.10). (4) 


The experimental values of the isotopic form 
factors as defined above are given in Table II. 
We observe that F;y and Fay are quite similar 
to each other while Fs has a different behavior. 
Fs is the least well-known isotopic form factor 
since it is associated with the small difference 
between the absolute values of F2, and F2y. We 
are investigating the possible errors in the values 
of the isotopic form factors but believe that the 
errors are not large enough to change the basic 
pattern of behavior exhibited in Table I. 

We wish to thank C. Buchanan, F. W. Bunker, 
B. Chambers, H. Collard, and E. Dally for their 
help in this experiment. Our special apprecia- 
tion is due Donald W. Aitken who helped us con- 
siderably during the early stages of this work. 
Finally we wish to acknowledge a number of 
valuable discussions with Dr. A. Goldberg and 


Table Il. Proton, neutron, and isotopic form factors. * 





(1) (2) (3) 
F F F 


(4) (5) (6) (7) (8) 
F F F F F 





q? (£7?) 1p 2p 1n 2n 1S 2s 1V 2Vv 
0 1.00 1.00 0.00 1.00 1.00 1.00 1.00 1.00 

5 0.62 0.62 0. 04 0.70 0.66 1.9 0.58 0.66 

10 0.48 0.39 0.11 0.51 0.59 2.2 0.37 0.45 
15 0.43 0.21 0.17 0.38 0.60 3.0 0.26 0.30 
20 0.42 0.08 0.19 0.32 0.61 3.9 0.23 0.20 





“The entries in columns (3)-(4) represent smoothed values taken from curves through the experimental data of 
this communication. Columns (5)-(8) are the result of calculations using Eqs. (1)-(4). The results for the other 


sets of intersections will be discussed in a forthcoming 
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paper. 
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We attempt to present in this paper a unified 
interpretation of the presently known experi- 
mental data on the electromagnetic form factors 
of two fundamental particles: the proton and the 
neutron. As we shall show, this interpretation 
is fully consistent with the idea that the two 
particles are two different aspects of a single 
entity-the nucleon. The third component of the 
isotopic spin of the nucleon is then used to dis- 
tinguish between the two fundamental particles. 
The new experimental material on the neutron 
form factors,’ which now completes a block of 
information on the proton? and neutron, has 
served as the stimulus for the attempted explan- 
ation. 

We would like to explain the main features of 
the experimental behavior of the Dirac form 
factors (F1,,F4,,) and Pauli form factors 
(F 2p F 2») of the proton (p) and neutron (1) as 
functions of the momentum-transfer invariant 
(¢). We propose to do this in a well-known way" 
by expressing each proton and neutron form 


factor as a sum of a scalar and vector contri- 
bution. This decomposition is rooted in the idea 
that the scalar and vector form factors are simp- 
ler and more basic than those of either the pro- 
ton or neutron. Accordingly we make the follow- 
ing definitions, which are standard except per- 
haps for the normalizing constants: 


Fis =F iy Pin? (1) 


Fiy Kip Fin’ (2) 
F,,=([1.79F , pt L-9DF 5, I/(-0.12), (3) 
Fy =[1-79F , » (-1.91)F ,, J/(3.70). (4) 


This choice of normalization has the advantage 
that at g=0 all four isotopic form factors take 
on the value of unity. 

We shall now attempt to find the four isotopic 
form factors from the experimental information 
given in Hofstadter et al.’»? for the values of 


293 








VoLuME 6, NUMBER 6 


PHYSICAL REVIEW LETTERS 


Marcu 15, 1961 





F ip: Fin, F2p, Fan. A comprehensive, though 

approximate, fit to all the experimental data can 
be represented by the following expressions for 

the fundamental isotopic form factors: 


0.56 


F g70-h+ 1+0.214q"’ (5) 


1.20 


F yy = ~0.20+ 0g” (6) 


1.20 


Foy = -0.20+ 700g’ (7) 
-3.0 


Fo5=4-9+ 7 oa14g" (8) 


These results have very few independent fitting 
parameters. The independent parameters of 
Eqs. (5), (6), and (7) are only 0.44, -0.20, and 
a root-mean-square radius, a=0.85x107 cm. 
All other numerical values in these equations 
are determined, once the above choice is made. 
The rms radius is obtained from the coefficient 
(-4a*) of ¢ in the expansion of Fjs5, Fjy, and 
Fey in powers of q. The quantity Fgg in Eq. (8) 
requires the additional fitting parameter 4.0. 
Fas is the least well-known quantity of the set 
of isotopic form factors, and we regard both its 
values and its form as somewhat indeterminate 
at the present time. 

Equations (5) to (8) are remarkably simple and 
have the same fundamental structure, namely, 
the Clementel-Villi (C-V) form.*~™ It is very 
satisfactory that this simple C-V form is also 
suggested by the dispersion relations idea’ that 
the approximate nucleon form factor is just the 


PROTON FORM FACTORS 


FORM FACTOR 
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NEUTRON FORM FACTORS 
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result of a pole plus a constant representing the 
core of a nucleon. The Fourier transform of the 
C-V form factor is a delta function at the center 
of the distribution (r=0) plus a Yukawa cloud. 
Thus the spatial interpretation of Eqs. (5) to (8) 
is very clear: Each form factor corresponds to 
a distribution in space of a simple Yukawa 
cloud and a point-like core. (Our present ex- 
periments are not capable of distinguishing be- 
tween a point core and a core of radius com- 
parable to a nucleon Compton wavelength.) 
Though we are aware that the spatial transform 
is not a completely consistent relativistic con- 
cept, we believe that the density distributions 
so obtained are approximately correct and cor- 
respond, at the same time, to dispersion theory 
models.* 

We may solve Eqs. (1)-(4) for F ip, Fin; F 2p, 
and F9, and substitute the values of Fis, Fiy, 
Fay, and Fgs given in Eqs. (5)-(8). Thus we 
obtain 














aes ony I herd (9) 
F 4, °SA8 +4 aoatig “FZ ard (10) 
<i nay 1 etre (11) 
Fo, ° 38-7 ony I ah og (12) 


Graphical representations of these equations 
are given by the solid lines in Fig. 1. Experi- 
mental points of references 1 and 2 are also 


FIG. 1. The experimental 
points of references 1 and 2 to- 
gether with theoretical curves 
(solid lines) representing Eqs. 
(9)-(12). Below g?=7 no 
points are shown but the theo- 
retical curves are in good 
agreement with previously 
published cross sections. ° 
The dashed lines refer to the 
empirical curves given in 
reference 2. An improved 
fitting of the data at q?>7 by 
the method of least squares 
is now in progress and can be 
carried through with slight 
adjustments of the constants 
in Eqs. (5)-(8). The solid 
and hollow triangles refer 


2 -2 
— to the neutron data of Sobottka.° 
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given in the figures. Even though a best adjust- 
ment of the free parameters has not yet been 
made, the fit between experiment and theory is 
already satisfactory over the entire range of 
values of g and embraces all the measurable 
quantities F1», Fyn, F2p, and Fay. The largest 
departures from the curves correspond to the 
points of Sobottka® which we have analyzed, but 
which were measured at an earlier time when 
our spectrometer was not stabilized by magnetic 
flux-coil methods. 

The coefficients of g* in the expressions (9)- 
(12) give immediately the rms radii of the Dirac 
and Pauli charge and magnetic moment distri- 
butions in the proton and neutron. We find that 


6,, 70.08 f; an = 0.00 f; 


a, =0.94f; a 


2p 170-76 £. (13) 


2 
These radii are consistent with known facts about 
these distributions. We note the important point 
that the root-mean-square radius of the neutron 
is zero, in agreement with the measurements 

of Fermi, Rabi, Hughes, Havens, and their 
collaborators on the neutron-electron inter- 
action. The rms magnetic radius of the neutron 
is nearly the same as that of the proton. One of 
the conditions employed in finding the parameters 
of Eqs. (5)-(8) was that a},,=0. Thus the long- 
standing problem of a small or zero neutron 
charge radius and a normal magnetic radius 
seems to be resolved. 

The splitting of F;, and Fg, at small values of 
¢ is perfectly in accord with known data on the 
proton cross sections, as may easily be verified 
by substituting such form factors into the Rosen- 
bluth formula. 

The choice of positive sign for values of Fy, 
was required for the above set of isotopic form 
factors. If negative values of Fy, are taken for 
the intersections of reference 1, a different set 
of isotopic form factors is obtained which seems 
difficult to understand in any simple way.”»”° It 
is possible, in principle, to find the sign of Fj, 
relative to Fj» by making elastic scattering 
measurements in the deuteron. The present ex- 
perimental evidence» is not definitive on the 
question. If the relativistic correction of Blanken- 
becler** is employed together with the results of 
references 11 or 12, the choice of Fy, should 
be positive, as we have suggested. But the ex- 
perimental errors do not permit a definite deci- 
sion on this point. 

We now find the Fourier spatial transforms of 
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FIG. 2. The proton and neutron charge density dis- 
tributions given by the Fourier transforms of Eqs. (9) 
and (10). The expressions for the spatial transforms 
are given in reference 6. The delta functions at the 
origin are omitted in this figure. 


Eqs. (9) and (10) and present the results graphi- 
cally in Fig. 2. It may be seen from Eqs. (9) and 
(10) that the neutron charge distribution is ob- 
tained from that of the proton essentially by 
flipping over one of the two Yukawa clouds. Thus 
the neutron and proton charge clouds are ina 
partial sense mirror images of each other. The 
fact that the cores are different (0.12 for the 
proton, 0.32 for the neutron) is probably a con- 
sequence of the inexact nature of our approxi- 
mation. It seems quite likely that the higher 
order terms in Eqs. (9) and (10), which are 
omitted in our analysis, might account for the 
actual differences between Fj y and Fay which 
we ignored in our approximate choice of isotopic 
form factors. Such higher order terms may well 
restore full symmetry between neutron and pro- 
ton. 

The magnetic moments of proton and neutron 
are found from the combinations (Fj, +1.79F 25) 
and (Fj, -1.91F9,), and have the approximate 
mirror symmetry expected of them. The de- 
tails of the magnetic moment clouds will be pre- 
sented in a subsequent communication. 

We call attention particularly to the prediction 
that the neutron charge cloud has a positive outer 
fringe.’* The positive sign of Fj, is connected 
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with the positive outer cloud. It would be inter- 
esting to seek other experimental evidence on 
the sign of the outer cloud. 

We also note the fact that the ranges of the 
component Yukawa charge clouds in the proton 
(or in the neutron) are different. The vector 
cloud has a range of approximately 0.32 fermi 
and the scalar cloud a range of approximately 
0.47 fermi. Thus, this evidence indicates that 
the three-pion resonance has a lower energy than 
the two-pion resonance.”*® An improved adjust- 
ment of the constants in Eqs. (5)-(8) might change 
this quantitative relationship between the two 
resonances. 

If the above considerations prove to be true, 
the scheme of construction of proton and neutron 
is simpler than might have been expected. Fur- 
thermore, the internal consistency of the results 
suggests that the techniques of quantum electro- 
dynamics are still valid at distances whose 
values lie between a nucleon Compton wave- 
length and a pion Compton wavelength. 

We would like to offer our special thanks to 
Professor S. Fubini of the University of Padua 
for his important suggestion to choose positive 
Fj» Values and for his illuminating comments on 
the dispersion-theoretic aspects of the C-V form 
factors. We also appreciate his gracious en- 
couragement during the course of analysis of 
the form factors. We wish to thank Professor 
L. I. Schiff also for similar suggestions and for 
his warm and constant encouragement during the 
course of this work. We are grateful to Profes- 
sor G. Breit for his kind critical comments. 
Finally we wish to express our appreciation to 
Mrs. Penny Bakey Seligman and Dr. Denos Gazis 
for their generous assistance with some of the 
calculations. 
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PION PRODUCTION BY NEGATIVE PIONS* 


Barry C. Barish, Richard J. Kurz, Paul G. McManigal, Victor Perez-Mendez,f and Julius Solomon 
Lawrence Radiation Laboratory, University of California, Berkeley, California 
(Received February 13, 1961) 


The reaction 17> +p—1-+n++n, for T, incident 
=365 Mev and 432 Mev, was studied. In a previ- 
ous experiment Perkins et al. investigated this 
reaction. Their total cross section was much 
larger than predicted by static- model theory. 
Rodberg suggested that the large difference be- 
tween experiment and theory was explainable by 
the inclusion of a pion-pion interaction.? Perkins 
also noted a possible deviation of the energy spec- 
trum of the observed pion from a statistical dis- 
tribution. The significant interactions at these 
energies are 


a +p+7 +p, (1a) 
~7°+n, (1b) 
—1~+n* +n, (1c) 
~~ +n7°+p, (1d) 
1° + 7° +n. (le) 


The 7+ from Reaction (1c) with 30 Mev <T_+< 180 
Mev was observed at @=20, 50, 80, and 110 de- 
grees. The energy of the 1+ was determined by 
a magnetic spectrometer. 

The m~ incident beam was produced by the pri- 





of the machine and yield a parallel beam outside 
the shielding. The 7~ beam was momentum-ana- 
lyzed by a 40-deg bend in an H magnet and focused 
at the target by a doublet quadrupole magnet. 
After collimation the cross section of the beam 
was } in. x13 in. The properties of the pion 
beams as determined by range curves and calcu- 
lations are listed in Table I. The liquid hydrogen 
target was a 4-in long 2-in. diameter Mylar cyl- 
inder whose axis lay along the beam direction. 
The target was constructed so that 300 deg could 
be viewed in the laboratory, obstructed by only 
0.031 in. of Mylar. The experimental arrange- 
ment of the target area is shown in Fig. 1. For 
the spectrometer, a C magnet with 13 x24-in. 
pole pieces and a 4-in. gap (maximum field=19 
kgauss) was used in a wedge configuration® to 
obtain focusing in the horizontal plane. The 7;, 
mj’ counters were located at the image of the tar- 
get for a momentum corresponding to the desired 
m+ energy. The energy resolution and solid angle 


Table I. Properties of 1- beams. 
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of the 7+ channels were determined to within +5% 
by wire-orbit measurements. A n* count required 
a coincidence (12347;7,’). Protons with momen- 
tum equal to the 7+ momentum were excluded by 
appropriate thicknesses of Cu absorber inserted 
between 7; and 7;’. For the monitor coincidences 
(12) a fast coincidence (4-nsec resolving time) 
followed by a 40-Mc/sec transistorized discrimi- 
nator-scaler circuit was used. The high 7~ beam 
intensities and cyclotron beam rf fine structure 
introduced problems in direct counter monitoring 
of the > beam. Therefore, a thin-walled ioniza- 
tion chamber was used as an absolute monitor, 
and was calibrated versus the (12) coincidences 
at lower beam intensities. The accuracy of the 
monitor was +2 % at 365 Mev and +0.5% at 432 
Mev. 

The 7*-detection system did not exclude posi- 
trons of the proper momentum resulting from Re- 
action (1b) either by conversion of y rays from 
the usual decay (n°~y+y) or directly through 
Dalitz decay (1°+y+e++e~). The number of 
conversion e+ was experimentally determined 
by inserting Cu converting material between the 
target and Counter 3 in multiples of the amount 
of converting material in the target (= 0.01 radia- 
tion length). The slope of the plot of yield vs con- 
verting material indicates the number of con- 
version e* entering the spectrometer magnet. 
These measurements were made at @ =20 deg 
for 365 Mev and 432 Mev. The measured values 
agreed with calculations of do/dEdQ for con- 
version e+ made using the differential cross- 
section data on 17> +p+7°+n by Caris et al.* 
d’°a/dEdQ for Dalitz e+ was also calculated by 
using the Caris data and the energy and angular 
distribution of the Dalitz e* in the 7° rest sys- 
tem.® The total positron contamination at @ =20 
deg was ~0.17 at 365 Mev and =0.12 at 432 Mev. 
The values decrease rapidly as 6 increases. The 
calculated values were used to correct the data. 
The efficiency of the 1* detection system was 
calculated. Losses of 1* + p++v (=2% to 7%), 
Coulomb scattering, and nuclear absorption of 
the 7+ were taken into account. 

The results for d0/dT*dQ*, do/dQ*, and oiotal 
are presented in Table II. Figure 2 illustrates 
a typical 7+ energy distribution along with a phase- 
space statistical distribution normalized to give 
equal do/d2*. At 365 Mev, do/dQ* is isotropic. 
At 432 Mev, do/dQ* was fitted to a curve of the 
form a,+a, cosé@*. The coefficients a), a, ob- 
tained by the least-squares method are presented 
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FIG. 2. (a) d’o/dT“dQ* versus T*,, for Ts oident 
= 365 Mev and 6,+ =50 deg. The solid curve is a 
phase-space distribution normalized to da/dQ*. (b) 
The ratio of d*c/dT*dQ* (observed) and normalized 
d*a/dT*dQ* (phase space) versus w, the total energy 
for the 7~ and neutron in their barycentric system. 
The errors indicated include that of d*c/dT*d2* and 
the uncertainty in the normalization of d*0/dT *dn* 


phase space. 


in Table III with those of Perkins et al.’ for com- 
parison. 

The total cross sections agree with Perkins 
et al.’ and therefore substantiate the experimental 
deviation from calculations based on static models.* 
The inclusion of a pion-pion interaction by Rod- 
berg” produces agreement with the measured total 
cross sections but predicts a pion energy spec- 
trum favoring higher energies than phase-space 
statistical distributions. The striking feature of 
the measured z* energy distributions is the low- 
energy peaking. The laboratory- system energy 
and angle of the observed n+ determines w, the 
total energy of the 7~, and n(7 =3/2) in their 
barycentric system.” The isobar model® assumes 
that final-state 7-N combinations in the resonant 
state (T=3/2, J=3/2, w=1230 Mev) are preferred. 
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Table Il. Differential cross sections with respect to 7+ energy and angle; 7+ angle; and total cross sections 
(the asterisk denotes total center-of-mass quantity). 











Ts me, il e* d*o/dT*an* do/da* or 
(Mev) 6 (deg) (Mev) (deg) (ub/sr-Mev) (ub/sr) (mb) 
365 +15 20 35 +4.2 33 2.26 +0.19 168 +14 2.07 +0. 09 
48 +6.0 31.5 2.04 +0.17 
65.5 +6.8 31 1.48 +0.13 
83.5 +8.2 30.5 1.00 +0. 08 
50 37.5 +4.3 78.5 2.17 +0.16 174 +13 
52.5 +4.8 75.3 2.05 +0.15 
70.2 +7.0 73 1.53 +0.11 
87.3 +8.6 71.5 1.27 +0. 09 
105 +11.5 70.5 0.77 +0. 06 
80 37.5 +1.2 116 1.87 +0. 23 148 +14 
53.7 +3.2 109 1.99 +0.15 
71.0+5.3 107.5 1.29 +0. 09 
96.2 +6.1 105 0.69 +0. 06 
115.0 +5.1 104 0.45 +0. 05 
110 53.0 42.5 138 2.29 0.31 168 +17 
76.1+4.3 135 1.82 +0.13 
96.0+7.2 133 0.98 +0. 08 
115.0 +4.5 132 0.73 +0. 08 
432 +15 20 19.5 +3.3 37 3.11 +0.33 297 +24 3.26 +0.14 
35.0 +4.3 34.5 3.81 +0.32 
53.0 +6.1 33 3.67 +0. 27 
69.5 +8.5 32 2.46 +0.20 
90.0 +8.5 31 1.51 +0.14 
50 35.0 +4.0 83 3.21 +0. 24 277 +22 
51.5 +5.0 78 2.74 +0.21 
68.7 +8.5 74.5 2.88 +0.20 
86.5 +8.3 73 2.21 +0.16 
115.3 +9.7 72 0.81 +0.08 
80 38.5 +2.3 116 2.20 +0. 27 239 +22 
61.0 +3.6 110 3.16 +0. 25 
79.3 +5.5 108 2.47 +0.18 
101.0 +7.5 106 1.76 +0.13 
117.3 +8.8 105 0.95 +0. 09 
110 58.2 +3.2 132.5 2.17 +0.21 232 +23 
81.0 +4.3 131.5 2.04 +0. 16 
101.5 +5.5 130.6 1.88 +0.14 
126.8 +9.0 130.5 0.85 +0. 07 
146.7 +5.0 130.5 0.56 +0. 09 
At our incident 7” energies the total center-of- resonance is pronounced. Anisovich® has pre- 
mass energy available requires a low-energy m* sented a model which explains the enhanced total 
if w is to approach 1230 Mev. The low-energy cross section for Reaction (1c) and the strong 
peaking agrees with the qualitative predictions influence of the (3,3) resonance on the n* differ- 
of the isobar models. In Fig. 2(b), the ratio of ential distributions, using only the resonant pion- 
d°o/dT*dQ* (observed) to d*o/dT*dQ* (phase nucleon final-state interaction. 
space) is plotted versus w for all experimental We should like to acknowledge the continued 
Points at both energies. The increase of this interest and support of Professor A. C. Helmholz 
ratio in the region of the (3,3) pion-nucleon and Professor Burton J. Moyer. We also wish to 
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Table Ill. Least-squares-fit coefficients for da/dn* 
=@)+a,cos@ at T,- incident = 432 Mev. 








@ (ub/sr) a;(ub/sr) 
This experiment 259 +11 47 +21 
Perkins et al. 268 +20 138 +36 





thank Mr. James Vale and the cyclotron crew for 
their assistance and cooperation during the course 
of the experiment. 





*Work done under the auspices of the U. S. Atomic 
Energy Commission. 

+Present address: Hebrew University, Jerusalem, 
Israel. 
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RESONANCE IN THE K-7z SYSTEM* 


Margaret Alston, Luis W. Alvarez, Philippe Eberhard,t Myron L. Good,? 
William Graziano, Harold K. Ticho,!! and Stanley G. Wojcicki 
Lawrence Radiation Laboratory and Department of Physics, University of California, Berkeley, California 
(Received February 16, 1961) 


In a continuation of the study of the interaction 
of 1.15-Bev/c K' mesons in hydrogen by means 
of the Lawrence Radiation Laboratory 15-inch 
hydrogen bubble chamber, we now report a study 
of the reaction! 


K~ +p~K° +27 +p. (A) 


Examples of this reaction were easily identified 
in those cases in which the K° decayed into 
charged pions and appeared in the chamber as 

a two-prong interaction associated witha V. A 
kinematic analysis isolated 48 events of reaction 
(A) from other events with similar topology.? In 
only one case was the identification not unique. 
Correcting for neutral decays of the K° and for 
escape from the chamber, we find a total cross 
section of 2.0+0.3 mb for Reaction (A). 

The events are shown on a Dalitz plot in Fig. 1. 
If the reaction were entirely dominated by phase 
space, the Dalitz plot would be uniformly popu- 
lated. Instead, a strong clumping around proton 
kinetic energy of 20 Mev is observed. This effect 
cannot be explained by an interaction matrix ele- 
ment that increases monotonically with decreas- 
ing proton energy. Whereas an extrapolation 
from the region 15 Mev < Tp < 25 Mev would lead 
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one to expect a minimum of 16 events in the re- 
gion T, < 15 Mev, only three are found there. No 
experimental bias against very low energy pro- 
tons in the K-p center-of-mass system can exist, 
since such protons have laboratory-system mo- 
menta of approximately 600 Mev/c, and are 
easily identified. The observed distribution can 
best be explained by a quasi-two-body reaction 
of the type 


K"+p-K* +p, (B) 
followed by a decay, 
K*"~K°+7°. (C) 


The 3-3 resonance of the pion-nucleon system 
would show itself on the Dalitz plot as a concen- 
tration of points along the diagonal line drawn 
through Fig. 1. The absence of any evidence for 
this resonance in our data can be explained if 
Reaction (A) proceeds primarily through the /=0 
channel, which cannot produce a p-7~ system in 
the J=3/2 state. Further, even in the J=1 chan- 
nel, the 3-3 resonance favors (n+ 7°) +K° over 
(p+77)+K°, and hence provides additional sup- 
pression of this resonance. 

The mass distribution of the K*~ is shown in 
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Fig. 2. The mean value is 885+3 Mev. After 
removing the number of background events esti- 
mated from the phase-space distribution and un- 
folding the experimental error on each of the re- 
maining 22 events (typically 3 to 4 Mev), we ob- 
tained a full width at half maximum of 16 Mev, 
corresponding to a lifetime of 4 x 10-?° second. 
The angular distribution for Reaction (B) is 

consistent with isotropy. Assuming that the K*” 
system is produced predominantly in the s state 
(which appears likely both because of the close- 
ness to the threshold and because of the isotropic 
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FIG. 2. Mass spectrum of the K®-1- system. The 
solid line represents the phase-space curve normalized 
to background events. 
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distribution), we can obtain an upper limit of the 
K* spin S. If S=0, then the reaction can be pro- 
duced only through the p,, ingoing channel; if 
S=1, through s,, andd,,, etc. In any case the 
decay angular distribution® is given by 


1(8) = la¥ o(% g) |? + DA 1 ¢) |?, (1) 


with |a|?+|b|?=1, where a=1 for S=0. Here 6 

is the angle of the K° in the K*~ rest system with 
respect to the incoming K™ direction. The mean 

value of cos*@, based on the distribution function 
(1), is 


(2S? + 2S - 3)+2\a|? 
4S* + 4S - 3 


Experimentally, (cos*@) for the 21 events lying 
in the K*” mass range between 870 and 900 Mev 
is 0.275. Using S=2 in Eq. (2), we find the ex- 
pected value of (cos*@) > 0.429, with a standard 
deviation of 0.051. The experimental result thus 
deviates from the range of values expected for 
S =2 by three standard deviations. For S> 2 the 
discrepancy is even greater. On the other hand, 
the experimental result is consistent (within 
errors) with S=0 or S=1. It is worth noting that 
an isotropic decay distribution is obtained both 
for S=0 and for S =1 if the d,, input channel does 
not contribute. For this reason, if S=0, experi- 
ments at several momenta will be required to 
settle this problem. 

Assuming that the K~ and K° are an isotopic- 


(2) 





(cos?@) = 
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spin doublet, and /=1/2 for the K*~ isotopic 
spin, the branching ratio R = (K*” ~K™ +7°)/ 
(K*- +K°+7~) equals 1/2; for J=3/2, R=2. In 
either case another charge state of K*, namely 
K*°, should exist and decay into K~ + 7+ or K°+7°. 
The value of this branching ratio, also, depends 
on the isotopic spin of the K*. 

To investigate the isotopic spin properties of 
the K*, we searched for examples of the following 
two reactions: 


K~ +p+K~+7°+p, (D) 
K~+p+K~ +nt+n. (E) 


These events appear as two-prong interactions 
and are much more difficult to identify than the 
events already discussed, since there are usually 
several possible interpretations for each inelastic 
two-prong event. In particular, there was a pion 
contamination of about 10% in our incident K™ 
beam, and the inelastic pion interactions are 
kinematically very similar to Reactions (D) and 
(E). Both the kinematic fits and the ionization of 
the tracks were used to identify the events, but 
these criteria were not always sufficient to dis- 
tinguish between various hypotheses. At present 
we have processed only about 2/3 of our two- 
prong interactions, but we feel that the data ob- 
tained are reasonably unbiased. 

In both Reactions (D) and (E) there are peaks 
in the nucleon kinetic energy distribution in the 
R* resonance region. On the basis of the number 
of events in the proton peak of Reaction (D), our 
present data allow us to make a crude estimate 
of the branching ratio: R=0.75+0.35. The data 
thus strongly favor the J/=1/2 state. 

The experimental production ratio of K* via 
Reactions (A) and (E) is about 1, and is thus con- 
sistent with the production of the K* through a 
pure isotopic spin state. 

An I=1/2 particle, called the K’, with negative 
parity with respect to the K meson, has been in- 
voked by Tiomno to explain the backward A peak- 
ing in associated production.* Gell-Mann postu- 





lated the existence of such a particle to permit 
the construction of a strangeness-violating weak- 
interaction axial current.’ The K* that we have 
observed has properties consistent with those 
postulated for the K’, but, as discussed above, 
the K* spin and parity remain to be established. 
As in our previous communication, we acknowl- 
edge gratefully the assistance of the many people 
who helped us to obtain and analyze these data. 
One of us (P.E.) is grateful to the Philippe’s 
Foundation, Incorporated, and to the Commi- 
sariat a |’ Energie Atomique for a fellowship. 
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CHARGED HYPERON PRODUCTION BY 16-Gev/c 7” MESONS 


J. Bartke,* R. Bick, R. Budde, W. A. Cooper, H. Filthuth, Y. Goldschmidt-Clermont, 
F. Grard,f G. R. MacLeod, A. Minguzzi-Ranzi, L. Montanet, W. G. Moorhead, D. R. O. Morrison, 


S. Nilsson, C. Peyrou, B. W. Powell, J. Trembley, and D. Wiskott 
CERN, Geneva, Switzerland 


I. Bertanza, C. Franzinetti, I. Manelli, and V. Silvestrini 
Pisa University, Pisa, Italy 


and 


G. Brautti, M. Ceschia, and L. Chervosani 
Trieste University, Trieste, Italy 
(Received February 27, 1961) 


The CERN 30-cm hydrogen bubble chamber 
was exposed to a beam of 16-Gev/c 7~ mesons. 
In a sample of 35000 pictures, 48 charged strange 
particles were found and the results obtained from 
them are given below. The events were meas- 
ured on digitized machines. Measurements on 
three views were assembled, and helices fitted 
to the tracks by a least-squares method, using 
a Ferranti Mercury computer.’»? A further pro- 
gram tested various two-body decay kinematics.* 
By this means, 7-p decays and p scattering 
events were easily recognizable, but only ina 
few cases was it possible to distinguish kinemat- 
ically between the decays £* +a*+n (“Ly”), 
Lt+p+n° (“Zp”), and K +7 or pu (no evidence 
for = or anti-= was found). It was possible to 
distinguish between pion and proton secondaries 
in all cases by means of ionization measure- 
ments using a mean gap length method. Ona 
normal track the identification is better than one 
standard deviation up to 1.2 Gev/c and better 
than two standard deviations up to 0.8 Gev/c. 

The composition of the sample is given in Table I. 

In many cases the primary track was short 
(<5 cm) and its momentum therefore poorly 
measured. The assumption was therefore made 
that all charged V’s were hyperons, unless they 
could be shown to be K mesons (two double 
events were excluded from the present analysis 
since the above assumption was not justified in 
these cases). The primary momentum could 


Table I. Composition of the sample. 





z= positive 28 
z negative 18 
K positive 2 
™p 25 
bp scattering 2 





then be calculated from the secondary momen- 
tum and angle. 

In 11 cases the calculated primary momentum 
had two energetically possible values. In the 
absence of a means of distinction, the further 
assumption was made that the two values were 
equally probable. 

The assumptions were justified by a lifetime 
determination which gave 75 - = (2.17}:?) x10 
sec and T+ = (0.607$-77)x107° sec, in agree- 
ment with accepted values. The angular distri- 
bution of the decays was compatible with isotropy 
in the c.m. system. 

In order to find the number of hyperons pro- 
duced, the probability of observation of each 
decay had to be determined. For this purpose 
3 cm were subtracted from the potential path 
length of the hyperon in the chamber (since de- 
cays with very short tracks were unrecognizable) 
and the probability of decay in this distance found. 
Each event was weighted by the inverse of this 
probability. The mean weight was 1.8 for =~ 
and 1.3 for =*. All results which follow are 
weighted in this way. 

Scanning bias was looked for by plotting the 
distribution of y, the angle between the normal 
to the decay plane and the normal to the front 
glass. Because the angle between the incident 
track and the hyperon was small, one expects 
cylindrical symmetry so that the distribution of 
y, rather than cosy, should be independent of y. 
No significant bias due to foreshortened decay 
angles was found on a plot for all events. Num- 
bers were too small to plot + and =~ separately. 

Scanning bias can also be checked by the num- 
bers of Z»* and =,* observed which should be 
equal. In fact it was found that Zp*/Z_*=1/4. 
There was therefore bias operating against the 
Zp mode. The decay angle in this mode is about 
6 times less, on average, than in the 2q mode. 
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Hence one would expect the Z, decays to be 
missed either when the decay angle is much fore- 
shortened (cosy<1/6, i.e., 80°<y~<100° which is 
too small a region to be visible on the y plot, with 
our limited statistics) or when the momentum 
is about 6 times higher than the average. 

These small biases will act equally on £,* 
and 2,” decays. It is therefore safe to take the 
ratio of £,*/Z,~ as unbiased. The charge ratio 
can then be found by assuming that as many £,* 
as D,* had occurred. Taking the 22 £,* and 
18 Z,° and applying the corrections, we find 


No. at production of =* 
No. at production of =~ 





+ + - + - +0.7 
= =2 =i. “3° 
(Z 7, Ve z. [z 1.8_5'5 


By assuming that the numbers of =* and =~ 
have Poisson distributions with the means 22 
and 18 (the numbers actually found), one can 
compute the expected distribution of the ratio 
=*+/Z-. The errors given above refer to the 
width at half height of that distribution. 

The numbers at production correspond to cross 
sections of 0.35+0.08 mb for £* and 0.19+ 0.05 
mb for =~. These values may be compared with 
the Dubna results at 8 Gev.*»® They find 0.1 
mb for the =~ production by 7~ mesons on pro- 
tons. The charge ratio is not given explicitly 
but seems to be less than 1. 

The momentum distribution of the hyperons at 


No. No. 
10 5 - 








production is best obtained from those which de- 
cay via the 2, mode since they are unbiased up 
to high momenta (i.e., higher than those of the 
observed hyperons). In the center-of-mass sys- 
tem of the colliding particles, the maximum 
momentum of either of a DK pair is 2.6 Gev/c. 
As shown in Fig. 1(a) and (b), the events analyzed 
had values ranging up to this maximum. The 
mean values were not significantly different for 
=~ and =*: 

(p*) for 2~=1.15+4 0.13 Gev/c, 

(p*) for £* =1.25+ 0.11 Gev/c. 
The mean transverse momenta were: 

(p,) for =~ =0.54 Gev/c, 


(p,) for =* =0.58 Gev/c. 


If the lower momentum value is taken in the 
eleven ambiguous cases referred to above, lower 
limits are found: 


(p,) for Z~ > 0.38 Gev/c, 
(p,) for =*>0.56 Gev/c. 


These are larger than the mean transverse mo- 
mentum of 7 mesons produced in the interactions, 
which is 0.33 Gev/c.® 

The angular distribution of the hyperon in the 
m-p center-of-mass system is shown in Fig. 2(a) 
for Z,~ and in 2(b) for =,* and is seen to be 
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FIG. 1. Momentum distri- 
bution at production in the 1™ +p 
ee c.m. system (a) for £,~, and 








(b) for £,*. —— corrected 
numbers; --- uncorrected num- 
bers. 
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FIG. 2. Angular distribu- 
tion at production in the 1" +p , 
c.m. system (a) for £,°, and ; a 
(b) for Z,*. —— corrected 4%. . 
numbers; --- uncorrected num- ; b L 
bers. u 
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peaked backward. This and the preponderance 
of positive hyperons are the most striking fea- 
tures of the results. Both are consistent with 
the interpretation of production via single meson 
exchange. In this case it is possible to make 
more detailed predictions. A comparison of 
these predictions with the present results, and 
those from further measurements on neutral and 
charged strange particles will be published in 
due course. 

No significant difference was found between 
the multiplicity of jets with a charged V decay 
(4.4+0.2) and of jets where no strange particle 
decay was visible (4.3 + 0.1). 

One of us (J.B.) wishes to thank the Ford 
Foundation for a fellowship during this work. 
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EXISTENCE OF THE w°® PARTICLE* 
R. G. Sachs and B. Sakita 


Department of Physics, University of Wisconsin, Madison, Wisconsin 
(Received February 20, 1961) 


Recent experiments by Abashian, Booth, and 
Crowe’ on the p-d reaction have been interpreted 
as an indication of the existence of an w® particle 
(boson) having an isotopic spin of J/=0 or J=1 
(probably the former), a mass in the neighbor- 
hood of 305 Mev, anda rather long lifetime for 
decay into pions (width less than 16 Mev). It has 
already been remarked’ that such a state of the 
two-pion system appears to be inconsistent with 
the pion energy distribution in 7 decay. We would 
like to point out that the existence of either a two- 
pion or a (bound) three-pion state* of angular 
momentum J=0 or J=1 at this energy would ap- 
pear to be inconsistent with observations on the 
K*-decay modes. The point is that the decay 
process 





K* = w° +n" (1) 


could occur and would have the appearance of an 
anomalous (27) decay mode of the K* when the w® 
decayed into neutral particles (such as a 7° and 
y). Estimates of the branching ratio for this 
process which are given below indicate that it 
is of the same order or larger than the branch- 
ing ratio for Kg3 decay; hence it should have 
been readily detected and has not been. 

In order to estimate the branching ratio for 
process (1), we first assume that the w° has 
spin 0 so that the decay into an S state is al- 
lowed. Since | A/| =4 is possible for (1), the 
process is very similar to the 27 decay mode 
of the neutral K meson. If it is assumed that 
the coupling is about the same for the two modes, 
the ratio may be estimated on the basis of the 
available phase space: 


W(K*t —w°+*) 


WK, =n*+0") = (q/p)= 0.45, (2) 





where W is the transition probability and g, p 
are the magnitudes of the momenta of the w°® and 
m~, respectively. Equation (2) indicates that the 
K* lifetime would be nearly as short as the K,° 
lifetime, which is quite clearly contrary to the 
observations. Since the lifetime ratio is about 
100 and the branching ratio for the anomalous 
mode of the K* is presumably much less than 
10%, the above estimate would need to be in 
error by a factor of at least 10° to permit the 
existence of the spin-0 state. Although the esti- 
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mate is crude, it does not seem likely that the 
errors are so large as that. 

If the w° has a spin of 1, the final state in re- 
action (1) is a P state and the centrifugal barrier 
may reduce the ratio from the value given by 
Eq. (2). In order to estimate the effect of the 
barrier, it is necessary to introduce a range 
v,=A7 for the interaction in the decay. We then 
find 


W(K* = w° +7") 
W(K,° -1* +77) 





= (q/p)(e uX . pry, 


where e,, is the polarization vector of the w° and 
K, the four-momentum of the K*. After sum- 
ming over directions of polarization, we find 
W(K* = w° +27) 
W(K,° -1* +77) 


= (m,/m w) (9/)(q/m as =0. 5(m_/r)?, (3) 





where mx, m,,, and m, are the masses of the 
indicated particles. 

Ogawa‘ has demonstrated that an interaction 
range corresponding to A= $m, is suggested by 
the branching ratios for nonleptonic decay of the 
strange particles. However, even if we are much 
more conservative and take A=M, the nucleon 
mass, the ratio Eq. (3) becomes 10~ which is 
still much too large to have escaped detection. 

Because in the case that w° has spin 1 the above 
comparison involves the relationship between a 
scalar (for the 7) and vector (for the w) coup- 
ling, it is desirable to make an estimate on an 
entirely different basis. Since the weak electron- 
neutrino current is a vector current, we have 
estimated the ratio of the rate of (1) to the Ke3 
decay rate on the basis of the diagrams shown 
in Fig. 1. The bubble in the w® line of Fig. 1(b) 
is taken to be a nucleon-antinucleon bubble and 
the weak coupling of this bubble to the black box 
is assumed to have the same form as the leptonic 
coupling. In this way we find 


W(K*t —1*+w) 
W(Kt —1°+e*+v) 


x (q/m (CG /4n) "S/F "Y, (4) 





=7.5x10°(5m /m,)° 


where 6m, is the mass shift of the w® due to the 
nucleon bubble, G is the strong coupling between 
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FIG. 1. Diagrams of the related processes for K* 
decay. The loop in Fig. 1(b) is a nucleon-antinucleon 
loop. 


the w° and nucleons, and f and f’ are the weak- 
coupling constants for nucleons and leptons, re- 
spectively. An estimate of G has been given by 
Sakurai® on the assumption that a neutral vector 
boson is responsible for the spin-orbit force in 
p-p scattering. The value obtained in this way 
depends on the mass of the boson. By extra- 
polating his results as a function of the mass to 
the value m,,, we estimate G?/41=2. To obtain 
imw, we make a Goldberger-Treiman® type of 
analysis of the matrix element of the weak inter- 
action current of the w°, including just nucleon- 
antinucleon states. The result is 5m,,?=m.,,’. 
Then Eq. (4) becomes 


W(K* ~1* +w") 
W(K* =1°+eT+v) 





s(f/f'y. (5) 


We may take f =f’ although there is some reason 
to believe’ that (f /f’)? is larger than one. Thus, 
we find that the anomalous decay mode should be 
competitive with the Kg3 mode, which again is 
much too large a rate. Since Fig. 1(b) describes 
Only one of the possible decay mechanisms, it 





seems likely that the estimate given by Eq. (5) 
is quite conservative. 

While we recognize that arguments of the kind 
presented here are subject to some doubt, they 
certainly suggest® that there does not exist an 
w° particle composed of pions and having the 
mass observed in the experiment of Abashian 
et al." 

However, these arguments do leave open the 
following possibilities: 

(1) An w® state of the pions exists but has a 
spin greater than 1. 

(2) An w® particle with spin 0 or 1 exists but 
its interactions with pions are not very strong. 

(3) The observed effect is merely a manifesta- 
tion of a normal energy dependence of the matrix 
element, as suggested by Tubis and Uretsky® 
and by Greider.’° 

A careful examination of the energy spectrum 
of the charged pion in (apparent) rt’ decay would 
be of interest in this connection since the proc- 
ess (1) would cause a narrow line to be super- 
imposed on this spectrum. 





*This work was supported in part by the U. S. Atomic 
Energy Commission and in part by the University of 
Wisconsin Research Committee with funds provided by 
the Wisconsin Alumni Research Foundation. 
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ANALYSIS OF THE ANOMALY IN DOUBLE MESON PRODUCTION IN p+d COLLISIONS 
AND THE S-WAVE PION-PION INTERACTION* 


Tran Nguyen Truong 
Laboratory of Nuclear Studies, Cornell University, Ithaca, New York 
(Received February 20, 1961) 


The purpose of this note is to show how the 
anomaly in double pion production in p+d colli- 
sions observed by Abashian, Booth, and Crowe’ 
may be explained by a nonresonant final-state 
interaction of the produced pions in the isospin 
state T=0, instead of assuming the existence of 
a neutral w® particle. 

The separation of the production process into 
two separate mechanisms, that of the primary 
interaction in which the particles are produced, 
and the final-state interaction in which the pro- 
duced particles interact with each other, is well 
known.” It is assumed in this note that the final- 
state interaction is responsible for any deviation 
from simple phase-space predictions. In the 
following calculation we neglect the effects of the 
deuteron and He® wave functions and the pion- 
nucleon final-state interaction, since they tend 
to smear out the spectrum and therefore are 
unlikely to produce a sharp peak. We limit our- 
selves to an S-wave (T=0) 1-7 interaction®; the 
next angular momentum D state, which is per- 
mitted by 7-m isospin T=0, is very unlikely to 
contribute.* A straightforward analysis based 
on charge independence yields 


o(p +d +He® +2* +27) +0(p+d +He® +7°+7°) 
=$1f I? + Ifol’, (1) 
o(p+d ~H° +2*+n°)=If,1?, (2) 


where f 7 is the production amplitude in a state 
in which the two pion charge states couple to- 
gether to make an isospin T. Experimentally 
(2) is much smaller than (1), which is what one 
expects for the production of the two pions in 
the nonresonant P state at this energy.‘ In the 
following analysis we neglect the contribution 
of the f, production amplitude to (1). The cross 
section is proportional to 


d*p, d*p, d*p, a ae 
QE, 2E, 26, Py +P2 +P, -Pp) 
_ 2 
XO(E, +E, +E, -E,)IT,,(q/2)!?. (8) 


The subscripts 1, 2, and 3 refer, respectively, 
to the two pions and the He*; qg/2 is the momen- 
tum of each pion in their own center-of-mass 
system. We have put the primary interaction 
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matrix to be constant, and 7;(q/2) is the part 
which is due to the final-state interaction. After 
carrying out the integration in (3), we have 


d’o/dp,dQ, =c(p,*/w4)(q/w ITI", (A) 


where w=’ +4m,” and the He* dynamical quan- 
tities are evaluated in the laboratory system; 
(p;” /w,)(q/wg) is the phase-space volume ele- 
ment. To obtain |74!* empirically, we divide 
the measured spectrum by the normalized phase 
space. Using the data given in reference 1, the 
plot of | 74!" is shown in Fig. 1. For a very 
small radius of the primary interaction, it is 
the square of the ratio of the pion-pion wave 
function with and without 7-7 interaction taken 

at the origin, and is nothing but the usual en- 
hancement factor.° We use an exponential poten- 
tial well of range d= 3f/uc to calculate this en- 
hancement factor; we also use an asymptotic 
m-m wave function with exponential cutoff for the 
range of the primary interaction to obtain the 
energy dependence of the matrix element, and 





—— de *0.54 (S*0.28) 
---- de *0.72 (S*0.35) 
—— de = 1.00 (S*0.42) 
—“— de = 1.27 (S*0.49) 











4 1.36 1.38 140 1.42 


Ps Bev/c 


FIG. 1. The empirical probability density function 
| T¢j|? as a function of the He? momentum. Experi- 
mental data. are taken from reference 1 at incident 
proton energy of 743 Mev. Fits to the data by the 
theory of pion-pion final-state interactions are shown. 
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find that to fit the data one needs a larger scat- 
tering length. A reasonable fit to the data is 
obtained with the well depth parameter s=0.3 to 
0.5 corresponding to a scattering length a,= $f/uc 
to #i/uc (attractive), a, being the scattering 
length for T=0. Since we take into account only 
the effect of m-7 interaction, the result should 
only be regarded as evidence for the T=0 1-7 
interaction being attractive and for its scatter- 
ing length being not too large. It is interesting 
to note that this conclusion is in agreement with 
that obtained by Ishida et al.* and Efremov et al.” 
in their work on the 6,, and 6,, phase shifts in 
pion-nucleon scattering. 

Mitra found that it is possible to fit the t-decay 
data spectrum with either or both T=0 and T=2 
resonances, and with the position of the T=0 
resonance at w,,”= 12m,” which is quite far 
away from the energy region we are consider- 
ing.® A resonance at w,,”=5m,,” would not fit 
the r-decay data. A low-energy T=0 m-n reso- 
nance would also imply a too large inelastic 
cross section for the process 1+N ~27+N near 
the pion-production threshold.’ The analysis of 
the t decay by Lomon, Morris, Irwin, and 
Truong® shows that the m-m interaction is attrac- 
tive for the state T =2 and probably repulsive 
for the state T=0; with both T=2 and T=0 at- 
tractive, and 7=2 more attractive than T=0, a 
fit to the data was also obtained although it was 
not as good. 

It is worthwhile noticing that the previous re- 
mark is not in disagreement with the results ob- 
tained by the dispersion relation methods of Khuri 
and Treiman.'° Combining the present analysis of 
the anomaly in the p+d experiment and the T- 
decay data, it may be possible to infer that the 
S-wave pion-pion interaction in the states T =2 
and T=0 are both attractive with T=2 more at- 
tractive. 

We would like to mention that our analysis does 
not exclude the possibility of the existence of the 
w° particle. Accurate experimental data at dif- 
ferent incident proton energies and He’® angles, 
could be used to tell whether the observed peak 
is due to the final-state interaction of the two 
pions or to the existence of the w® particle. 

We would like to thank Professor H. A. Bethe 
for his advice, encouragement, and interest in 


this work. A number of helpful discussions with 
Dr. R. Schult, Professor T. Kinoshita, and Pro- 
fessor E. Lomon are very much appreciated. 
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and T=0 repulsive with T=0 more repulsive, a) <a,<0; 
or T =2 weakly attractive and T=0 repulsive, a,>0, 

a) <0; or T=2 more attractive than T=0, 0<a)<ap. 

If it was required that a,>0, then only the last solu- 
tion holds. These solutions could be understood quali- 
tatively from the enhancement factor argument. A 
similar calculation to that in reference 5 is being 
carried out under less severe restriction on the 
curvature. 
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TWO-PARAMETER APPROXIMATION TO S-WAVE SCATTERING* 


Daniel M. Greenberger and B. Margolis 
Physics Department, Ohio State University, Columbus, Ohio 
(Received February 13, 1961) 


The purpose of this note is to introduce a 
simple two-parameter approximation to scatter- 
ing amplitudes calculated from the Mandelstam’ 
representation, which is more accurate than 
the usual simple pole approximation. The tech- 
nique is illustrated for the singlet s-wave nucleon- 
nucleon amplitude. 

In the Mandelstam representation the analytic 
structure of the partial-wave scattering ampli- 
tudes are completely determined by the unitar- 
ity requirement. In the case of nucleon-nucleon 
scattering, the s-wave amplitude in the v=¢ 
plane has two branch lines? along the real axis, 
one for v 20 and the other for v < -} [see Fig. 
1(a)]. The usual procedure in solving for the 
amplitude h(v) is to write’ h(v) =N(v)/D(v), where 
N(v) is analytic except for a branch cut for v < -4 
and D(v) is analytic except for a branch cut for 
v>0Q; Cauchy’s theorem then states that [per- 
forming one subtraction in D(v) at v9] 


a 4 P 
N(v) == ‘ole, (1) 


-2 
Vv! N(v’) 
m (v'-v,)(v'-v)’ 

(2) 
where f(v)=Imh(v), v <-3, and where we have 
assumed elastic scattering. The simplest approx- 





D(v)=1-7="2 { “av! 
0 


(a) v-plane 





(b) w- plane 








+ 


FIG. 1. Location of singularities of scattering 
amplitude for singlet s-wave nucleon-nucleon scatter- 


ing. 
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imation for solving these coupled integral equa- 
tions is to replace N(v) by a pole, whose strength 
and location are left as free parameters, which 
leads to the effective-range approximation.” 

Our idea is merely to introduce a conformal 
mapping which “localizes” the extended singular- 
ity of N(v) and then to approximate it by a pole. 
This will not be equivalent to merely finding the 
best placement of the original pole. Consider 
the successive mappings v’=(v + 3)”, -1<argy <7 
and w=1/(v’-a), a>}. The first “opens out” 
the left-hand branch cut, while the second trans- 
forms it into a circle [Fig. 1(b)]. Thus, N(w) 
is now singular only along the circle (radius = 1/2a) 
and the low-energy physical region begins at the 
point D [w=-2/(2a-1)]. The point y= has been 
mapped into the origin. 

The procedure next is to replace N(w) by a pole 
at the center of the circle and to calculate its 
effects in the physical region. Thus, 


2 
1 
niw)=r(—1, -<) . (3) 


This choice makes N(w) vanish as w’ at the 
crigin [i.e., N(z)+1/z as z+]. To see that 
this is quite different from a similar approxi- 
mation in the z plane, consider the case of a 
line charge near an L-shaped conducting plate 
with semi-infinite legs. Here, after the mapping 
w =z" into the half-plane, the problem requires 
one image charge, while in the original con- 
figuration it requires three. 

We now have 


N(v) =T'/[(v + 9) +a)’, (4) 


with T and a as our two parameters. Equation 
(2) can now be integrated and the phase shift is 
given by 


(v? /m)cot6 = Re D(v)/N(v). (5) 


The calculation is simplified by the choice v, =0. 
Rather than compare directly with experiment, 
the parameters were fitted to the Jackson-Blatt* 
straight line in order to be able to compare with 
the extrapolation of Cini, Fubini, and Stanghel- 
lini.’ (See Fig. 2.) This two-parameter fit is a 
fair approximation to their curve even at nega- 
tive energies.* For v<-3, w is right on the 
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FIG. 2. Results of the present calculation (GM), 
compared to those of Cini et al. (CFS). The straight 
line (JB) is that of Jackson and Blatt.‘ 


circle so the approximation is not expected to 
be valid. However, by removing the restriction 
that the pole be at the center of the circle, one 
could exactly fit their point at vy = -4 and still 
have a two-parameter fit. Note that by intro- 
ducing higher multipoles one could develop a 
systematic approximation scheme. 

The shape parameter has the same sign as that 


of Cini et al.,° and also as that of Noyes and 
Wong,’ which may well be a consequence of the 
values of the effective range and scattering length 
and the location of the singularities in the Man- 
delstam representation, even though an experi- 
mentally meaningful shape dependence might be 
critically affected by relativistic effects, Coulomb 
effects, and many other diverse effects.* One 
further advantage of this solution is that it pre- 
serves the left-hand branch line, which the pole 
approximation does not, and h(v) is real only 
between -}<v<0. 

Frazer® has developed a conformal mapping 
method for use in the Mandelstam representa- 
tion, though we are unfamiliar with the details 
of this work. 





*This work was supported in part by the U. S. Atomic 
Energy Commission. 

'S, Mandelstam, Phys. Rev. 112, 1344 (1958). 

2G. F. Chew, Lawrence Radiation Laboratory Re- 
port UCRL-9289, 1960 (unpublished). 

5G. F. Chew and S. Mandelstam, Phys. Rev. 119, 
467 (1960). 

‘J. D. Jackson and J. M. Blatt, Revs. Modern Phys. 
22, 77 (1950). 

M. Cini, S. Fubini, and A. Stanghellini, Phys. 
Rev. 114, 1633 (1959). 

®The pole approximation lies far above the straight 
line at negative energies. 

"H. P. Noyes and D. Y. Wong, Phys. Rev. Letters 
3, 191 (1959). 

*We would like to thank Dr. H. P. Noyes for calling 
these effects to our attention. 

*w. R. Frazer, Bull. Am. Phys. Soc. 6, 81 (1961). 





SINGULARITIES OF THE COSMOLOGICAL SOLUTIONS OF GRAVITATIONAL EQUATIONS 


I. M. Khalatnikov, E. M. Lifshiftz, and V. V. Sudakov 
The Institute of Physical Problems, Academy of Sciences, Moscow, U.S.S.R. 
(Received February 28, 1961) 


In cosmological applications of the general 
relativity theory extensive use is made of the 
well-known (Friedmann’s) solution of the Ein- 
stein gravitational equation which is based on 
the assumption of complete homogeneity and 
isotropy of the space distribution of matter. 

This assumption is far reaching in its mathe- 
matical aspects, not to mention that its fulfill- 
ment in the actual universe could at best be only 


of approximate nature. Hence the question arises: 


To what extent does the important property of the 
resulting solution—the existence of the time sing- 


ularity, depend on this specific assumption? 
The solution of this problem, which is of pri- 
mary importance for the entire cosmology, re- 
quires an investigation of the situation arising 
for a quite arbitrary distribution of matter and 
gravitational field in space. A short summary 
is given here of the results of such an investiga- 
tion. 

The natural choice of the reference system in 
dealing with this problem turns out to be a sys- 
tem, subject to the conditions -g,,=1, 2oq =9, 
a@=1,2,3 (we shall call such a system synchro- 
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nous, since it allows of the synchronization of 
clocks along the entire space). It was long ago 
pointed out by Landau that due to one of the 
gravitational equations (the 00 equation) the 
metric determinant g must inevitably become 
zero in a finite time. However, this result 
(which was recently found independently also by 
other authors’) does by no means prove—con- 
trary to the opinion expressed in the literature— 
the inevitability of the existence of a real (physi- 
cal) singularity in the metric, which cannot be 
excluded by any transformation of the reference 
system. The singularity can turn out to be fic- 
titious, nonphysical, being connected merely 
with the specific nature of the chosen reference 
system. 

An answer to this question emerges from the 
geometrical analysis of the space-time proper- 
ties in the synchronous system of reference. 

It is easily seen that in a synchronous refer- 
ence system the lines of time are geodesics in 
the 4-space. This property can be used for a 
geometrical construction of such a system in 
any space-time. We choose an arbitrary space- 
like hypersurface and construct a set of geodesics 
normal to this hypersurface. If one defines now 
the time coordinate as the length of a geodesic 
between a given world point and the intersection 
with the hypersurface, one arrives, as it is easy 
to see, at a synchronous reference system. 

But geodesic lines of an arbitrary set in gen- 
eral intersect each other on some envelope hyper - 
surfaces —the four-dimensional analogs of the 
caustic surfaces of geometrical optics. Thus 
there exists a geometrical reason for the appear- 
ance of a singularity, which is due to specific 
properties of the synchronous reference system 
and is therefore obviously of a nonphysical nature. 
It is to be emphasized, however, that an arbitrary 
metric of a 4-space in general allows also for the 
existence of nonintersecting sets of time-like 
geodesics. But the above-mentioned property 
of the gravitational equations means that the 
metric admitted by them excludes the possibil- 
ity of the existence of such sets, so that the lines 
of time necessarily intersect each other in any 
synchronous reference system. 

This means, from the analytical point of view, 
that in a synchronous system of reference the 
Einstein equations have a general solution with 
a fictitious singularity with respect to time. 

Thus any foundation is removed for the exis- 
tence, along with this general solution, of yet 
another, which would also be a general one but 
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would have a real singularity. The criterion of 
the generality of the solution is the number of 
arbitrary functions (of the space coordinates) it 
contains. Among these functions there are in 
general also such, whose arbitrariness is due 
merely to the freedom in the choice of the refer- 
ence system admitted by the equations. What is 
essential is only the number of the “physically 
different” arbitrary functions, which cannot be 
decreased by any specific choice of the refer- 
ence system. For the general solution this 
number must be eight; these functions must 
provide for the possibility to put arbitrary initial 
conditions, determining the initial space distri- 
butions of the density and the three velocity 
components of the matter, and of the four quan- 
tities which determine the free gravitational 
field. (The latter number can be arrived at, 
e.g., by considering weak gravitational waves; 
since these waves are transverse, their field 

is characterized by two quantities which obey 
differential equations of the second order, and 
therefore the initial conditions for this field 
must be given by four space functions.) 

The above geometrical considerations do not 
exclude, of course, the possibility of the exis- 
tence of narrower classes of cosmological solu- 
tions with a real singularity. Indeed an extensive 
search (carried out by two of us”) for such solu- 
tions has shown that the widest of them contain 
only seven physically different arbitrary func- 
tions, i.e., one less than it is required for a 
general solution; hence even this solution in 
spite of its wideness is only a special case. In 
other words, this solution is unstable; there 
exist small perturbations which lead to its dis- 
sipation. Since in the synchronous reference 
system the singularity cannot disappear entire- 
ly, this means that it must go over, as a result 
of the perturbation, into a fictitious one. 

Thus we are led to the fundamental conclusion 
that the existence of a physical time singularity 
is not an obligatory property of the cosmological 
models of the general relativity theory. The 
general case of an arbitrary distribution of mat- 
ter and gravitational field leads to an absence 
of such a singularity. 

This result is formally equally valid for the 
singularities towards both directions of time. 
However, physically these directions are of 
course not equivalent and there is an essential 
difference between both cases already in the 
statement of the problem itself. The singularity 
in the future can have a physical meaning only 
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if it is admitted by quite arbitrary conditions 
given at any previous moment of time. On the 
other hand, it is perfectly clear that there are 
no reasons at all for the distribution of matter 
and field attained in the course of the evolution 
of the universe to comply with the specific con- 
ditions which are necessary for realization of 
the special solution with a physical singularity. 
Even if one admits the realization of such a 
specific distribution at some moment of time, 
it will inevitably be violated in the following time 
already as a result of the unavoidable fluctua- 
tions. ‘Therefore the above results exclude the 
possibility of the existence of a singularity in 
the future; this means that the contraction of 
the universe (if it is at all to come) must after- 
wards change again to an expansion. As to the 
singularity in the past, an investigation based 


only on the gravitational equations, can only 
impose certain restrictions on the admissible 
character of the initial conditions, the complete 
elucidation of which is impossible in the frame- 
work of the existing theory. 

A detailed account of this work will be pub- 
lished in the Journal of Experimental and Theo- 
retical Physics (U.S.S.R.). 

Our sincere thanks are due to Professor L. D. 
Landau for his constant interest in our work and 
for numerous discussions. 





‘A. Komar, Phys. Rev. 104, 544 (1956). 

2E. M. Lifshitz and I. M. Khalatnikov, J. Exptl. 
Theoret. Phys. (U.S.S.R.) 39, 149 and 800 (1960) 
(translations: Soviet Phys.—JETP 12, 108 (1961) and 
to be published]. 





DISPERSION RELATIONS FOR PRODUCTION AMPLITUDES 


Y. S. Kim 
Palmer Physical Laboratory, Princeton University, Princeton, New Jersey 
(Received February 27, 1961) 


Dispersion relations for production amplitudes 
were first proposed by Polkinghorne.’ His work 
was extended later by Kibble.” Logunov and other 
authors*~° have published a series of papers in 
which they obtain single-variable dispersion re- 
lations for various production processes. They 
presented a rigorous proof for the double Comp- 
ton effect.° On the other hand, for the process 
of pion production in pion-nucleon collision, their 
discussion is based on a rather delicate theorem 
on analyticity in the pion mass. In this note we 
shall demonstrate, by means of a concrete ex- 
ample from perturbation theory, that the con- 
jectured dispersion relation is in fact not valid, 
i.e., that the production amplitude has complex 
singularities when regarded as a function of the 
Logunov variable. 

For the process 1+N~-1+7+WN, we take the 
five independent variables in the form: 


E=-k-(p+p’), 

n= (k'-k")-(p+p')/(e +k”)-(p +p’), (1) 
x,=-k'-(p-p’), 

x,=-k"-(p-p’), 

v=(p-p’'f, 


where p and p’, respectively, denote the four- 
momenta of the incoming and outgoing nucleons; 
k refers to the incoming pion; k’ and k” refer to 
the two outgoing pions. 


p?=p?=-m?, k® =k? =k" =- p?; (2) 


m and yp represent, respectively, nucleonic and 
pionic masses. 

Logunov and his co-workers define the follow- 
ing two four-vectors in the Breit system where 
p+p’=0:4 


A=4[a7 (1 - tk’ +a'7(1+ Ek”), 
B=} [a7 '*k’- "|, 
where 


a =k,'/k,”, =-p?(a~! - w)/(2B"), 


and choose A,, A?, Bp, a, and p’ as the five 
independent variables. Then they investigate 
analyticity in A, and assert a cut-plane repre- 
sentation for the pion production amplitude for 
fixed physical values of the other four variables. 
We shall show in the following that a simple anal- 
ysis in perturbation theory leads to the conclu- 
sion that the existence of complex singularities 
is inevitable. 
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One can show, after some complicated algebra, 


that these five variables are related to the covari- 


ant forms E, n, x,, 2, andv by 


a =(1+n)/(1- n), 


v2 1+7\ 


--1(1- 1-n) 2_if/i+7 
A? =-i(1 (Ze Ge ea (4) 
-4(1- &)[ x, +x,+30- $y"), 


where 


-2u?n 
8-30? - (P(x, +ax,+0-3p)’ 





and 
1- 12 
auni( dat)". 


It is seen in the expressions (4) that the vari- 
ables A?, B-p, a, and p are independent of E 
and that Logunov and his co-workers, in effect, 
investigate analyticity in E for fixed n, *x,, x2, 
and v. 

Among many interesting Feynman diagrams, 
those with simple loops can be treated by the 
elegant method developed along with proofs of 
the Mandelstam representation.* We shall con- 
sider, in particular, three diagrams of Figs. 1 
and 2. These diagrams are characterized by 
triangular and square loops with unstable verti- 
ces. 

A standard analysis’ shows that coincident 
singularities of the triangle diagram lead to two 
complex branch points located at 


E =+([-3(y? +v) 
-i{ (4m? - up?) (u? +2x-v)(2x-v-3 u?)/4 7} J, 


(5) 








P Pp Pp Pr’ 


FIG. 1. Triangular diagrams due to the direct pion- 
pion interaction. 


314 


FIG. 2. Square diagram due to the direct pion-pion 
interaction. 


where 
X=U+X,+Xp- 


This conclusion is consistent with the Feynman 
“je” prescription. The quantity inside of the 
square root sign is positive, that is, the singu- 
lar points are complex if 


-(k’+k”P > 4y?. (6) 


The quantity (k’+k”)? depends only on the fixed 
variables, and the inequality (6) is always satis- 
fied for every set of physical n, x,, x,, andv. 

The square diagram also introduces the branch 
points of (5), corresponding to two end-point 
singularities. In addition to these, there are 
four real singularity points due to the leading 
coincident singularities. They need not be men- 
tioned here. 

The preceding discussion establishes the 
counter-example to the assertion of reference 3. 
Regarding the double Compton effect, y+N- 
y+y+N, we may mention that in addition to the 
dispersion relation of reference 5, other single- 
variable dispersion relations can also be proved. 

The results will be published elsewhere. 

The author would like to thank Professor S. B. 
Treiman for his suggestion and guidance, Pro- 
fessor M. L. Goldberger for his continued inter- 
est and encouragement, and Dr. S. Brown, Dr. 
L. F. Cook, and Dr. J. Tarski for many fruitful 
discussions. 
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In this section are printed the abstracts of Articles 
that have been forwarded to The American Institute of 
Physics for publication in THE PHYSICAL REVIEW. 
In quoting information obtained from this section be- 
fore the appearance of the corresponding Article, ref- 
erence should be made to “Physical Review (to be 
published)” rather than to this Journal. 


THEORY OF STRONGLY COUPLED MANY- 
FERMION SYSTEMS. I. CONVERGENCE OF 
LINKED CLUSTER EXPANSIONS. Leon N. 
Cooper, Brown University, Providence, Rhode 
Island (Received December 27, 1960). 


A strongly coupled system —the limiting case 
of a highly degenerate many-fermion system for 
which the variation of the kinetic energy is neg- 
lected, and the interaction restricted to a region 
of momentum space neighboring the Fermi sur- 
face—has been analyzed in a manner not de- 
pendent upon agsumptions about the convergence 
of power series expansions or on partial sum- 
mations of infinite series. The vacuum expecta- 
tion value of the resolvent operator, (1/(H -z)),, 
is expressed as the Laplace transform of the 
exponential of a function linearly dependent on 
the volume of the system. It is shown that the 
linked-cluster expansion of the vacuum expecta- 
tion value of the resolvent operator has a zero 
radius of convergence as a power series in the 
coupling constant. The most serious physical 
consequence of this is that a nontrivial inter- 
action never results in a “normal” system. 


GENERALIZED HARTREE-FOCK METHOD. 

J. G. Valatin, Department of Mathematical Phys- 
ics, University of Birmingham, Birmingham, 
England and Institute for Advanced Study, Prince- 
ton, New Jersey (Received November 21, 1960). 


A variational principle is formulated to de- 
termine the single-particle states, their pairing, 
and the occupation number distribution for a 
trial state vector of the Bardeen, Cooper, and 
Schrieffer type. The equations which are derived 
generalize those of the Hartree-Fock method ob- 
tained with a Slater-determinant trial wave func- 
tion. It is pointed out that in a suitable repre- 
sentation the vacuum state of a general quasi- 
particle transformation has such a trial form 
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which exhibits directly the pairing of single- 
particle states. Another variational principle 
determines the excitation energies. 

Two coupling cases are distinguished: the 
commutative case in which the self-consistent 
densities and energies are related to quantities 
which all commute, and the more general non- 
commutative case. The latter is of importance 
in critical-field phenomena. The equations for 
the commutative case can be written in a matrix 
form which retains its validity in the more gen- 
eral noncommutative case. The simple matrix 
commutator equations appear as direct general- 
izations of the density matrix form of the Hartree- 
Fock equations. 

The equations for small oscillations have an 
equally simple form. Their connection with a 
diagonal representation of the quasi-particle 
energies is exhibited in a way which remains 
valid in the general coupling case. The “unphys- 
ical” solutions are excluded by the supplementary 
condition. The contact with the Green’s function 
approach is established. The generalized matrix 
form of the Green’s function equations shows 
especially clearly the symmetry properties of 
the method. 


THEORY OF NEGATIVE IONS IN LIQUID HELI- 
UM. C. G. Kuper,* Department of Physics, Uni- 
versity of Illinois, Urbana, Illinois (Received 
January 11, 1961). 


It is shown that Atkins’ electrostriction model 
gives reasonable values for the mobility of posi- 
tive ions in liquid helium. However, that model 
cannot account for the observed difference be- 
tween positive- and negative-ion mobilities. 
Arguments are advanced in support of the “bub- 
ble” model; the negative ion is believed to be a 
free electron in a cavity of radius about 12 A. 
The bubble model leads to a mobility in fair 
agreement with experiment. 


*On leave of absence from Department of Theoreti- 
cal Physics, The University, St. Andrews, Scotland. 


LINKED- PAIR EXPANSIONS IN QUANTUM 
STATISTICS. Franz Mohling,* Department of 
Physics, Columbia University, New York, New 
York (Received July 21, 1960). 


In the quantum statistical method of Lee and 
Yang, the cluster functions of quantum statistics 
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are expressed in terms of the cluster functions 
of Boltzmann statistics, which in turn are com- 
puted in terms of certain two-body functions. In 
the present paper, following a detailed study of 
the Boltzmann cluster functions, it is shown that 
the symmetric representation can be used for 
the two-body functions and that large classes of 
diagrams can be summed. This leads to the in- 
troduction of linked-pair graphs to describe the 
functions of quantum statistics. The two-body 
functions are expressed in terms of the reaction 
matrix, and are therefore well-defined for hard- 
core repulsions. For weak potentials the method 
is shown to be equivalent to the theory of Bloch 
and DeDominicis. 


*Now at Cornell University, Ithaca, New York. 


VERY LOW TEMPERATURE FERMIGAS. Franz 
Mohling,* Department of Physics, Columbia Uni- 
versity, New York, New York (Received July 21, 
1960). 


The momentum distribution in a low-tempera- 
ture Fermi gas is investigated using the methods 
of quantum statistics developed by Lee and Yang 
together with the linked-pair expansions of the 
previous paper. It is shown that in order to de- 


termine the momentum distribution at very low 
temperatures two coupled integral equations 
must be considered, one in momentum variables 
and due to Lee and Yang, and the other in tem- 
perature variables. It is also shown that the 
dominant low-temperature behavior of the mo- 
mentum distribution can be extracted in terms 
of a certain function v’(k). For a low-density 
Fermi gas with strong, short-range, two-body 
interactions, it is shown to third order in the 
scattering parameters of the interaction that at 
T=0 the function v’(k) is equal to the free- 
particle momentum distribution. Also, the en- 
ergy and other thermodynamic quantities are 
expressed in terms of v’(k), so that the theory 
permits a generalization of perturbation-theo- 
retic results to nonzero temperatures. The 
ground-state energy, momentum distribution, 
and thermodynamic potential are calculated to 
third order in the scattering parameters. 


*Now at Cornell University, Ithaca, New York. 


CYCLOTRON EMISSION FROM PLASMAS WITH 
NON-MAXWELLIAN DISTRIBUTIONS. G. Bekefi, 
Jay L. Hirshfield, and Sanborn C. Brown, De- 


partment of Physics and Research Laboratory of 
Electronics, Massachusetts Institute of Technol- 
ogy, Cambridge, Massachusetts (Received Jan- 
uary 3, 1961). 


Cyclotron emission from high-energy plasmas 
is calculated for two classes of electron distri- 
bution functions: (a) those that decrease monot- 
onically with increasing electron energy, and (b) 
those that have one or more maxima displaced 
from zero energy. In (a) the emission does not 
differ greatly compared with the emission from 
a Maxwellian plasma of the same mean energy. 
In (b) the emission can grow exponentially with 
distance traversed in the plasma, resulting in 
a greatly enhanced loss of radiant power. 


RADIATION FROM FAST PARTICLES MOVING 
THROUGH MAGNETIC MATERIALS. T. B. Day, 
Physics Department, University of Maryland, 
College Park, Maryland (Received December 23, 
1960). 


The problem of the generation of a changing 
magnetic field due to the interaction of a fast 
particle with a magnetic medium is studied. 

This combined Cerenkov-spin wave effect is 
shown to give rise to a “ringing” of the spin sys- 
tem under certain conditions of frequency and 
angle of observation, at least within an approx- 
imate evaluation of the general Green’s function 
for the problem. Some striking differences from 
the usual Cerenkov effect are discussed and pos- 
sibilities of using this effect as a neutral magnetic 
moment detector or as a probe of magnetic ma- 
terials are mentioned briefly. 


POINT-DEFECT MIGRATION AND BINDING 
IN METALS. A. Sosin, Atomics International, 
Division of North American Aviation, Canoga 
Park, California (Received January 3, 1961). 


The kinetics of decay of an excess defect con- 
centration in metals is examined with special 
attention to the initial stages of decay. A particu- 
lar case, excess vacancy migration to sinks in 
a slightly impure metal, is treated in detail; 
analog computer plots of isothermal and con- 
stant-tempering-rate recovery studies are pre- 
sented and analyzed. Initial recovery is deter- 
mined by the migration energy only; final re- 
covery is determined by an energy generally 
less than the sum of the migration plus vacancy- 
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impurity binding energy but more than the migra- 
tion energy alone. Initial and final recovery are 
easily resolved into two annealing stages. The 
intermediate recovery range may actually give 
rise to a resistivity increase. 


EFFECT OF PLASTIC DEFORMATION ON THE 
LOW-TEMPERATURE IONIC CONDUCTIVITY 
OF POTASSIUM CHLORIDE. B. S. H. Royce and 
R. Smoluchowski, Princeton University, Prince- 
ton, New Jersey (Received January 10, 1961). 


Measurements were made on the ionic con- 
ductivity of KCl in the temperature range be- 
tween 25° and 250°C. Both “as-received” crys- 
tals and specimens having up to 24% plastic 
deformation were examined and the activation 
energies for ionic motion determined. It was 
found that up to 4% plastic deformation produced 
no change in the activation energy over the “as- 
received” specimens, whereas deformations 
above 10% produce a decrease of about 0.2 ev. 
After such deformation, structure normally 
present in the ln(IT) vs T~ plot is absent and 
the decrease in activation energy can be attri- 
buted to the generation of ion vacancies not 
associated with divalent impurities. Measure- 
ments on a deformed and annealed specimen 
support this view. 


NUCLEAR MAGNETIC RESONANCE SATURA- 
TION IN NaCl AND CaF,. Walter I. Goldburg, 
The Pennsylvania State University, University 
Park, Pennsylvania (Received January 3, 1961). 


Experiments were performed to test the 
hypothesis, due to Redfield, that a nuclear spin 
system in a sufficiently large rf field, H,(v), is 
properly described by a spin temperature re- 
ferred to a frame of reference rotating about 
the Zeeman field with the frequency v of the rf 
field. Measurements were made on the Na™ 
spins in NaCl and the F’® spins in CaF,. A com- 
bination of steady-state and pulse techniques 
was used to measure the magnetization M, as 
a function of the frequency and amplitude of 
H,(v). When H, is sufficiently large, the data 
show that Redfield’s theory is correct within 
ten percent for NaCl but appreciably in error 
for CaF,. The same type of measurement per- 
formed at low amplitudes of H, shows that the 
theory of Bloembergen, Purcell, and Pound is 
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in agreement with the NaCl measurements but 
measurably in error for CaF,. It is suggested 
that neither theory is correct for CaF, because 
the fluorine spins relax by means of paramag- 
netic impurities, and therefore do not relax 
independently with characteristic time T,, as 
is assumed in both theories. In an intermediate 
range of saturating field amplitudes neither 
theory is expected to apply; this was found to 
be the case experimentally. 


STUDY OF SUPERCONDUCTORS BY ELECTRON 
TUNNELING. Ivar Giaever and Karl Megerle, 
General Electric Research Laboratory, Sche- 
nectady, New York (Received January 3, 1961). 


If a small potential difference is applied be- 
tween two metals separated by a thin insulating 
film, a current will flow due to the quantum 
mechanical tunnel effect. For both metals in 
the normal state the current-voltage character- 
istic is linear, for one of the metals in the super- 
conducting state the current voltage character- 
istic becomes nonlinear, and for both metals 
in the superconductive state even a negative re- 
sistance region is obtained. From these changes 
in the current-voltage characteristics, the change 
in the electron density of states when a metal 
goes from its normal to its superconductive state 
can be inferred. By using this technique we 
have found the energy gap in metal films 1000 A- 
3000 A thick at 1°K to be: 2€ pp =(2.68 + 0.06)x 10° 
ev, 2€gy =(1.11+ 0.03) x10 ev, 2e;,, =(1.05 + 0.03) 
x107% ev, 2€4) =(0.32 + 0.03) x10™° ev. 

The variation of the gap width with tempera- 
ture is found to agree closely with the Bardeen- 
Cooper-Schrieffer theory. Furthermore, the 
energy gap in these films has been found to de- 
pend upon the applied magnetic field, decreasing 
with increasing field. 


SCATTERING OF ELECTRONS BY PHONONS 
AND IMPURITIES IN SEMICONDUCTORS. 

H. Reiss and A. I. Anderman, Research Depart- 
ment, Atomics International, Division of North 
American Aviation, Canoga Park, California 
(Received January 3, 1961). 


A theory is developed for the mobility of an 
electron in an n-type semiconductor under the 
combined scattering of phonons and impurities. 
An attempt is made to combine the two processes 
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in a rigorous manner by treating the effects of 
impurities as though due to an applied external 
field. Although in principle a solution for any 
value of tT, the phonon scattering relaxation 
time, should be possible, this paper only ar- 
rives at a limiting law for the mobility at small 
values of 7. In the region of applicability satis- 
factory agreement with experiment is achieved. 


PARAMAGNETIC RESONANCE OF Cr** IN 
YTTRIUM OXIDE. J. W. Carson, D. P. Devor, 
and R. H. Hoskins, Hughes Research Labora- 
tories, Malibu, California (Received January 9, 
1961). 


The paramagnetic resonance spectrum of Cr** 
in single-crystal yttrium oxide has been ob- 
served at microwave frequencies from 9 kMc/ 
sec to 71 kMc/sec at temperatures from 4.2°K 
to 300°K. The spectrum is described in terms of 
the spin Hamiltonian = g8H-S+D[S,? -4S(S + 1)], 
where g=1.97+0.01, S=3/2, and 2D=72.9+0.2 
kMc/sec at 300°K, and 72.7+ 0.2 kMc/sec at 
77°K and 4.2°K. 


FERROMAGNETISM IN DILUTE SOLUTIONS 
OF COBALT IN PALLADIUM. R. M. Bozorth, 
P. A. Wolff, D. D. Davis, V. B. Compton, and 
J. H. Wernick, Bell Telephone Laboratories, 
Murray Hill, New Jersey (Received January 9, 
1961). 


Measurements of the magnetic properties of 
solid solutions of cobalt in palladium show that 
ferromagnetism exists in solutions as dilute as 
0.1 atomic percent of cobalt, when the cobalt 
atoms are about 10 atomic diameters apart. 

The ferromagnetic Curie point is established 

in this alloy by three different methods as 7°K. 
Bohr magneton numbers are determined for the 
whole series of alloys, and the moment associ- 
ated with one cobalt atom is found to increase 
with dilution from 1.7 Bohr units in pure cobalt 
to a limit of 9 to 10 units in the greatest dilution. 
Lattice constants measured at room temperature 
show the effect of the onset of ferromagnetism 

at about 10 atomic percent cobalt. A simple 
model explains the magnitude and variation of 
the cobalt moment. It is also suggested that the 
exchange of spin waves may be responsible for 
the long-range interaction that causes ferro- 
magnetism in the dilute solutions. 


MECHANISMS OF DOUBLE RESONANCE IN 
SOLIDS. J. Lambe, N. Laurance, E. C. 
McIrvine, and R. W. Terhune, Scientific Lab- 
oratory, Ford Motor Company, Dearborn, 
Michigan (Received November 10, 1960). 


A study of electron-nuclear double resonance 
(ENDOR) in ruby and other solids demonstrates 
the existence of the “distant-ENDOR” effect, 
which involves a change in the EPR signal caused 
by the depolarization of “distant” nuclei (nuclei 
having negligible hyperfine interaction with the 
paramagnetic centers). In order to obtain inter- 
pretable data on the mechanism it proved nec- 
essary to perform most of the experiments 
without modulation, observing not the deriva- 
tives but the functions y’ and y’’ themselves, the 
dispersive and absorptive parts of the spin sus- 
ceptibility. The former shows a large decrease 
upon application of rf power at a nuclear transi- 
tion frequency; the latter shows a moderate 
increase. Both the distant ENDOR (Al’’ nuclear 
Zeeman frequencies) and local ENDOR (Cr* 
hyperfine frequencies) affect the EPR with a 
response time comparable to the spin-lattice 
relaxation time of the distant aluminum nuclei. 
Nuclear-nuclear double-resonance experiments 
show that applied rf corresponding to Cr®™ nu- 
clear transitions depolarizes Al*’ nuclei. Both 
of these observations are consistent with a 
mechanism involving dynamic nuclear polariza- 
tion. A theoretical analysis of this mechanism, 
based on forbidden transitions involving distant 
nuclei, gives good agreement with observed 
nuclear polarizations and with the observed be- 
havior of x’, but predicts small increases in x’’. 
The increased absorption signal may be explained 
by enhanced spectral spin diffusion or by a spin 
packet considerably wider than assumed. Dis- 
tant ENDOR is expected to occur quite generally. 


NUCLEAR MAGNETIC SPECIFIC HEAT IN TWO 
FERROMAGNETIC IRON ALLOYS. C. T. Wei, 
Michigan State University, East Lansing, Michi- 
gan, and C. H. Cheng and Paul A. Beck, Uni- 
versity of Illinois, Urbana, Illinois (Received 
January 3, 1961). 


From the nuclear magnetic specific heat, 
measured at 1.6° to 4.2°K, Hog, at the Co nuclei 
in Co, ,Fe,,, was calculated to be 312 koe, while 
H oe at the V nuclei in V, ,,Fe,.,, is 61 koe, or 
less. Both of these alloys are bcc and ferro- 
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magnetic. The large difference in the H g¢r 
values may be associated with the fact that in 
Co, 3Fe,,, the Co™® nucleus is located in an atom 
with appreciably polarized 3d electrons, while 
in V,.ssF €o.67 the V™ nucleus is the only abundant 
nuclide with a nuclear magnetic moment and the 
atomic moment of V is very small or zero. 
Since in ferromagnetic alloys the polarization 
of the core s electrons is expected to be much 
stronger in those atoms which do have polarized 
d electrons than in adjacent atoms which do not, 
the above results suggest that, in the alloys 
investigated, the dominant contribution to H g¢¢ 
arises through Fermi contact interaction from 
the polarization of the core s electrons, as 
found for iron by Hanna et al. 


OPTICAL CONSTANTS OF GERMANIUM IN THE 
REGION 0-27 ev. Om P. Rustgi,* J. S. Nodvik, 
and G. L. Weissler, Department of Physics, 
University of Southern California, Los Angeles, 
California (Received October 21, 1960; revised 
manuscript received February 13, 1961). 


The reflectivity, R, has been measured for Ge 
crystal over the energy range 7.6 to 18.0 ev and 
the transmissivity, T, for an evaporated Ge 
film over the range 16.6 to 27.4 ev. The reflec- 
tivity data near normal incidence, combined with 
existing data at lower energies, were used in the 
Kramers-Krénig dispersion relation to evaluate 
the index of refraction, m, and extinction coeffi- 
cient, k, of Ge over the range 0 to 18 ev. Elec- 
tron characteristic losses are predicted at 10 
and 15 ev, consistent with those observed at 11 
and 16 ev. The plasma frequency of bulk Ge as 
determined by the optical data is Rw =10 ev. 

*Present address: Smithsonian Astrophysical Ob- 


servatory, 60 Garden Street, Cambridge 38, Massa- 
chusetts. 


MAGNETIC RESONANCE IN CANTED FERRI- 
MAGNETS. Perry Miles, Laboratory for Insula- 
tion Research, Massachusetts Institute of Tech- 
nology, Cambridge, Massachusetts (Received 
October 19, 1960; revised manuscript received 
February 13, 1961). 


The classical theory of the uniform (k=0) modes 
of a four-sublattice canted ferrimagnet is de- 
veloped. Two of these modes should lie in the 
microwave range for reasonable values of ap- 
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plied field, anisotropy, and exchange constants. 
One is the normal low-frequency mode familiar 
in collinear ferrimagnets; the other is excited 
by longitudinal rf fields at a frequency depending 
on anisotropy and angle of cant. Observation of 
this mode should allow analysis of ferrimagnetic 
structures and phase changes. 


LOW-TEMPERATURE THERMAL RESISTANCE 
OF n-TYPE GERMANIUM. Robert W. Keyes, 
International Business Machines Research Labo- 
ratory, Poughkeepsie, New York (Received 
November 21, 1960). 


It is proposed that the scattering of phonons by 
donors in germanium at low temperatures re- 
sults from the large effect of strain on the en- 
ergy of an electron in a hydrogenlike donor state. 
A calculation of the thermal conductivity with 
this scattering mechanism is presented. Reason- 
able agreement with the following features of the 
observed thermal conductivity is obtained: the 
very large scattering power of donors, the dif- 
ference between the scattering powers of anti- 
mony and arsenic, a temperature dependence of 
thermal conductivity stronger than T°, anda 
dependence of the scattering on number of occu- 
pied donors rather than on the total impurity 
concentration. 


ENERGY, SPECIFIC HEAT, AND MAGNETIC 
PROPERTIES OF THE LOW-DENSITY ELEC- 
TRON GAS. W. J. Carr, Jr., Solid State Physics 
Department, Westinghouse Electric Corporation, 
Pittsburgh, Pennsylvania (Received April 14, 
1960; revised manuscript received July 18, 1960). 


A perturbation expansion in powers of 7, = 
has been used to investigate the ground-state en- 
ergy of a dilute electron gas, the results being, 
in rydberg units per particle, E =-1.792r,~* 
+2.667,°%* + br, * + O(r,-™) +terms falling off 
exponentially with r,“*. The dimensionless 
parameter 7, is the radius of the unit sphere in 
Bohr radii. The term in r,™ is the energy of a 
body-centered cubic lattice of electrons as cal- 
culated by Fuchs; the r,~™” term is the zero- 
point vibrational energy of the lattice, for which 
the normal modes have been calculated, the re- 
sult differing only slightly from the value given 
by Wigner; and br,~ is the first-order effect 
of anharmonicities in the vibration. The con- 
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stant b has been estimated, its magnitude being 
small compared with unity. 

The vibrational part of the specific heat has 
been calculated, and a first-order approxima- 
tion has been obtained for the exponential terms 
inthe energy. Part of this energy comes from 
exchange, which leads to the result that, except 
for very low densities (ry, 2270), the electron 
spins are antiferromagnetically aligned. An 
order of magnitude for the Néel temperature 
has been calculated. 


EFFECT OF PRESSURE ON THE ABSORPTION 
EDGES OF SOME III-V, II-VI, AND I-VII COM- 
POUNDS. A. L. Edwards and H. G. Drickamer, 
Department of Chemistry and Chemical Engineer - 
ing, University of Illinois, Urbana, Illinois (Re- 
ceived December 13, 1960). 


The effect of pressure to 160 kilobars was 
measured on the absorption edges of the III-V 
compounds AlSb, GaSb, InP, and InAs, the 
II-VI compounds CdS, CdSe, and CdTe, and the 
I-VII compounds CuCl, CuBr, and Cul. It is 
possible to discuss the band structure of the 


II-V compounds with reasonable assurance 
relative to known group IV and III-V structures. 
For the II-VI and I-VII compounds, ionic and 
other effects must be important. A number of 
new phase transitions were noted at high pres- 
sure. For CuCl and CuBr the T-P curves of 
some of these transitions were established. 


EFFECT OF FISSION SPECTRUM NEUTRONS 
ON x- TYPE GERMANIUM. Daniel Binder, Nu- 
clear Measurements Department, Nuclear Elec- 
tronics Laboratory, Hughes Aircraft Company, 
Culver City, California (Received January 13, 
1961). 


The electron removal rate for n-type germa- 
nium irradiated with fission spectrum neutrons 
is 8+ 1 per neutron at room temperature. This 
value is compared with the results of monoen- 
ergetic neutron irradiations from 2 to 5 Mev. 
The fact that the removal rate is roughly con- 
stant is explained by the constancy of the energy 
dissipated in elastic collisions. 


NUCLEAR SPIN-SPIN INTERACTION ENERGY 
IN THE HYDROGEN MOLECULE. J. P. Auffray 


and J. W. Cooley, Institute of Mathematical 
Sciences, New York University, New York, New 
York (Received January 5, 1961). 


An accurate theoretical estimate of the inter- 
action energy of the two proton magnetic mo- 
ments in the v=0, J=1 vibrational-rotational 
level of the electronic ground state of H, is ob- 
tained. Agreement with the experimental value 
for the nuclear spin-spin interaction energy is 
found to be within 1 part in 10°. This is of the 
order of magnitude of the experimental error. 


FREQUENCY SHIFTS IN HYPERFINE SPLIT- 
TING OF ALKALIS: A CORRECTION. R. Herman 
and H. Margenau, Yale University, New Haven, 
Connecticut (Received December 27, 1960). 


The effect of the deformation of the wave func- 
tions by Van der Waals interactions, previously 
ignored, upon the hyperfine shifts of alkalis caused 
by rare gas atoms is computed. It is found to be 
large and clearly in need of consideration. When 
applied to the experimentally observed shifts, the 
model proposed earlier, with the new values of 
the interaction constants, leads to “interaction 
radii” somewhat greater than before and more 
nearly equal to gas kinetic radii. 


STARK BROADENING OF HYDROGENIC ION 
LINES IN A PLASMA. H. R. Griem, University 
of Maryland, College Park, Maryland and U. S. 
Naval Research Laboratory, Washington, D. C., 
and K. Y. Shen,* University of Maryland, Col- 
lege Park, Maryland (Received December 5, 
1960). 


The frequency distributions of the ionized he- 
lium lines at 4686 A and 3203 A broadened by 
the local fields of both ions and electrons in a 
plasma are calculated in the classical path ap- 
proximation, which is shown to be always ap- 
plicable. General formulas are developed for 
the Stark profiles of lines from multiply ionized 
hydrogenic systems, and the validity domains 
of the impact and quasi-static approximation for 
electron and ion broadening are delineated. The 
results are compared with the Holtsmark theory 
and an approximation for high series members. 


*Now at the Theoretical Institute of Physics, 
Naples, Italy. 
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ALPHA-PARTICLE IONIZATION IN POLY- 
ATOMIC GASES AND THE ENERGY DEPEND- 
ENCE OF W. William P. Jesse, St. Procopius 
College, Lisle, Illinois (Received January 4, 
1961). 


An extended series of measurements by three 
different methods has been carried out to de- 
termine the variation of W, the average energy 
to make an ion pair in the gases N, and C,H, as 
a function of the energy of the ionizing alpha 
particle. In one method the ionization ratios 
were determined in the two gases for single 
alpha particles from two collimated polonium 
sources, the particles from one source being 
reduced in energy by passage through a suc- 
cession of interchangeable mica windows. The 
corresponding energy ratios were determined 
by auxiliary measurements in pure argon, the 
W values in argon being assumed constant. W 
values for alphas of initial energy up to 9 Mev 
were also obtained by a comparison in C,H, and 
N, of the relative ionization from Po and ThC’ 
single alpha particles. Data from all these ex- 
periments indicate a continuous decrease in W 
values in C,H, and N, with increasing alpha 
energy over a range from 1 to 9 Mev. The dif- 
ferential w (defined as the ratio of energy in- 
crement AE to ionization increment A/ at any 
point on the alpha path) also decreases similarly 
and seems to approach but never quite reach, 
within the limits of alpha energy so far investi- 
gated, the corresponding W value for beta parti- 
cles. No explanation of this continuous decrease 
with energy of the W values can at present be 
advanced. 


DISSOCIATION OF MOLECULAR IONS BY 
ELECTRIC AND MAGNETIC FIELDS. John R. 
Hiskes, Lawrence Radiation Laboratory, Uni- 
versity of California, Berkeley, California 
(Received December 5, 1960). 


A general discussion of the dissociation of 
diatomic molecules and molecular ions by elec- 
tric fields is presented. These calculations 
pertain primarily to the ground electronic states 
of the molecular systems. The H,* ion is treated 
in considerable detail; the required fields for 
the dissociation range from 10° v/cm for the 
uppermost vibrational state to 210° v/cm for 
the ground state. The many-electron homo- 
nuclear ions are treated in successive charge 
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states. The HD*, HT*, HD, LiH*, and LiH** 
heteronuclear ions are considered. The disso- 
ciation of homonuclear ions and heteronuclear 
ions exhibit distinctly different features. The 
HD* and HT* ions are more susceptible to dis- 
sociation than is H,*. The extent to which the 
dissociation by an electrostatic field and by the 
Lorentz force, ev xB, are equivalent is con- 
sidered. The rates of induced dipole transitions 
to lower vibrational states can be made negligibly 
small compared with the dissociation rates. The 
application of this work to particle accelerators 
and to the injection problem for fusion devices 
is discussed. 


SIMULTANEOUS EFFECT OF DOPPLER AND 
FOREIGN GAS BROADENING ON SPECTRAL 
LINES. Louis Galatry,* Sloane Physics Labora- 
tory, Yale University, New Haven, Connecticut 
(Received December 29, 1960). 


Using the classical Fourier integral theory, an 
expression is given for the shape of a spectral 
line, broadened by phase changes due to colli- 
sions and by the actual changes in velocity of the 
emitting particles resulting from collisions. The 
result is not a simple Voigt-type folding of an 
exponential into a dispersion distribution; it ex- 
hibits the contraction noted by Dicke and leads 
to the usual formulas when the time interval 
between path-deflecting or phase-disturbing col- 
lisions becomes very great. 


*Present address: Centre National de la Recherche 
Scientifique, Bellevue (Seine et Oise), France. 


COHERENCE EFFECTS IN RESONANCE FLUO- 
RESCENCE. M. E. Rose, Oak Ridge National 
Laboratory, Oak Ridge, Tennessee, and R. L. 
Carovillano, Boston College, Chestnut Hill, 
Massachusetts (Received December 13, 1960). 


The crossed-level method of atomic spectro- 
scopy is discussed and the angular distribution 
formula for both the incoherent and coherent reso- 
nance scattering is derived. The form of this 
distribution function, as given here, explicitly 
displays the geometric factors depending on radia- 
tion propagation vectors. With the application to 
hydrogen in mind, the distribution function is ex- 
pressed explicitly for single electron transitions 
with external fields possessing axial symmetry. 
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The properties of the distribution function are 
discussed with emphasis on the case of unpolar- 
ized radiation. For the case of hydrogen there 
are two possible applications of major interest. 
The first concerns the possibility of a precision 
measurement of the 2p fine-structure splitting 
and, hence, a determination of the fine-structure 
constant. Explicit results for the shape and other 
properties of the resonance line with a uniform 
magnetic field are obtained. The other applica- 
tion is concerned with the possibility of meas- 
uring the 2s-2p Lamb splitting. This requires 

an electric field parallel to the magnetic field. 
Unfortunately, the level crossings which are 
sensitive to the Lamb splitting cannot radiate 
sufficiently rapidly while those which do radiate 
appreciably occur at field strengths which are 
extremely insensitive to the Lamb splitting. 


COULOMB EXCITATION OF THE SECOND 

2+ STATES IN W, Os, AND Pt NUCLEI. F. K. 
McGowan and P. H. Stelson, Oak Ridge National 
Laboratory, Oak Ridge, Tennessee (Received 
January 5, 1961). 


The location of a second 2+ state has been 
established for six even-even nuclei by means 
of Coulomb excitation produced by 4- to 5-Mev 
protons. The relatively weak excitation of these 
states is detected by a measurement of the 
gamma-ray yields from singles spectra and 
from coincident measurements of the cascade 
gamma rays. The B(E2)’s for decay of the 
second 2+ state to ground state by the cross- 
over transition exhibit some uniformity for the 
even-even isotopes of W and Os, being about 6 
times the single-particle value. The cascade/ 
crossover ratio for the decay of the second 2+ 
state is known for these nuclei. The upper cas- 
cade B(E2)’s exhibit enhancements of 10 to 60 
times the single-particle value. The ratios of 
the B(£2)’s for decay of the first and second 2+ 
states are compared to the predictions of several 
collective models. For five of these nuclei the 
E2/M1 ratio is known for the upper cascade 
transition. The B(M1) values obtained are ex- 
ceedingly small compared to the single-particle 
estimate. This result is in qualitative agree- 
ment with the collective models which predict 
that M1 radiation is forbidden in the decay of 
vibrational excitations. 


DECAY OF Er'” AND Tm’”. C. J. Orth and 
B. J. Dropesky, Los Alamos Scientific Labora- 
tory, University of California, Los Alamos, 
New Mexico (Received January 6, 1961). 


Er’” has been produced by double neutron 
capture in enriched Er’”. The beta decay of 
50.4-hr Er’” and its daughter, 63.7-hr Tm’”, 
has been studied with a solenoidal beta spectrom- 
eter and beta and gamma scintillation spectrom- 
eters. The highest energy group of the Tm!” 
beta spectrum has an end-point energy of 1.83 
Mev; this group represents the Tm'” + Yb!” 
ground-state beta transition. The beta decay 
of Tm’” is accompanied by gamma rays of the 
following energies: 0.079, 0.180, 0.422, 0.495, 
0.915, 1.095, 1.29, 1.39, 1.41, 1.47, 1.51, and 
1.59 Mev. A decay scheme for Tm?” is pro- 
posed with excited states in Yb’” at 0.079, 0.259, 
1.174, 1.47, 1.55, 1.59, and 1.67 Mev. The beta 
decay of Er'” is accompanied by gamma rays 
of the following energies: 0.050 (Tm K x ray), 
0.108, 0.126, 0.408, and 0.610 Mev, the last 
of which represents a transition to the Tm?” 
ground state. The beta spectrum measured in 
coincidence with the 0.610-Mev gamma ray has 
an end-point energy of ~0.26 Mev, which estab- 
lishes a decay energy of 0.87 Mev for Er’”. 


(d,p) REACTIONS IN DEFORMED HEAVY NU- 
CLEI. G. B. Holm,* J. R. Burwell,f and D. W. 
Miller, Department of Physics, Indiana Univer- 
sity, Bloomington, Indiana (Received January 
6, 1961). 


The spectra of protons from (d,p) nuclear re- 
actions initiated by 11- Mev deuterons incident 
on Th?™, U*55, y*55, U25*, and Pu?*® have been 
obtained using separated isotopic targets. The 
protons were analyzed in a double-focusing mag- 
netic spectrometer. States with assignments 
known from radioactive-decay studies are identi- 
fied and assignments for some other observed 
states are discussed using the Nilsson model 
and Satchler’s theory of (d,p) reactions in de- 
formed nuclei. From an analysis of the results 
the following binding energies of the last neutron 
are obtained: U*™, 6.834 0.11 Mev; U***, 4.74 
+ 0.06 Mev; and Pu™°, 6.494 0.05 Mev. 


*On leave from the Nobel Institute of Physics, 
Stockholm, Sweden. 

TNow at the University of Oklahoma, Norman, 
Oklahoma. 





VoLUME 6, NUMBER 6 


PHYSICAL REVIEW LETTERS 


Marcu 15, 1961 





EXCITATION FUNCTIONS FOR THE (a, an) AND 
(a, 2pn) REACTIONS ON Ce’. B. M. Foreman, 
Jr.,* Chemistry Department, Brookhaven Na- 
tional Laboratory, Upton, New York (Received 
January 11, 1961). 


Cross sections have been measured radio- 
chemically for the reactions Ce™**(a, an)Ce™ 
and Ce***(a@, 2pn)Ce** in the helium-ion energy 
range 16.8-40.1 Mev. The cross section for the 
(a, an) reaction begins to rise sharply at about 
25 Mev and reaches a value of 69+ 5 mb at 40.1 
Mev. The cross section for the (a, 2pn) reaction 
begins to rise at about 32 Mev and reaches a 
value of 2.5+ 0.4 mb at 40.1 Mev. An upper limit 
of ~0.1 mb for the cross section for the reaction 
Ce’**(a, ap)La™ in the energy range covered by 
this study has also been obtained. The results 
for the (a, an) and (a, ap) reactions are dis- 
cussed in terms of these possible mechanisms: 
compound nucleus formation and decay, knock- 
on, and direct inelastic scattering followed by 
neutron evaporation. The results seem to be 
most consistent with the last mechanism. The 
existence of a measurable cross section for the 
(a, 2pn) reaction in this energy region suggests 
that the reaction proceeds mainly by He*® emis- 
sion, probably by a stripping mechanism. The 
data reported are consistent with the hypothesis 
that in this energy range at most one particle 
is emitted as a result of direct interaction. 


*Present address: Department of Chemistry, 
Columbia University, New York, New York. 


REDUCED WIDTHS AND ISOTOPIC SPIN IM- 
PURITIES OF }+ STATES OF N’°. J. B. French,* 
Sarah Mellon Scaife Radiation Laboratory, Uni- 
versity of Pittsburgh, Pittsburgh, Pennsylvania, 
and Fysisch Laboratorium der Rijksuniversiteit, 
Utrecht, The Netherlands, Syurei Iwao,f Uni- 
versity of Rochester, Rochester, New York, 

and Erich Vogt, Atomic Energy of Canada 
Limited, Chalk River, Ontario, Canada (Re- 
ceived January 3, 1961). 


The nuclear reactions C™( p, n)N** and 
C**( p, y.)N"*, for protons of less than 1.6 Mev, 
involve the interference of two adjacent states 
having the same spin and parity ($+) but differ- 
ent isotopic spin (7 =1/2 and 3/2). By taking 
properly into account the effect of other, more 
distant, $+ levels on the cross section near the 
interfering pair of levels, we are able to fit 
well the (p,m) cross section from the neutron 
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threshold up to a proton energy of 1.6 Mev and 
to obtain reliable estimates for the reduced- 
width amplitudes of the interfering pair, as well 
as for the physically significant phases of the 
amplitudes. Since the neutron decay of the 

T =3/2 state is “forbidden,” the neutron re- 
duced widths of the close-lying pair lead toa 
direct measure, 4%, of the isotopic spin im- 
purity of the pair. The results found in the 
cross-section analysis are compared to shell- 
model calculations based on the N** wave func- 
tions of Halbert and French, and reasonably 
satisfactory agreement is found. 


*On leave from the University of Rochester, Roches- 
ter, New York. 

tPresent address: Physics Department, Univer- 
sity of Syracuse, Syracuse, New York. 


SYMMETRIC AND ASYMMETRIC FISSION. 
Henry W. Newson, Duke University, Durham, 
North Carolina (Received January 9, 1961). 


Fission yields have been calculated assuming 
as a first approximation that they are proportion- 
al to the product of the level densities of a pair 
of binary fission products. The level densities 
have been calculated with the simplified shell- 
model methods of Newson and Duncan. The 
calculations predict a single symmetric peak 
when the mass of the fissioning nucleus, A,<N, 
+Z,+N,+Z,=50+28 +82 +50=210 (in agreement 
with the observed fission yields for bismuth and 
lighter elements), and three maxima in the fis- 
sion yield curves for heavier compound nuclei. 
The peak corresponding to approximately equal- 
size binary fission products is very much higher 
than is observed experimentally. This is un- 
doubtedly due to the fact that in asymmetric 
fission a core corresponding to 82 neutrons and 
50 protons remains intact in the heavier fission 
product, whereas for symmetric fission this 
core is disrupted at the cost of several Mev. 
Since correction for this energy effect involves 
a number of unknown factors, the calculated 
yields for symmetric fission have been reduced 
by the same empirical factor in all calculations. 
An additional parameter, n, is introduced in 
correcting for excitation energy of the fission 
products and for possible departures from equi- 
librium. These calculations, which involve only 
two free parameters, explain most of the fission 
yield data for all five known cases where the 
compound nucleus is within a Mev or so of the 
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fission threshold: Pu*®, U?%*, 755, y*53, Th?5?, 
but it is necessary to treat nm as a free param- 
eter for each curve to fit the small steep regions 
on each side of mass number A,/2. The calcu- 
lated fission yields of the more highly excited 
compound nucleus, Ac”’, predict three equally 
prominent maxima in qualitative agreement with 
observation. 


DECAY OF Zn®. J. B. Cumming and N. T. 
Porile, Chemistry Department, Brookhaven 
National Laboratory, Upton, New York (Re- 
ceived January 11, 1961). 


The decay of 38.4-min Zn® has been investi- 
gated using scintillation and beta spectrometer 
techniques. 84% of the decays are to the Cu® 
ground state. Gamma rays having energies of 
0.67, 0.96, and 1.42 Mev are present with in- 
tensities of 9.0, 6.7, and 0.9% of the g* transi- 
tions. Low-intensity gamma rays (<0.2%) were 
observed at energies of 1.55, 1.83, 2.04, 2.34, 
2.55, 2.77, and 3.10 Mev. Coincidence meas- 
urements established positron feeding of the 
Cu® levels at 0.67, 0.96, and 1.42 Mev but no 
y-y coincidences were observed. Internal con- 
version coefficients of the 0.67- and 0.96-Mev 
gamma rays are in agreement with predominant- 
ly M1 assignments to both transitions. A decay 
scheme is presented which differs significantly 
from that previously reported. It is inferred 
that the spin of Zn® is 3/2~. M1 and E2 transi- 
tion probabilities between the various levels of 
Cu® and Cu® are discussed in terms of the 
“center-of-gravity” model for states in these 
nuclei. Both agreements and disagreements 
with the model predictions are observed. 


SYMMETRY OF NEUTRON-INDUCED U** FIS- 
SION AT INDIVIDUAL RESONANCES. G. A. 
Cowan, Anthony Turkevich, * C. I. Browne, and 
LASL Radiochemistry Group, Los Alamos 
Scientific Laboratory, University of California, 
Los Alamos, New Mexico (Received December 5, 
1960). 


Neutrons in the resonance energy region from 
a nuclear explosion were resolved by time-of- 
flight and used to induce fissions in U*** attached 
to a revolving wheel. The symmetry of fission 
at individual resonances from approximately 10 
to 60 ev was examined by radiochemical means. 


As measured by the ratio Ag" /Mo™, the prob- 
ability of symmetric fission decreased at some 
resonances by a maximum of 10% compared to 
thermal fission of U*** and at other resonances 
increased by a maximum of 40%. With varying 
degrees of assurance, nine resonances are 
identified with an increase in symmetry; five 
more regions of increased symmetry are asso- 
ciated with resonances or a background effect. 
Twenty resonances are identified with a decrease 
in symmetry. 

In a sample containing 500 levels in the reso- 
nance region, there was no level with a Ag™™ 
yield even one-thirtieth as great as Mo”. It is 
thus very improbable that there are any neutron 
resonances in U*** that lead to predominantly 
symmetric fission. 


*Permanent address: Enrico Fermi Institute for 
Nuclear Studies, University of Chicago, Chicago, 
Illinois. 


PSEUDOSCALAR INTERACTION AS THE ORIGIN 
OF ANOMALOUS BETA-DECAY SPECTRUM 
SHAPES. J. M. Pearson, Université de Montréal, 
Montréal, Canada and Western Reserve Univer- 
sity, Cleveland, Ohio (Received December 19, 
1960). 


The (1+5/W) anomaly factors which have been 
found by Langer and co-workers in the spectrum 
shapes of all the Gamow-Teller allowed and 
unique forbidden transitions that they have meas- 
ured are here interpreted in terms of the ex- 
istence of a pseudoscalar interaction in nuclear 
beta decay. Following previous treatments of 
the pseudoscalar interaction, we consider the 
gradients of both the lepton covariant and the 
nuclear potential. It then follows that a pseudo- 
scalar interaction will distort the shapes of 
these spectra more than those of the AJ=0 (yes) 
decays, which have been frequently, but incon- 
clusively, examined in this connection. The 
shape factors for the transitions in question are 
derived and found to display the observed low- 
energy anomaly. In particular, it is found that 
the fact that b has the same sign for both elec- 
tron and positron emitters can be accounted for. 
The single nuclear matrix element ratio that is 
involved depends only on the nuclear coupling 
scheme and in particular is independent of the 
nuclear radial wave functions. Confining atten- 
tion to those nuclides for which pure j-j coupling 
can be assumed, the required coupling constant 
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ratio is gp/g,4 =-180. In addition to producing 
low-energy distortions of spectral shapes, it is 
found that a pseudoscalar interaction should be- 
come effective again at very high energies, 
W>aZ/R. Such effects may be expected in the 
B’, N’*, and N** decays, unless the momentum 
transfer dependence of gp is such as to lead at 
high energies to values of &p/&, considerably 
smaller than the aforementioned low-energy 
value. 


GAMMA-GAMMA DIRECTIONAL CORRELA- 
TIONS IN Nd'*’. Atam P. Arya,* Department of 
Physics, The Pennsylvania State University, 
University Park, Pennsylvania (Received May 6, 
1960; revised manuscript received August 1, 
1960). 


Directional correlation measurements have 
been made on the 320-92 kev and 280-320 kev 
gamma-ray cascades in Pm**’ following the 
decay of 11.1-day Nd’*’ with a coincidence 
scintillation spectrometer using Nal detectors. 
The observed correlation functions are W(é@) 
=1 - (0.1030 + 0.0298)P,(cosé) + (0.0107 + 0.0099) 
xP,(cosé@), and W(@)=1+(0.0710+ 0.0162)P,(cosé) 
- (0.0126 + 0.0103)P,(cosé@), respectively, for the 
two cascades. The spin and parity assignments 
for energy levels of Pm‘*’ at ground state, 92 
kev, 410 kev, and 690 kev were found to be 3*, 
3*, #*, and §*, respectively. It was found that 
the 92-kev gamma ray has a mixture of (95+ 2)% 
M1 and (5+2)% E2 with 6,, =+0.229+ 0.143, the 
320-kev gamma ray has a mixture of 1% M1 and 
99% E2 with 6,..=+9.95+ 0.11, and the 280-kev 
gamma ray has a mixture of 99% M1 and 1% E2 
with 6,,,=-0.11+0.11. 


*Permanent home address: 3623 Chawri Bazar, 
Delhi-6, India. 


SPIN AND HYPERFINE STRUCTURE OF ARSE- 
NIC-76. R. L. Christensen,* D. R. Hamilton, 

H. G. Bennewitz,t J. B. Reynolds,} and H. H. 
Stroke, ll Palmer Physical Laboratory, Princeton 
University, Princeton, New Jersey (Received 
January 13, 1961). 


Hyperfine structure in the *S,, ground state of 
the radioactive atom As” has been investigated 
by the method of magnetic resonance in an atomic 
beam produced by microwave discharge dissocia- 
tion of arsenic vapor. AF =0 resonances were 


326 


observed within both the F =5/2 and F=7/2 
atomic levels at several values of magnetic field 
up to about 5 oersteds, indicating that the spin of 
the As”* nucleus is 2. An analysis of multiple 
quantum transition spectra within the same F 
states gave a measurement for two of the hfs 
intervals: Avyg 2 =+(117+4) Mc/sec and Arg, x» 
=+ (69+ 16) Mc/sec, with the same sign for both. 
From the value of the hfs constant A, the magni- 
tude of the magnetic field at the arsenic nucleus 
is (1.33+ 0.15) x10° oersteds, in reasonable agree- 
ment with the variation in this field among similar 
atoms. The value ofg,; has been found to be 1.994 
+ 0.003 for arsenic. 


*Now at IBM Research Center, Yorktown, New York. 

tPermanent address: University of Bonn, Bonn, 
Germany. 

t Now at the Department of Physics, Washington Uni- 
versity, St. Louis, Missouri. 

|| Now at the Department of Physics, Massachusetts 
Institute of Technology, Cambridge, Massachusetts. 


FAST NEUTRON ACTIVATION CROSS SECTION 
OF Au’. S. A. Cox, Argonne National Labora- 
tory, Argonne, Illinois (Received November 3, 
1960). 


The neutron activation cross section of gold 
was measured in the neutron energy range from 
30 kev to 1500 kev. The absolute value of the 
cross section was based on the U*** fast-fission 
cross section which was used for absolute neutron 
flux measurements from 200 kev to 1500 kev. 
For measurements below 200 kev, the B’%(n, a) 
cross section was used for monitoring the neu- 
tron flux. The relative cross section from 30 
kev to 200 kev was then normalized at 200 kev 
to the absolute measurement. The results agree 
well with recent measurements except for some 
pulsed beam-—liquid scintillator measurements 
and spherical shell transmission measurements 
which yield much lower cross-section values. 


NUCLEAR ENERGY LEVELS OF Na™ IN THE 
REGION FROM 350 TO 630 kev. Carl T. Hibdon, 
Argonne National Laboratory, Argonne, Illinois 
(Received May 9, 1960). 


The neutron cross-section data from 350 to 
630 kev show 71 peaks, consisting of a relatively 
small number of large peaks and many small 
peaks. Each of the previously known large peaks 





VoLUME 6, NUMBER 6 


PHYSICAL REVIEW LETTERS 


Marcu 15, 1961 





was resolved into two or more components. The 
analyses show a few s-wave levels, a small 
number of p-wave levels, and a large number 

of d- and f-wave levels. For all of the levels 

of Na™ up to 630 kev, a plot of the number of 
levels having energies <E, as a function of the 
neutron energy E, shows an essentially linear 
distribution. The distribution of the angular 
momenta is in agreement with the theoretical 
distribution for a value of o=1.8. The level 
spacings appear to agree with an exponential 
distribution. For the reduced neutron widths, 
the results appear to agree equally well with the 
exponential and Porter-Thomas distributions. 
The strength function obtained from the reduced 
widths has an average value of 0.045 for both 
values of J for /=0 and an average value of 0.37 
for all values of J for /=1. For higher values 
of /, the strength function is too large. An ex- 
pression developed for the distribution of the 
levels above the ground state tends to agree with 
the data for a value of 0.50 Mev for 6, the aver- 
age level spacing of the nucleons in the nucleus. 


ELASTIC AND INELASTIC SCATTERING OF 
ALPHA PARTICLES BY N"*. W. D. Ploughe, * 


Purdue University, Lafayette, Indiana (Received 
January 12, 1961). 


The absolute differential cross section for the 
elastic and inelastic scattering of 19.2-Mev 
alpha particles by N** was measured as a func- 
tion of the scattering angle. The forward part 
of the ground-state angular distribution was 
fitted with the black-nucleus diffraction model 
of Blair using an interaction radius of 5.89 x10-* 
cm. The angular distribution for the Q = -3.95 
Mev group was fitted with [j,(qR)]* using an 
interaction radius of R =5.9x107 cm. With 
the same interaction radius, a best fit to the 
data for the unresolved doublet, Q = -4.91 and 
-5.10 Mev, was obtained using the sum [j,(qR)]? 
+[j,(q¢R)]?. The best fit by odd-order spherical 
Bessel functions is consistent with the assign- 
ment of negative parity to both the 4.91- and 
5.10-Mev levels. 


*Now at Ohio State University, Columbus, Ohio. 


COMPARISON OF THE SCATTERING OF POSI- 
TRONS AND ELECTRONS FROM NUCLEAR 
CHARGE DISTRIBUTIONS. George H. Rawitscher, 


Yale University, New Haven, Connecticut, and 
C. Rutherford Fischer, New Mexico State Uni- 
versity, Las Cruces, New Mexico (Received 
January 3, 1961) 


Elastic-scattering cross sections of 183-Mev 
positrons and electrons are calculated for vari- 
ous charge distributions of the Ca and Au nuclei. 
It is shown that the combined use of positron 
and electron scattering measurements can lead 
to a determination of the nuclear charge distri- 
bution which is more accurate than that derived 
from either one of the scattering cross sections 
when used by itself. The scattered particles 
obey Dirac’s equation and the nuclei are as- 
sumed to be static spherically symmetric charge 
distributions, whose radial dependence is given 
in terms of a three-parameter family of curves. 


1*-PROTON SCATTERING AT 990 Mev. J. K. 
Kopp, Brookhaven National Laboratory, Upton, 
New York, A. M. Shapiro, Brown University, 
Providence, Rhode Island, and A. R. Erwin, 
University of Wisconsin, Madison, Wisconsin 
(Received January 13, 1961). 


Positive-pion scattering at 990+30 Mev has 
been examined in a 6x3 x2 in. hydrogen bubble 
chamber without a magnetic field. The elastic 
and inelastic scattering cross sections were 
found to be 15.3+ 1.5 mb and 12.6+3.3 mb, re- 
spectively. The inelastic cross section includes 
0.1949:°8 mb of 2+ -K* production and 0.78 + 0.14 
mb of 1*2*n~p production. A simple pion-pion 
model which predicts the branching ratios for 
double-pion production in 1~-p collisions is 
found to be inconsistent with the double-pion 
production observed in this experiment. The 
relation of the experiment to 1~-p experiments 
in the region of the second and third resonances 
is discussed. 


PHOTOPRODUCTION OF PIONS IN CARBON. 
T. R. Palfrey, B. M. K. Nefkens,* L. Mortara, 
and F. J. Loeffler, Department of Physics, 
Purdue University, Lafayette, Indiana (Re- 
ceived December 23, 1960). 


Positive and negative pions were produced by 
photons on a carbon target and observed at 
laboratory angles of 35°, 73°, and 121°. At each 
angle the yield of mesons of constant energy 
was observed by a magnetic spectrometer as 
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a function of peak bremsstrahlung energy. Seven 
values of the latter, ranging from 205 to 335 
Mev, were used. Yields and 1~/n* ratios, cor- 
rected for various systematic experimental 
errors, are presented. By using a photon dif- 
ference method, the bremsstrahlung spectra 
were unfolded from the yield curves to give 
meson cross sections versus photon energy at 
fixed pion energies. These functions are com- 
pared with predicted yields which consider the 
internal momentum distribution of the target 
nucleons; the agreement is adequate. 


*Fulbright Travel Grant recipient from Ryksuni- 
versiteit, Utrecht, The Netherlands. 


LONG-RANGE INTERACTION IN K-NUCLEON 
AND K-NUCLEON ELASTIC AMPLITUDES. 
Fabio Ferrari,* Graham Frye,! and Modesto 
Pusterla,} Lawrence Radiation Laboratory, Uni- 
versity of California, Berkeley, California (Re- 
ceived November 4, 1960). 


A method of calculating, for the K-nucleon in- 
teraction, the long-range force arising from the 
exchange of a pion pair and of a possible three- 
pion resonant state is formulated. It is shown 
that the long-range force can be related with the 
electromagnetic structure parameters of the 
nucleon and K meson. Finally, relations between 
K-nucleon and K-nucleon elastic amplitudes are 
discussed. 


*Now at Istituto di Fisica, Universita di Padova, 
Padova, Italy. 
TNow at Physics Department, University of Wash- 


— Seattle, Washington. 
Now at Istituto di Fisica Teorica, Universita di 


Napoli, Napoli, Italy. 


ENERGY DEPENDENCE OF THE LOW-ENERGY 


K~-PROTON AND K*- PROTON CROSS SECTIONS. 


Fabio Ferrari,* Graham Frye,! and Modesto 
Pusterla,} Lawrence Radiation Laboratory, Uni- 
versity of California, Berkeley, California (Re- 
ceived November 4, 1960). 


The K -proton and K*-proton S-wave scatter- 
ing is analyzed by using a relativistic effective- 
range formula derived by studying the analytic 
properties of partial-wave scattering amplitudes. 
The influence of the pion-pion interaction on the 
elastic and reaction cross sections is discussed. 


*Now at Istituto di Fisica, Universita di Padova, 
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Padova, Italy. 

TNow at the Physics Department, University of 
Washington, Seattle, Washington. 

tNow at Istituto di Fisica Teorica, Universita di 
Napoli, Napoli, Italy. 


EMISSION OF PHOTONEUTRINOS AND PAIR 
ANNIHILATION NEUTRINOS FROM STARS. 
Hong-Yee Chiu, Institute for Advanced Study, 
Princeton, New Jersey, and Robert C. Stabler,* 
Laboratory of Nuclear Studies, Cornell Univer- 
sity, Ithaca, New York and Aeronutronic Divi- 
sion, Ford Motor Company, Newport Beach, 
California (Received December 28, 1960). 


Field-theoretic calculations of the cross sec- 
tions for the photoneutrino process, y+e~ +e~ 
+v+v, and the pair annihilation process, e+ +e~ 
+vy+v, are performed in order to obtain the 
neutrino luminosity of very hot stars (T,25x10 
°K). The energy loss rate which is obtained for 
the latter process is sufficient to determine the 
rate of evolution of the stellar core when T; 
=5x10° °K. 

*Now at Aeronutronic Division, Ford Motor Company, 
Newport Beach, California. 


POLARIZATION OF THE RECOIL PROTON 
FROM 7° PHOTOPRODUCTION IN HYDROGEN. 
J. O. Maloy, California Institute of Technology, 
Pasadena, California, G. A. Salandin, Istituto 
di Fisica dell’Universita, Padova, Italy and 
Istituto Nazionale di Fisica Nucleare, Sezione 
di Padova, Italy, A. Manfredini and V. Z. Peter- 
son, Istituto di Fisica dell’ Universita, Roma, 
Italy and Istituto Nazionale di F'sica Nucleare, 
Sezione di Roma, Italy, and J. L Friedman* and 
H. Kendall, Stanford University, Palo Alto, Cali- 
fornia (Received January 9, 1961). 


The D,, nature of the second resonance in neu- 
tral single pion photoproduction, y+p—-p+n°, 
suggested by Peierls, has been confirmed by 
additional experimental observations of the polari- 
zation of the recoil proton over a range of photon 
energies. The photon energy dependence of the 
polarization at 90° c.m. is in substantial disagree- 
ment with alternative models suggested by Stop- 
pini and Pellegrini, and Landovitz and Marshall 
if the observed angular distributions are also 
considered. An experimental method using nu- 
clear emulsion as scatterer-detector, in con- 
junction with a magnetic spectrometer, is shown 
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to have both good energy resolution and reason- 
able counting rate. 


*Present address: Massachusetts Institute of Tech- 
nology, Cambridge, Massachusetts. 


INTERNAL STATE OF A GRAVITATING GAS. 

G. E. Tauber and J. W. Weinberg, Department 
of Physics, Western Reserve University, Cleve- 
land, Ohio (Received January 15, 1958; revised 
manuscript received December 28, 1959). 


The significance of a theory of gravitational 
equilibrium of concentrated masses is discussed 
in connection with possible general relativistic 
effects in white dwarf stars. The covariant form 
of phase space and Liouville’s theorem is de- 
veloped, using the canonical equations for a 
particle under gravitational and electromagnetic 
forces. The dynamical isotropy of the ideal 
fluid is formulated, and the associated equations 
of state and allowed streaming patterns are 
found. A covariant kinetic theory yields general 
relativistic forms for the Maxwell and Fermi 


distributions in the case of thermal equilibrium, 
and limits their streaming to rigid motion. 
Rotating fluids are studied in comoving coordi- 
nates, and the problem of determining their 
gravitational equilibrium is reduced, in most 
cases of physical interest, to a simple standard 
form with constant density and vorticity. 


KINEMATICAL AND DYNAMICAL RESONANCES. 
A. O. Barut and K. H. Ruei,* Department of Phys- 
ics, Syracuse University, Syracuse, New York 
(Received January 16, 1961). 


A method is given to distinguish between the 
solutions of the dispersion relations correspond- 
ing to kinematical and dynamical resonances. It 
consists of studying the resonance energy as a 
function of the coupling constant. The method is 
illustrated for potential scattering, for charged 
scalar meson theory, and for resonances due to 
unstable particles. 


*Now at Essex College, Windsor, Ontario, Canada. 





